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Atom Smashing Versus Humanity Smashing 


R. Arthur H. Compton was reported re- 
cently as declaring that the progress made 
in breaking up atoms and putting a portion of 
their energy to other uses is of far more impor- 
tance to the world today than the war in Europe. 
To those who are familiar with the rich and con- 
tinuing rewards of research in all fields of science 
such a reminder is scarcely needed. Yet it must 
be repeated again and again in the daily press if 
the population as a whole is ever to recognize 
that promoting research is just as much a part 
of its responsibility as promoting big armies and 
big navies for defense purposes. 
The need for emphasizing the importance 
of research was never 


ous that, if a nation expands its geographical fron- 
tiers, it can do so only at the expense of other 
people. How different the situation is if its intel- 
lectual frontiers are expanded. There is no con- 
ceivable limit to the expansion which is possible. 
No one will lose anything and all will gain. 

In his book, Atoms in Action, Dr. G. R. Harri- 
son estimates that, if we could discover a way 
to use the energy of the sun in an efficient 
manner, we could increase the wealth (in terms 
of available energy which is the only sensible 
unit) of every person on earth by a factor of 
200,000 times. Some way must be used to bring 
to the attention of the so-called ‘‘man in the 

street’? such examples 


greater than at the pres- 
ent time. We have be- 


of the potentialities of 
research. Strong state- 


fore us the strange spec- 
tacle of nations which 
have themselves reaped 
countless benefits from 
science forsaking the 
methods of research for 
those of bloodshed and 
destruction. Their pur- 
pose is said to be to up- 
hold national honor and 
\o give their people a 
better living. Yet with 
the very finite area of 
‘his world on which we 
live it is perfectly obvi- 
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ments coming from 
leading scientists like 
Dr. Compton will un- 
doubtedly do a great 
deal of good. If the bene- 
fits which can accrue 
from scientific research 
become generally known, 
there will be no chance 
for the rise of dictato- 
rial governments bent 
upon the destruction of 
humanity for the sake 
of a small geograph- 
ical expansion. 
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Ferromagnetic Domains and the Magnetization Curve 


By WILLIAM FULLER BROWN, JR. 


Palmer Physical Laboratory, Princeton University, Princeton, New Jersey 


The magnetization curve and hysteresis loop of a ferromagnetic 
material used to be explained qualitatively through the inter- 
action of “‘molecular magnets.”’ It is now known that inter- 
action of the atoms produces, directly, only magnetic saturation, 
and that the “‘magnets” whose interaction is responsible for the 
details of the hysteresis curve are much larger things—mag- 
netically saturated regions big enough to be seen under a 
microscope. These “spontaneously” 
known as “‘“domains,”” 


magnetized regions are 
and with them progress is being made 
toward a quantitative theory of the magnetization process. The 


basic ideas of the domain theory are described in this article. 


Introduction 

HE magnetization curve of a typical ferro- 

magnetic specimen, Fig. 1, may for theo- 
retical purposes be divided into three parts: 
saturation (IIT), reversible magnetization at high 
fields (II), and the low field range (1) in which 
hysteresis occurs. It is with Part I that we shall 
be particularly concerned in this article, but in 
order to arrive at an understanding of it some- 
thing must first be said about Parts II and ITI. 


Saturation (111) 


The large magnetic moments observed in 
ferromagnetic materials at saturation are at- 
tributed, in modern theory, to the “spins” of 
atomic electrons. The forces alining these spins 
Are ones whose origin can be explained only by 
quantum mechanics, but whose effect can be 
taken into account, for our purposes, by at- 
tributing to every pair of neighboring atoms a 
classical potential energy 


u=—}I cos ¢, (1) 
where J is a positive constant, and @ is the 
angle between the directions of the two magnetic 
moments: see Fig. 2. Since u is a minimum when 
@=0, the tendency of these forces is to produce 
perfect alinement of all the atomic moments. 
At the absolute zero of temperature there is 
nothing to oppose this tendency, and there 
results a magnetic moment per unit volume 


Jo=nyo, (2) 
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where is the number of atoms per cubic 
centimeter and yo is the magnetic moment per 
atom. 

Temperature agitation tends to break up this 
perfect alinement of atomic moments; the 
“spontaneous” magnetization J, therefore de- 
creases as the temperature rises, and _ finally 
vanishes at a critical temperature @, the Curie 
point. This occurs when the thermal agitation 
energy kT becomes of the same order of rmagni- 
tude as the energy per atom associated with the 
alinement forces; or 


(3) 


From the value of 6 (=631°K for nickel) and of 
Boltzmann's constant k erg per 
degree) it follows that J for nickel is of order of 
magnitude 10-" erg, and the corresponding 
energy per unit volume is of order of magnitude 
In=10-" XK 10 = 10" ergs per cm*. This is much 
greater than the energy density associated with 
interatomic magnetic forces, which is given in 
order of magnitude by 327/22 X (500)? =5 X10° 
ergs cm*. Magnetic forces, as is well known, are 
much too weak to account for ferromagnetism. 

At room temperature J, is only slightly smaller 
than the value Jo at absolute zero, and the 
problem may be simplified by supposing that the 
atomic moments are perfectly alined but have a 
slightly decreased value uw, such that 


J,= nu. (4) 
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The alinement of spins is accompanied by a 
distortion of the lattice (saturation magnetostric- 
tion). The magnetic forces exerted by the atoms 
on each other are capable of producing such a 
distortion. A calculation by classical theory, 
however, leads to a result that is correct only in 
its order of magnitude. No calculation by 
quantum mechanics has been carried out. The 
saturation magnetostriction must at present be 
evaluated from experimental data. 


Reversible Magnetization at High Fields (II 


The atomic theory of spontaneous magnetiza- 
tion predicts that the magnetization should have 
the saturation value J, even in zero applied field. 
The theory of Parts Land IT of the magnetization 
curve must therefore take account of factors 
that were neglected in the simple atomic theory. 

The most obvious of these is the crystalline 
nature of the material. Ordinary polycrystalline 
material is made up of crystal grains with their 
crystal axes oriented in different directions. 
Within such a crystal, we may expect the spon- 
taneous magnetization to orient itself along 
certain axes of symmetry, unless a sufficiently 
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iG. 1, Typical magnetization curve. I] = magnetic field; 
‘=magnetic moment per unit volume, or magnetization. 
‘IT, saturation; 17, slow and reversible magnetization; /, 
rapid and largely irreversible magnetization; O, demag- 
ictization. Arrows indicate the directions in which the 
curves are traversed. At points P and Q the magnetization 
s the same, but the field is different. The ‘‘reversible sus- 
eptibility” is the slope of PP’ or QQ’. 


strong field is present to pull the magnetic 
noments into the field direction. In Part II of 
‘he magnetization curve of a_ polycrystalline 
specimen, the magnetic moments of the indi- 
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vidual crystals are gradually turned out of 
directions determined by the crystalline forces, 
and into the direction of the field. The com- 
ponents of magnetization normal to the field 
direction cancel each other for the specimen as a 
whole, so that only the component in the field 
direction contributes to the observed mag- 
netization. 


Fic. 2. Energy u of the alining forces between the magnetic 
moments yo of atoms that are nearest neighbors. 


In single crystal specimens, the normal com- 
ponent should be observable; moreover, if the 
field is along a direction in which the magnetiza- 
tion tends to point because of crystalline forces, 
Part II of the magnetization curve should be 
absent, and saturation should be attained in 
much smaller fields. Both of these conclusions 
are verified experimentally: the normal com- 
ponent of magnetization is observed, and certain 
directions in the crystal are found to be “‘direc- 
tions of easy magnetization.” 

The theory that has been successful here is a 
formal one that says nothing about the nature 
of these atomic forces, but deduces formulas for 
their macroscopic effects from the requirements 
of thermodynamics and of crystalline symmetry. 
Such a theory involves a number of constants 
which cannot be determined from the theory 
but whieh can be evaluated from the experi- 
mental data. It is then the problem of atomic 
theory to explain the values of these constants. 
This portion of the atomic theory is still incom- 
plete, although the general nature of the atomic 
forces involved is understood. The formal theory 
is sufficient for our purposes, and we need 
consider only two very simple cases. 

The first is crystalline cobalt. At room tem- 
perature the direction of easy magnetization is 
the hexagonal axis, and the crystal is mag- 
netically isotropic about this axis. The magnetiza- 
tion vector therefore lies in the plane determined 
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if 
x 
 u=-gl cos ¢ 
a 
a 
. 
| 
: 


by the direction of easy magnetization and the 
field direction, as is shown in Fig. 3(a). It may be 
regarded as subject to a torque Lx, due to 
“anisotropy” forces, tending to pull it into the 
direction of easy magnetization, and to a torque 
Lu =H J, sin (¢— 6), due to the field, tending to 
pull it into the field direction. It is simpler, 
however, to consider not torques, but an energy 
function W(@), the thermodynamic “‘free energy” 
per unit volume, from which the torque Z acting 


(a) (b) Ly 
Ae 
Hexagonal axis Tension & 


snto 


Fic. 3. Two simple cases of magnetic anisotropy: (a) 
cobalt crystal; (b) isotropic material under tension (mag- 
netostriction positive). The spontaneous magnetization 
(J,) is pulled toward the direction of the field (7) by a 
torque Ly, and toward the direction of easy magnetization 
by a torque Lx or Lz. The torque Lx or Lg may be derived 
from an energy function or Wg. 


Wy = K' sin’'@ + 


on the magnetic moment of unit volume can be 
derived by means of the relation L= —dW. 08. 
The condition for stable equilibrium at a given 
IT is that @ shall have such a value as to make W 
a minimum. The part of W due to the external 
field is the potential energy of a magnet of 
moment J, oriented at an angle (¢—4@) to the field 
II, namely Wy = —JIJ, cos (@— 84). The part of W 
due to the crystalline forces must by symmetry 
be of the form 


Wr=K’ sin? K” (5) 


Good agreement with experimental magnetiza- 
tion curves is obtained by setting K’=5.1X 10°, 
K" =2.2X10° (both ergs/cm*), and neglecting 
higher powers of sin @. 

The second case, illustrated in Fig. 3(b), is 
that of an isotropic material subjected to uniform 
and constant tension ¢ dynes, cm? by the action 
of forces of external origin. In this case Wx is 
replaced by W,, the potential energy of the 
weight producing the tension. The saturation 
magnetostriction A, (fractional elongation in the 
direction of the spontaneous magnetization) will 
be assumed to be positive ; in this case the length 
of the specimen in the direction of the tension 
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is greatest when 6=0 or x. Thus the weight is 
lowest, and in the absence of a field the minimum 
of W is attained, for 6=0 or x—that is, the 
direction of the tension is a direction of easy 
magnetization. It can be shown that for any @, 
the potential energy per unit volume due to 
the tension is 

sin® 0. (6) 


Thus W, is of the same form as the first term 
in Wr. 


The Demagnetized State (O) 


Since the atomic theory requires the existence 
of spontaneous magnetization even in zero 
applied field, it might be supposed that in the 
demagnetized state the microcrystals are mag- 
netized to saturation, and that the zero large- 
scale magnetization observed in ordinary meas- 
urements results from the random orientation of 
the different crystals in the polycrystalline 
specimen. But the experiments on single crystals 
show that they, too, can be demagnetized, and in 
fact Part I of the magnetization curve is about 
the same for them as for ordinary specimens. 
The conclusion is that even within a crystal, the 
uniform magnetization required by the atomic 
theory extends only over limited regions or 
“domains,” and changes its direction from one 
direction of easy magnetization to another at 
boundaries between such domains. This involves 
less offense to the quantum and anisotropy forces 
than would a complete lack of alinement of the 
atomic magnets; for only the comparatively few 
atoms at the boundaries between domains will 
feel the difference between this and a completely 
uniform magnetization. 

According to this picture, Part I of the mag- 
netization curve begins (O, Fig. 1) with equal 
volumes magnetized along each of the directions 
of easy maghetization, and ends (X) with the 
entire crystal or microcrystal magnetized along 
the one of these directions that is nearest to the 
field direction. A further increase of field rotates 
the magnetization vector closer to the field 
direction, giving Part II of the curve, unless the 
angle between the field and the nearest direction 
of easy magnetization happens to be zero, in 
which case the crystal is already saturated at 
the end of Part I. 
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We shall consider later the process by which 
ihe erystal passes from an initial state of demag- 
netization to a final state in which only favored 
directions, of easy magnetization are occupied. 
First, we shall consider the following questions: 
(1) What direct experimental evidence is there 
that such domains exist, and what information 
do the experiments give us about their nature? 
(2) What theoretical reason is there why the 
spontaneous magnetization should extend only 
over such limited regions, rather than over the 
entire specimen? (3) Can theory explain the 
observed size and shape of the domains? 


(1) EXPERIMENTAL EVIDENCE FOR DOMAINS 


The domains are not merely abstractions 
required by theory: they can be both heard 
and seen. 

They can be heard by surrounding the speci- 
men with a test coil connected to an amplifier 
and phones, and gradually applying a magnetic 
field. Loud cracks are heard, indicating sudden, 
discontinuous changes of the magnetic flux 
through the test coil—as might be expected if 
small chunks of the specimen change their 
directions of magnetization one by one as the 
field is increased. This is the well-known Bark- 
hausen effect. 

They can be seen by a refinement of the 
familiar “iron filings’? experiment. A properly 
prepared magnetic powder or colloid, placed on 
a polished surface of a ferromagnetic specimen, 
forms patterns (visible in a microscope) that not 
only reveal the presence of magnetized domains 
in the macroscopically demagnetized material, 
but give considerable information abont their 
size, shape, and arrangement. 

The original crude method of observing the 
Barkhausen effect, by listening to noises made by 
ihe domains, can be replaced by a more sophisti- 
cated procedure in which the phones are re- 
placed by a photographic record and quantitative 
information is obtained. The information relates 
‘o the regions that change their directions of 
Magnetization discontinuously; these are not 
necessarily identical with the spontaneously mag- 
netized domains, since it may occur that only 
part of a domain is involved in such a discon- 
(inuous change. The conclusions reached by 
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(a) (b) (c) 


FG. 4. Magnetic colloid patterns obtained by Elmore on 
a cobalt surface approximately perpendicular to the basal 
plane with applied normal field & outward; (b) zero; 
(c) inward. The lines run parallel to the hexagonal axis. 
Magnification approximately 65 x. 


Bozorth and Dillinger from experiments of this 
sort are: (1) The effect occurs in single crystals 
as well as in polycrystalline material. (2) The 
size of the Barkhausen regions does not vary 
much with the material or with its mechanical 
or thermal history ; they may be larger or smaller 
than the crystal grains. (3) In wires 60 cm long 
and 1 mm in diameter, the average volume of a 
Barkhausen region varies from 1.2 for 
iron to 45 X10~-* cm’ for a 50 percent iron-nickel 
alloy. These experiments leave open the question, 
which had not been raised at the time they were 
performed, of a possible dependence of the 
domain size on the dimensions of the specimen. 

The surface patterns, discovered by Bitter, 
have been studied in great detail by Elmore. 
One of his photographs is shown in Fig. 4. This 
is a photograph of patterns obtained on a cobalt 
surface; the direction of easy magnetization is 
approximately in the plane of the figure, and the 
lines of the pattern run parallel to that direction. 
In the middle pattern, the average spacing of the 
lines is about 6X10-* cm. The spacing varies 
somewhat from one crystal grain to another; 
the lines are closer together in the smaller 
grains. 


(2) ORIGIN OF THE DOMAIN STRUCTURE 


The quantum forces described by Eq. (1) tend 
to produce a spontaneous magnetization extend- 
ing throughout the whole specimen. The aniso- 
tropy forces of Eq. (5) or (6) can orient this 
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magnetization in certain preferred directions, but 
have no tendency to orient different regions in 
different ones of these directions. 

With magnetic forces the situation is different. 
Magnetic atoms, regarded as small current 
circuits, tend to line up so as to link as much 
magnetic flux as possible. Fig. 5(a) shows the 
effect of this tendency in the case of two atoms. 


(a) 


Fic. 5. Influence of magnetic forces. (a) Field lines of an 
atomic magnet of moment u located at O. The arrow with 
head and tail points in the direction of the magnetic 
moment; the circular arrow shows the direction of the 
atomic currents producing it. A second atomic magnet at 
A or at B tends to set itself as shown. (b) Arrangement 
favored by quantum forces: a uniform magnetization of 
the specimen. (c) Arrangement favored by magnetic forces: 
long, narrow domains. (d) Arrangement opposed by mag- 
netic forces: short, wide domains. 


One is located at O; its field lines are shown. 
The second atom, if placed at A, alines itself in 
the same direction; but if placed at B, it alines 
itself in the opposite direction. In the case of a 
specimen containing many atoms, the uniform 
magnetization shown in (b) is less acceptable 
magnetically than the arrangement shown in (c), 
in which the flux is kept in closed paths inside 
the material rather than being allowed to 
escape to the outside. 

This rough argument may be made more 
rigorous by considering the field h, at a lattice 
point, produced by all the atoms at other lattice 
points. Stability requires that this be in the 
direction of the magnetic moment wp of the atom 
at the lattice point in question, and as large as 
possible, in order that the mutual magnetic 
energy of all the pairs of atoms, —}Zh-p, may 
be a minimum. The overwhelming majority of 
the atoms are far enough from any surface of 
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discontinuity so that this field is given by the 
Lorentz formula 


h=H+3}rJ, (7) 


where H is the macroscopic field and J is the 
vector magnetization in that vicinity. J has the 
same direction as w and a fixed magnitude J,, in 
(c) as well as in (b); the energy per unit volume 
contributed by this part of h is —jrJ#, a 
quantity already mentioned in the discussion of 
saturation. But H is oppositely directed in (b), 
being the demagnetizing field due to the poles 
at the ends of the specimen; while in (c), HZ is 
zero, since there are no poles. Thus (c), though 
not necessarily the best arrangement, is at least 
better than (b) as far as magnetic forces are 
concerned. On the other hand (d) would be 
worse than (b). 

The requirements of the quantum forces are 
not satisfied by arrangement (c) at the bound- 
aries between domains. Consider for simplicity a 
simple cubic lattice with lattice constant a, and 
suppose first that the boundary between oppo- 
sitely magnetized domains is sharp, as indicated 
in Fig. 6(a), where each arrow represents an 
atomic magnet. The mutual energy of neighbor- 
ing atoms on opposite sides of the boundary is 
increased, according to Eq. (1), from —J/2 to 
+I/2 by having the magnetic moments anti- 
parallel instead of parallel; there are 1/a? such 
pairs of atoms per unit area of the boundary, 
and therefore a sharp boundary involves addi- 
tional energy of amount J/a*? per unit area. 
A gradual transition (Fig. 6(b) and (c)) requires 
less energy. Here the angle 6 between the mag- 
netic moments of the atoms and the direction 
of easy magnetization is a slowly varying func- 
tion of x, and the angle ¢ between the magnetic 
moments of successive atoms is the small 
quantity ad@/dx. Thus their energy exceeds the 
value for parallel alinement by 4/(1—cos ¢) 
=1]¢?=}la*(d0/dx)*. If we now imagine the 
transition made more gradual by stretching out 
the curve in (c) over a greater thickness d, the 
number of pairs of atoms affected increases as d, 
but the excess energy of each pair decreases as 
1/d*, so that the net effect is a decrease as 1/d. 
Thus the quantum forces tend to make the 
transition as gradual as possible. In order of 
magnicude d6/dx=7/d, the number of pairs of 
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atoms per cm*® of boundary 
=4(1/a*)d, and hence the 
energy per unit area=4Ja’ 
=I /ad, a de- 
crease by a factor of about a/d 


as compared with the value 
for an abrupt transition. 

But now the anisotropy 
forces come into play. We 
suppose that the anisotropy 
energy per atom is of the form 
A sin® 6; this corresponds to 
the first term in Eq. (5) or to 
Iq. (6), with A equal to K’a* 
in the one case and to $),ca* 
in the other. Then the atoms 
in the transition layer have 
excess energy because their 


moments are turned out of the 
direction of easy magnetiza- 
tion. If the transition is made 
more gradual, the number of 
atoms affected increases as d, 
and since in this case there is 
no compensating decrease in 
the energy per atom, the total anisotropy energy 
per cm*® of boundary increases as the boundary 
thickness d. In order of magnitude this energy is 
A(sin® (1/a*)d=2A (d/a*). 

Thus the total energy associated with unit 
area of the interdomain boundary consists of a 
term y, proportional to 1/d, due to the quantum 
forces, and a term yq4 proportional to d, due to 
the anisotropy forces. The former tend to 
thicken the boundary layer and the latter to 
thin it. The actual thickness d will be that for 
which the sum of these terms is a minimum 
(see Fig. 6(d)). 

It is easily shown that the function 

y(d) =ad+B/d 
is a minimum when the two terms are equal; 
then d=(B/a)' and y=2(a8)!. Here /a* 
and BI] /a. Thus, except for numerical factors 
of order of magnitude unity, we have 


Actu. 


d=a(I/A)}, 


I 2a 
y (8) 
a d 


lf J=10-" and A = K’a*210-", we get d=100a, 
‘hat is, a transition layer about 100 atoms thick 
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Fic. 6. Factors determining thickness of interdomain boundaries (a—d) and of 
domains (e). (a) A sharp boundary: too large quantum energy. (b) A gradual 
transition. (c) Graphical representation of (b). (d) Dependence of boundary 
energy per unit area, y, on thickness of the boundary, d. y,, quantum energy; 
Ya, anisotropy energy; y, total energy, a minimum when y,.=7,. (e) Dependence 
of total boundary energy, W>, and end energy, W., on thickness of the domains, w. 


between domains. The total energy per cm? of 
boundary is less than the value for a sharp 
boundary in the ratio 2a/d, or 1/50. 

A more rigorous analysis verifies the results 
just obtained, and gives also the form of the 
function 0(x). 


(3) S1IzE AND SHAPE OF THE DOMAINS 


Returning to Fig. 5(c), we may now endeavor 
to determine the thickness w of the domains 
themselves. In the expectation that this is much 
greater than the thickness d of the boundaries 
between them, we may take the boundaries into 
account sufficiently well by treating them as 
surfaces of discontinuity with a surface energy 
density y given by (8). This energy by itself 
would tend to make the boundary area as small 
as possible, and thus to decrease the number of 
domains and increase their thickness. In fact 
there are c/w boundaries each of area bl, and 
the total boundary energy of the specimen is 
ybcl,'w, decreasing as w increases. 

Thick domains, however, involve large regions 
at the ends (see Fig. 6(e)), in which the mag- 
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‘ netization is turned squarely out of the direction 


of easy magnetization and the maximum aniso- 
tropy energy density, A/a’, is realized. The two 
triangular prisms at opposite ends of a boundary 
together form a rectangular prism of length 6 and 
square cross section (w/v2)*=w* 2; and since 
there are c/w such pairs of prisms, the total 
anisotropy energy is 
A bw* c 1A 


=— —bcw, 
2 w 2a’ 


increasing as w increases. 

Thus the boundary energy tends to make the 
domains thick, the end energy to make them 
thin. The actual thickness is determined by the 
minimum of the total energy, bc(,Aw a*+yl/w), 
which is the same type of function encountered 
before. The equilibrium value of w is found by 
equating the two terms, whence 


2a*yl\! f2a*l 2(1A)*} 
w=( ) -| | (9) 
A A a? 


According to this the size of the domains depends 
on the size of the specimen, a result that is to be 
expected if surface as well as volume energies are 
present. For specimens of linear dimensions 
=1 cm, with d=100a=10~* cm, we get w10-* 
cm, the order of magnitude observed in surface 
patterns. The volumes observed in the Bark- 
hausen effect, 10-* cm’, are consistent with this, 
if they are interpreted as chunks of domains 
whose linear dimensions are of the same order of 
magnitude as the thickness of the domains. 

We began by objecting to the large magnetic 
energy accompanying a uniform magnetization ; 
have we ended up with an energy enough 
smaller to justify our more complicated arrange- 
ment? Insertion of (9) in the formula for the total 
energy, and division by the volume bel, gives 
for the mean energy per unit volume 


2A 108 
W=— =——- ergs cm’*. (10) 
Vl 


The magnetic energy of a uniformly magnetized 
specimen is given approximately in terms of a 
demagnetization factor D by the formula 


Wuag= DJ 2210°D ergs ‘cm?, (11) 
For 17=100 cm, W< provided D>0.001. 
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This requires that the diameter of the specimen 
be greater than 0.4 cm. The theory therefore 
fails to give a satisfactory description of the 
domain structure of long thin specimens. 

In other respects, too, the analysis is incom- 
plete. It has been shown that a certain arrange- 
ment of the magnetic moments in domains leads 
to a lower energy than some other arrangements; 
but it has not been shown that this is the lowest 
energy attainable. A complete analysis would 
determine the magnetization direction as a func- 
tion of position in such a way as to minimize the 
total energy. The problem can be stated mathe- 
matically in this form, but the differential 
equations that result are too complicated to give 
much information. Moreover the details of the 
domain arrangement are probably determined 
very largely by accidental factors, such as 
irregular strains and impurities, which have not 
been included in the simplified theory. 


Low Fields (1) 


We now understand how a specimen can be 
demagnetized and yet possess spontaneous mag- 
netization throughout domains. By what process 
does the specimen pass from this state to one of 
uniform magnetization throughout each of its 
constituent crystals? And how does the mag- 
netization become reversed when the field is 
decreased from a high value to zero and then to 
a high negative value? A detailed theory has 
been worked out only for specimens in which 
the conditions are particularly simple. We shall 
consider these cases first, and then discuss their 
bearing on the properties of ordinary materials. 


(1) MATERIALS WITH PosITIVE MAGNETOSTRIC- 
TION UNpbER TENSION 


Iron-nickel alloys of a certain range of com- 
position have positive magnetostriction. If a 
specimen of such a material is subjected to a high 
longitudinal tension, its hysteresis curve assumes 
the rectangular shape shown in Fig. 7(a). There 
is 100 percent remanence, which persists until an 
opposing field /7, is reached; then the mag- 
netization reverses in a single Barkhausen jump, 
AB or A'B’. If, on the other hand, the specimen 
is demagnetized by heating and a field is then 
applied, the magnetization increases or decreases 
reversibly with the field along the curve OP, as 
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long as the field remains less than a critical value 
/I, which is considerably smaller than H,. As 
soon as Hf becomes equal to //,, the irreversible 
jump PQ occurs. 

The discontinuous processes AB, A’B’ have 
been investigated in great detail, especially by 
Sixtus and Tonks. In order to produce the jump 
A'B’ it is not necessary to apply a field H, 
throughout the whole specimen: it is sufficient 
io have a field this great in a small part of the 
specimen, provided the field in the remainder of 
the specimen exceeds the critical value /7,. With 
this arrangement, shown in Fig. 7(b), a nucleus 
of reversed magnetization first forms at a point 
where the field is /7,, and then grows until it 
swallows up the whole specimen. The boundary 
of the nucleus of reversed magnetization moves 
forward with a finite speed v of order of magni- 
tude 104 cm, sec. 


(a) J 
A 
P 
He! 
Y 
<> > 
—e 
Hs 


(b) 


000000000 


FiG. 7. Magnetic behavior of alloys with positiye mag- 
netostriction under strong longitudinal tension o. (a) 
\lvsteresis loop (B’ABA'’B') and virgin curve (OPQB’). 
//,="‘starting field,”” H.=‘‘critical field.” (b) Formation 
ind propagation of a nucleus of reversed magnetization. 
! he main field H exceeds H, but is less than H,; an auxiliary 
oil increases the field to H, over a short part of the 
pecimen, 


We may therefore interpret the ‘‘starting 
ield” IT, as the field necessary to produce a 
‘ucleus of reversed magnetization, and the 
‘critical field’”’ J7, as the field at which the re- 
ersal, once started, will be propagated through 
he remainder of the specimen. When the 
pecimen has been demagnetized by heating, 
' contains many domains magnetized in each 
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direction, instead of a single one magnetized in 
the wrong direction as is the case on BA’. The 
nuclei are therefore already present, and the 
field need only reach the value 7, in order that 
the rightly magnetized domains may swallow 
up the wrong ones (PQ). 

The speed of propagation v is found, in the 
experiment of Fig. 7(b), to be a linear function 
of the excess of /7 over the critical value /7/,: 


v= B(H—H.,). (12) 


This is easily explained, if we assume merely the 
possibility of a displacement of the boundary S 
(Fig. 7(b)) and the presence of internal forces 
that interfere with this displacement when the 
field is less than //,. If the boundary moves 
forward, the magnetic induction changes, e.m.f.’s 
are induced, and eddy currents flow. The mag- 
netic field of these eddy currents is in such a 
direction as to oppose the motion that produced 
it, and its magnitude is proportional to the 
speed v with which the boundary moves. The 
effective field acting at the boundary is the re- 
sultant of the applied field H and this eddy- 
current field, and may thus be written //— Cv, 
where C is a constant. If the eddy currents were 
not present the boundary would move forward 
almost instantly, but because of them it cannot 
exceed the speed at which the effective field is 
reduced to //,; thus the actual speed is deter- 
mined by 


H—Cv=H.,, 


which is (12) in another form. 

The constant B or C can be computed theo- 
retically, in approximate agreement with experi- 
ment. It does not depend on the tension o, but 
the criticak field /7, does. The significance of this 
critical field may be understood better by con- 
sidering the nature of the reversible change OP. 

At O the specimen consists of domains, pre- 
sumably of the approximate size and shape 
described before; but the exact position and 
shape of any boundary will be determined by 
the irregular internal stresses that were neg- 
lected before. Following Kersten, let us idealize 
this internal stress system to a succession of 
alternate regions of longitudinal tension and 
pressure, the variation being along a direction x 
perpendicular to the axis of the specimen and to 


(13) 
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the domain boundaries. When this is superposed 
on a large externally applied tension oo, the result 
is a tension o fluctuating with x as shown in the 
o vs. x graph in Fig. 8(a). The energy y per cm* 
of an interdomain boundary is given by Eq. (8) 
with A =}\,oa*. y therefore varies with position 
of the boundary because of the variation of o 
with x. y is a minimum when the boundary is 
located at a minimum (x=0) of o, as shown in 
(b), and increases if the boundary is displaced 
from this position a distance x as shown in (c). 
The forces acting in the boundary layer are 
equivalent in their effect to a hydrostatic pressure 
p, tending to push the boundary toward the 
position of minimum energy, x=0. In fact, if 
y(x) represents the potential energy per cm? of 
boundary expressed as a function of boundary 
displacement x, then p,=dy(x)/dx; p, is zero at 
the minimum point x«=0, and reaches its maxi- 
mum value at the point x, where the slope of 
the o vs. x curve is greatest. 
. In zero field, the boundary (a) 
will be at x=0. If a small field 
IT is applied in the direction 
indicated in (c), the favora- 
bly oriented domain A tends 
to grow at the expense of its 


neighbor B. potential 
energy in a given field J/ is de- 
creased by 2J,/Idx per cm? of 
the boundary if the boundary 
is displaced dx; the field is 
therefore equivalent to a hy- 
drostatic pressure py =2J,/T 


pushing the boundary down- rm 
ward in the diagram. The 2 
equilibrium position of the gq 


boundary under a given field is 
thatat which the two pressures 
p, and py just balance each 
other, as shown in (c); this 
requires 


dy(x)/dx=2J,11. (14) 


If the field is slowly increased, the boundary 
is gradually displaced, adjusting itself to the 
equilibrium position at each value of //, and 
returning through its previous positions if // is 
again decreased, as long as it does not reach the 
position x,.. This gives us the reversible mag- 
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netization curve OP of Fig. 7. But when J) 
reaches a field H, corresponding to displacemen 
x., the equilibrium becomes unstable: for a very 
slight further displacement decreases instead of 
increasing the opposing pressure p,, and the 
boundary moves on without opposition, at a 
speed limited only by the eddy currents, until B 
is completely absorbed. This gives us the irre- 
versible jump PQ. 


(2) MATERIALS WITH NEGATIVE MAGNETOSTRIC- 
TION, UNDER H1GH TENSION 


An entirely different magnetization process 
occurs if longitudinal tension is applied to a 
wire whose saturation magnetostriction \, is 
negative, for instance nickel. In this case the 
tension causes all the domains to be magnetized 
in a plane perpendicular to the axis of the speci- 
men, though at random orientations in this plane. 
With relation to a longitudinal field the domains 


(b) H=0 


(c) (d) —sH=H, 


Fic. 8. Reversible and irreversible boundary displacement. (a) Stress distribu- 
tion resulting from superposition of constant applied stress oo on irregular internal 
stresses. (b) Position of the boundary with no field. (c) Reversible displacement 
of the boundary for 7 <H,. (d) Barkhausen jump for /7 = H,: the boundary, upon 
reaching x,, is propagated irreversibly through the remainder of region B. 


are all equivalent, and there is therefore no 
tendency for one to absorb another; the only 
effect of the field in this case is to rotate the 
magnetic moments of all the domains simul- 
taneously so that they are at an angle 6<7/2 to 
the field direction; the resulting magnetization, 
J=J,cos 6, can be found as for crystals by 
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minimizing the total energy. The conclusions of 
the theory have been verified experimentally by 
Becker and Kersten, who also extended the rota- 
tion theory to the case of domains in which the 
stresses are of internal origin and vary from one 
domain to another. Although this rotation theory 
leads to correct orders of magnitude for many 
quantities, it now seems unlikely that the rota- 
tion process plays more than a minor role in 
ordinary ferromagnetic materials. 


(3) ORDINARY MATERIALS 


There seems little doubt that the two processes 
described in (1)—reversible and irreversible dis- 
placements of the boundaries between domains— 
are the ones by which magnetization proceeds in 
the steep part (J) of the magnetization curve of 
ordinary ferromagnetic materials. Boundaries 
are displaced reversibly in regions where a 
stable equilibrium is possible under the joint 
action of external and internal forces. They 
move irreversibly, with a speed limited by the 
eddy currents, when a point of instability is 
reached. Both these processes may be observed 
in the surface patterns of Fig. 4 when a field is 
applied with a component in the direction of the 
strips. As the field is increased, alternate strips 
widen and the ones* between them narrow. 
The process consists partly of a slow displace- 
ment that can be made to run backward by 
decreasing the field again, and partly of sudden 
jumps after which the boundary does not return 
to its previous position if the field is decreased. 

Remanence and coercive force depend on the 
irreversible jumps. The problem is much more 
complicated than the simple case considered in 
(1), and the quantitative theory is still in its 
infancy. 

Reversible processes are easier to handle 
theoretically, and here more quantitative results 
have been obtained. There are some properties 
of ferromagnetic materials that are, to a first 
approximation, dependent only on the mag- 
netization J and not on the procedure by which 
this value of J was reached, so that they are 
the same for instance at points P and Q, Fig. 1. 
This is true of the ‘reversible susceptibility,” 
which is the ratio AJ/AH measured in a small 
change AH carried out in the opposite direction 
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to that by which the original state was ap- 
proached—as indicated by the small arrows PP’ 
and QQ’, Fig. 1. This susceptibility is called 
“reversible’’ because such a change is com- 
paratively free from discontinuous Barkhausen 
jumps. In dealing with properties of this sort, 
one can obtain fairly satisfactory theoretical 
results by ignoring the irreversible properties and 
considering a hypothetical reversible specimen 
to which the actual specimen is, as far as these 
properties are concerned, equivalent. 

In the idealized specimen of Fig. 8 we may 
expand the function y(x) in a Taylor’s series, 


(15) 


with the first-order term absent and with the 
coefficient of x? positive because y is a minimum 
at x=0. Then (14) gives 


=2J,H, (16) 


and in a small reversible displacement 6x pro- 
duced by a change of field 677, 


(17) 


where the omitted terms are small if x is small. 
Hence for x not too large, 


bx=k-2J,6H, (18) 


where k is approximately constant. If the total 
boundary area in unit volume is S, the change in 
magnetization is 


6J Six =4kST 25H. (19) 


We shall now follow Kondorsky in assuming 
that, because of the irregularity of the internal 
stresses, the shapes of the domains are irregular, 
and Eq. (19) holds with a constant average 
value of k throughout the whole magnetization 
process. 

The area S will now depend on the relative 
volumes v4 and vg occupied by domains of the 
two types A and B, Fig. 8(c). It vanishes when 
either of these vanishes (saturation in either 
direction), and may be expected to be a maxi- 
mum when they are equal (demagnetization). 
Kondorsky assumes the relation 


S=svavpg, S=constant, (20) 
which, though impossible to justify rigorously, is 


simple and satisfies the conditions just stated. 
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Since the magnetization is given by J= J,(v4—vp), 
and since v4+vg,=1, we have v4=(1+J/J,)/2, 
vpe=(1—J, J.) 2, and Eqs. (19) and (20) lead to 
the following relation between the reversible 
susceptibility 6/, 6/7 and the magnetization J: 


2). (21) 


By a similar procedure in the more general 
case, Kondorsky is able to obtain formulas for 


08 
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Fic. 9. Reversible susceptibility x, of polycrystalline 
cobalt. Experimental points (M. Samuel): open circles, 
virgin curve; solid circles, “ideal” or ‘‘anhysteretic”’ curve; 
circles black at bottom, descending branch of hysteresis 
curve; black at top, ascending branch. Theoretical curves: 
1, boundary displacement, directions of easy magnetization 
isotropically distributed ; 2, boundary displacement, correct 
crystalline symmetry; 3, artificial model with permanent 
internal fields, but with correct crystalline symmetry. xo is 
the initial susceptibility. 

the reversible susceptibility of various types of 
crystal. No experimental data are available with 
which to test the legitimacy of Kondorsky’s 
assumptions. It is possible, however, to extend 
the analysis to polycrystalline material by 
making use of the concept of the ‘equivalent 
reversible specimen.” 

For a reversible specimen, J is a single-valued 
function of /7 which can be found by integrating 
(21); this gives for the specimen of Fig. 8 or for 
a cobalt crystal 

J=J,tanh ksJ,/1. (22) 


If J/ is inclined at an angle @ to the direction of 
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easy magnetization, JJ must be replaced by its 
component /7 cos @ along this direction, and the 
whole expression must be multiplied by cos @ to 
get the component of magnetization in the field 
direction; then for a polycrystalline specimen, 
the result must be averaged over all orientations 
of the crystals. This leads to an expression 
J=f(l1) for the magnetization curve of a 
fictitious reversible polycrystalline specimen of 
cobalt. Its derivative x=f'(//) is the suscepti- 
bility. These may be regarded as parametric 
equations determining x as a function of J. 
According to the hypothesis of the equivalent 
reversible specimen, the reversible susceptibility 
x, of an actual specimen is this same function 
of J, although the parameter // no longer has 
any simple physical meaning. 

In Fig. 9 the predictions of this theory 
(curve 2) are compared with experimental data 
for polycrystalline cobalt. Although the agree- 
ment is far from perfect, the curve fits the data 
distinctly better than does curve 1, which is 
calculated from a formula proposed long ago by 
Gans. Gans’ formula represents a limiting case, 
a crystal with many directions of easy mag- 
netization distributed uniformly over the unit 
sphere; it holds fairly well for polycrystalline 
iron and nickel. 

Curve 3 is calculated by a different theory, 
in which the symmetry of the crystals is taken 
into account but a different mechanism of 
magnetization is assumed. The close agreement 
between curves 2 and 3 and the wide deviation 
between them and curve 1 show that the crystal- 
line nature of the material has a more decisive 
influence on the magnetic properties than does 
the specific mechanism of magnetization. Two 
conclusions follow from this. The first is the 
rather discouraging one that an agreement 
between theory and experiment does not neces- 
sarily prove the correctness of the mechanism 
assumed in the theory. The second is the more 
cheering one that a knowledge of the exact 
mechanism of magnetization is not necessary in 
order to carry out an approximate calculation of 
magnetic properties. 

A number of properties have been calculated 
fairly successfully by a theory in which no 
mechanism is assumed, but the domain structure 
is taken into account in a very general way. Let 
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us suppose that we observe a ferromagnetic 
specimen through a large number of peep-holes, 
each small compared with the domains, and far 
enough apart so that what we see through one 
does not help us much in guessing what we shall 
see through another. By averaging over the 
peep-holes we shall get a good representation of 
the average magnetization, magnetostriction, 
and so on of the specimen as a whole. If the 
domain arrangement is governed largely by 
irregular internal forces, statistical methods may 
be used to determine the relative popularity of 
the various directions of easy magnetization as 
sampled through the peep-holes, and thus to 
calculate some of the magnetic properties. His- 
torically this method is older than Kondorsky’s, 
but all the results obtained by it may also be 
obtained by combining his assumptions with the 
hypothesis of the equivalent reversible specimen. 
Fig. 10 shows how well this theory predicts the 
dependence of magnetostriction on magnetiza- 
tion in polycrystalline nickel. Agreement as 
good as this is exceptional; in most cases small 
systematic deviations indicate that there are 
agencies at work that have not been taken into 
account in the theory. 


Conclusion 


An attempt has been made here to outline 
the methods of attack that have proved useful 
in dealing with a very complex phenomenon. 
To that end, the essential ideas have been 
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Fic. 10. Variation of magnetostriction \ with mag- 
netization, for polycrystalline nickel at 6°C. The circles 
are experimental values (D. Kirkham). The curve is cal- 
culated from a formula that can be derived either by the 
boundary displacement theory or by a statistical method. 
Xo is the saturation magnetostriction. 


presented and some of the simpler relations 
derived without pretense of completeness or 
rigor. Unfortunately a comprehensive presenta- 
tion of everything that has been done would still 
be sadly lacking, at certain pints, in both 
those qualities. It seems quite certain, however, 
that the basic concepts are correct: though a 
mathematical theory of the hysteresis curve has 
not yet arrived, domains are here to stay. 
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Glass as an Electrical Insulator’ 


BY C. J. PHILLIPS 


Introduction 


T does not seem to be generally recognized 

that glass is one of man’s most versatile 
creations. In almost every phase of our modern 
civilization glass plays a role of ever-increasing 
importance. The glass windows in our homes 
permit a combination of light with shelter which 
is insensitive to the caprices of nature. The in- 
candescent bulb has banished darkness, has 
given us a twenty-four-hour day and the means 
to double the constructive hours of our lives. 
\Ve receive our milk in glass bottles; we preserve 
our food in glass jars. We drink from glass 
tumblers and snuff out our cigarettes in glass 
ash trays. Glass is also the tool of the scientist 
and makes possible his test tube and microscope, 
his thermometer and spectroscope. With their 
aid he explores the innermost fastnesses of matter 
and the farthest reaches of space. 

Among the most important uses of glass are 
those in which it functions as an electrical insu- 
lator. To the incandescent lamp, as one example, 
can be added neon tubing, radio tubes, radio 
insulators, fuse plugs, lightning arrester bodies, 
and insulators for the transmission of electrical 
power over relatively long distances at high 
voltage. In some of these applications the virtue 
of transparency is added to that of insulating 
ability. In others, great dielectric strength or 
low coefficient of expansion or availability in a 
wide range of colors is of paramount importance. 
In every case, not one but several distinctive 
properties of glass are utilized, and to obtain 
these it may be necessary to employ glasses of 
widely differing compositions. 

This last point is significant. To the layman, 
glass is glass. Much of his experience is with 
windows, bottles, or tumblers; and he naturally 
issociates the word ‘“‘glass’’ with those articles. 
Since these account for much of the total tonnage 


* Modified from addresses to Pittsburgh Section, 
\. 1. E. E., and Electrical Section, E. S. W. P., Pittsburgh, 
‘ennsylvania, February 14, 1939, and to Wisconsin 


“ngineering Conference, Milwaukee, Wisconsin, March 15, 
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manufactured today, and since these glasses are 
indeed very similar in composition and charac- 
teristics, it is not surprising that one is inclined 
to think of all glasses as being similar. Yet, this 
is not the case. Corning Glass Works, for ex- 
ample, now melts on a commercial scale more 
than three hundred glasses, each different from 
all of the others in composition and physical 
characteristics. These glasses are deliberately 
formulated by the chemist so that each may be 
available for use where its particular combination 
of properties is most useful. It is important to 
have some knowledge of these properties ; other- 
wise, a discussion of insulating ability alone 
would be casual and incomplete. Accordingly, 
the first part of this paper is devoted to a brief 
discussion of the chemical and physical charac- 
teristics of various glasses, with particular em- 
phasis upon the range of values obtained from 
different glass compositions. The second part of 
the paper is concerned particularly with the 
electrical properties of glass, in bulk, as dis- 
tinguished from fiber glass. The last section dis- 
cusses the outstanding characteristics of fiber 
glass electrical insulation. 


I. General Properties of Glass 


Almost all of the glasses in use today contain 
silica as the primary constituent. In fact, if we 
could easily melt silica itself into a viscous liquid, 
free from bubbles, and could then fabricate 
articles from it, we would have an extremely 
useful material. Unfortunately, it is very difficult 
to do this. It is necessary to flux the silica with 
other materials. One of the more important of 
these is sodium oxide. So far as melting and 
working are concerned, a very good glass can be 
made from silica and soda alone, but this glass 
is too soluble in water to be of any commercial 
value. Various things can be added to overcome 
this defect; lime is the cheapest of these. Win- 
dows, tumblers, milk bottles, are all essentially 
soda-lime-silica glasses. There are many other 
materials which can be added as fluxes—any of 
the alkali metal oxides, for example. 
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After the silica has been properly fluxed, other 
materials may be introduced in a wide variety 
of proportions to change radically the properties 
of the simple silicate mixture. Thus, it is found 
that the addition of soda and lime has increased 
the thermal coefficient of expansion of fused 
silica by a factor of nearly twenty. On the other 
hand, borosilicate glasses often contain no lime 
and very little alkali, have one-third the coef- 
ficient of expansion of soda-lime glasses, and 
have very good resistance to water attack or 
attack by acids and alkalis. Other properties 
can be changed by a similar procedure, and the 
range of possible glass compositions is so large 
and so varied that it is frequently possible to 
vary one characteristic and have all others 
essentially unchanged. Despite the many per- 
mutations and combinations available to him, 
the glass chemist must at the same time contrive 
to keep the batch cost low, must mix materials 
which will melt at the lowest possible tem- 
perature, and must be sure that his glass will 
not crystallize upon cooling. The compositions 
which can be used are strictly limited by the 
tendency of all glasses to crystallize (devitrify), 
thus ceasing to be glass and losing the properties 
which distinguish glass from other materials. In 
general, glass may therefore be considered a 
homogeneous solution, usually obtained by 
melting together a large quantity of SiO. with 
smaller quantities of other alkaline oxides, alka- 
line earth oxides, or trivalent oxides, this solution 
being so cooled that devitrification does not take 
place. | 
, The range of properties obtainable from vari- 
ous glasses is shown in Table I. A few words 
about some of these properties may be illumi- 
nating. 

Refractive index.-Water is somewhat below 
the lower figure in the range 1.47 to 1.96, and the 
diamond is far above the higher limit. However, 
some glasses probably exceed the precious stone 
in dispersive power. 

Specific gravity.—The lightest glasses have a 
specific gravity lower than that of fused quartz. 
The ordinary soda-lime glasses weigh about as 
much as metallic aluminum. The heaviest weigh 
more than gray cast iron. 

Elasticity coefficient——The resistance to dis- 
tortion under stress (Young’s modulus) can be 
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TABLE I. Range of physical values obtainable from 
various glasses. 


PROPERTY RANGE OF VALUES 
Refractive index 1.467—1.962 
Specific gravity 2.125-8.120 


Elasticity coefficient 
Compressive strength 
Tensile strength 
Thermal conductivity 
Expansion coefficient 
Softening point 
Annealing point 
Volume resistivity 
Dielectric constant 
Power factor 


6,500,000—12,500,000 Ib. /in.? 
90 ,000-- 180,000 Ib. /in2 
4,000— 1,000,000 Ib. /in.? 

0.0018-0.0028 cal./cm °C/sec. 

8x 10°-7-140 10-7 cm/em °C 

500°C-975°C 

350°C-725°C 

10°-10'® ohms/cec 

3.7-16.5 

0.000117-0.00664 


varied in nearly a 2 to 1 ratio, from a value lower 
than that of most metals to a value higher than 
that of aluminum. 

Com pressive strength. ~The resistance to crush- 
ing possessed by glass is tremendous—-greater 
than that of brick, concrete, or cast iron. In 
general, glass will maintain its shape under 
pressures which cause metals to flow like putty. 
This property can often be put to good engineer- 
ing use. 

Tensile strength..-The great range shown in 
Table I depends not so much on composition as 
on the condition of the exposed glass surface. 
The work of Griffith' on glass and of Joffé? on 
crystals indicates the probability that the true 
tensile strength of a material is determined by 
the molecular cohesion, and that the lower values 
found in practice are caused by surface imper- 
fections, some of which may be inherent in the 
material. For glass, the modulus of rupture also 
depends upon the size and shape of the object, 
although both of these may be dependent vari- 
ables, the kind and degree of surface imperfection 
being the independent variable. For example, 
glass in the form of threads or fibers of small 
diameter is tremendously strong, ranging upward 
from 180,000 Ib./in for fibers 0.0004” in 
diameter. More is said of this in a later section 
of this paper. 

Glass can be tempered in a manner analogous 
to the tempering of steel, although the physical 
processes occurring in the glass due to this 
treatment are quite different from those in the 
metal. By tempering, it is possible to prepare 
glasses which have a tensile strength of 50,000 
lb. /in.*, about equal to the strength of wrought 
iron. 
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Expansion coefficient.—Ordinary soda-lime glass 
has an expansion of about 9010-7 cm/cm °C, 
whereas the values for iron and platinum are 
120X10-7 cm/em °C and 9010-7 cm/em °C, 
respectively. ‘Pyrex’’ brand electrical glass has 
an expansion of 3210-7 cm/cm °C. Other 
“Pyrex”” brand glasses have expansions of 
46-50) X 10-7 cm/cm °C and seal to tungsten 
or molybdenum, whose values are 45X1077 
°C and 5210-7 cm/cm °C. The new, 
low expansion, 96 percent silica glass developed 
by Corning Glass Works has a coefficient of ex- 
pansion of 810-7 em/em °C and compares 
favorably with (5.5-5.85) cm/cem °C for 
fused quartz. 

Annealing point.—The annealing temperature 
of glass—the temperature at which strain is 
relieved —is 350°C to 725°C, considerably higher 
than that for steel, which is between 200°C and 
500°C, 


cm, cm 


A comparison of some of the electrical proper- 
ties of various materials is given in Table II. 
These electrical properties will now be discussed 
in greater detail. 


Il. Electrical Properties of Glass 
DIELECTRIC ABSORPTION 


The electric current which flows through glass 
under a given impressed voltage is not a simple 
function of the temperature, physical dimensions, 
and conductivity of the sample. Instead, it 
requires appreciable time before the glass soaks 
up and becomes saturated with the electric 
charge. Upon discharge, the reverse may occur. 
The current flowing at any particular time is, 
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Fic. 1, Resistance of dry glass as a function of the time 


after applying voltage. (After Bronson.) 


therefore, an undetermined function of the total 
elapsed time since the initial application or 
removal of the voltage. This phenomenon, gener- 
ally called ‘‘dielectric absorption,’”’ is not confined 
to glasses. It appears to be present to some 
degree in all solid dielectrics. Both theory and 
experiment show that the alternating losses in 
such materials are due almost entirely to this 
effect. 


TABLE II. Electrical proportion of insulating materials. 


THICKNESS 
MATERIAL 


DIELECTRIC STRENGTH f 


DIELECTRIC VOLUME RESISTIVITY 


| 

(MM) | (KV/CM) CONSTANT | (OHM-—CM) 
|. Asbestos paper | 1.2 42 — 1.6 10" 
2. Bakelite, Micarta—213 — 314 5.0 | 5.0 10" 
3. Cellulose acetate 0.019 480 | 4.0— 6.2 10'2-10" 
+. Empire cloth, muslin 0.38 480 3.5- 5.5 -- 
5. Ethyl cellulose | _- 2.7- 3.3 10" 
6. Fiber, vulcanized 3.2 108 5.0 101° 
7. Glass, various kinds 0.50 700 3.7-16.5 108-10'8 
Methyl methacrylate | 3.3- 4.5 3.0 10" 
9. Mica | 0.60 280 5.0- 7.0 2.0 x 10"7 
‘0. Paraffine (Parowax) | — 115 | 2.1 10'¢ 
|. Porcelain, unglazed 2.5 130 5.7- 6.8 3.0 10" 
2. Rubber, hard 0.50 | 700 2.0- 3.5 10'8 
3. Sulphur — — 2.9- 3.2 
‘4. Wood—maple, paraffined 15.20 | 46 | 4.1 3.01018 
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Fic. 2. The surface conductivity of Pyrex resistant glass 
at 25°C and at 50°C for various frequencies as a function of 
relative humidity. (After Yager and Morgan.) 


In making measurements on the electrical 
resistance of glass, the dielectric absorption must 
be taken into account. The absorption current 
may predominate over the conduction current 
when the resistivity is high, and measurements 
must be taken at different time intervals until 
the absorption current has become negligible. 
The necessity for this is shown in Fig. 1.* 


SURFACE RESISTIVITY 


kor many years it has been understood that 
part of the electrical conductivity of glass is due 
to a film of moisture condensed on the surface. 
The existence of this film was conclusively 
demonstrated by Warburg and Ihmori.* They 
showed that with the better glasses the film 
could be removed in a dry atmosphere, whereas 
with others it could be removed only by heating. 
Many other observers have shown that an 
increase in the alkali content is reflected in an 
immediate decrease in the surface resistivity. It 
is problematical whether the film is continuous 
or spotty and whether it is wholly adsorbed or 
partly absorbed. It appears certain, however, 
that the film is an aqueous, alkaline solution 
whose properties may vary widely. At high 
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humidity or with glasses rich in alkali, the 
surface conductivity may greatly exceed the 
volume conductivity. On the other hand, it is 
entirely negligible for alkali-poor glasses with 
absolutely clean surfaces if the glass is properly 
polarized. Almost any clean glass can be used 
without serious surface leakage if the tem- 
perature is raised a few degrees above that of 
the surrounding air. 

Yager and Morgan® made a special study of 
the effect of both frequency and humidity on 
the surface conductivity of ‘Pyrex’ brand 
laboratory glass. They found that with low 
humidities, constant values could be obtained in 
two hours; but at 50°C and 95 percent humidity 
it was impossible even to approach equilibrium. 
Their results are plotted in Fig. 2. 


VOLUME RESISTIVITY 


The volume resistivity of glass depends upon 
and varies rapidly with changes in composition 
and temperature. At all temperatures glasses are 
electrolytic conductors, and most investigators 
ascribe all of the conduction to the sodium ion. 
This is probably largely, but not wholly, correct. 
It does not, for example, explain the role of 
calcium at higher temperatures, nor the mecha- 
nism of conduction in alkali-free glass. Much 
fundamental work remains to be done to explain 
these and other perplexing questions. 

Curtis® found the volume resistivity of fused 
silica to be about 5X10'* ohm-centimeters at 
22°C. At ordinary temperatures, ‘‘Pyrex’’ brand 
electrical glass has a resistivity of about 10", 
plate glass about 10", and certain special glasses 
will be as low as 10%. It will be seen that at 
ordinary temperatures some glasses are amongst 
the very best insulating materials we possess. 

Figure 3, reproduced from Fulda,’ shows the 
change in specific resistance of a Thuringian 
glass with temperature. This is the only case in 
which the resistance has been determined from 
room temperature to 1200°C. Littleton* pub- 
lished data on the resistivity of several widely 
different glasses over a considerable range of 
temperatures, using data obtained by Babcock.* 
His results are shown in Fig. 4. These scattered 
data are sufficient to show the tremendous effect 
on volume resistivity ot changes in glass com- 
position and temperature. 
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DieLectRic CONSTANT, POWER FACTOR, 
DieELEcTRIC Loss 


Glasses are unique among solid dielectrics in 
that they possess a range of dielectric constants, 
K, varying by more than four to one, from about 
3.7 to 16. This property is of more than purely 
theoretical interest. Where maximum capacity 
for a given size of condenser is required, a high 
dielectric constant is needed. Where high voltage, 
high frequency condensers are required, then a 
low dielectric constant and also a low power 
factor, 6, are indicated. The ordinary high 
voltage power insulator is another advantageous 
application of low dielectric constant material, 
because it can be shown that the smaller K is, 
the greater will be the voltage, V, necessary to 
form corona in the air in series with the insulator. 

For most glasses, K decreases only slightly 
with increasing frequency—4} percent as the 
frequency was varied from 500 to 10° cycles, 
according to MacLeod. K increases with in- 
creasing temperature ; more so at low frequencies 
than at high. 

The power loss per unit volume of dielectric 
material is directly proportional to the fre- 
quency, the square of the voltage gradient, and 
the product K6. For conditions of high frequency 
and or high voltage, the power loss is appreci- 
able, and Ké should be kept as small as possible. 
lor glasses, both K and 6 are to a certain extent 
separably alterable by changes in glass com- 
position, and the product Ké can be varied in 
a ratio of more than 200 : 1. Two general classes 
of glass are available: (1) high K and low 6; and 
(2) low K and low 6. As an example, for Corning 
glass 707 at 10* cycles, McDowell, and Begenian® 
find: K=4.1, 6=0.00077, and KéxX10'=31. 


Temperature in C — Scale in 


SY 


Specific Resistance — Logarithmic Scale 


Fic. 3. The change in specific resistance of a glass with 
temperature. (After Fulda.) 


This low loss glass is a low alkali borosilicate of 
low expansion and high surface stability—an 
ideal insulating material for most purposes. 
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The results on the variation of power factor 
with frequency present several anomalies and, 
therefore, are not in wholly satisfactory accord. 
However, it appears that for most glasses the 
power factor decreases with increased frequency. 
Power factor increases rapidly with rise in tem- 
perature. As the temperature increases, dielectric 
heating also increases, with eventual instability 
and breakdown. This condition of instability is 
normally reached only in the case of high field 
strength associated with high frequency. 
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Fic. 4. The change with temperature of the resistivity of 
various glasses. (After Babcock.) 


DIELECTRIC STRENGTH 


The dielectric strength of any material is the 
voltage gradient required to puncture that ma- 
terial. It is usually expressed as kilovolts per mil 
or kilovolts per centimeter. If the dielectric 
strength is a true physical constant, it should be 
independent of the thickness of the material. 
Nevertheless, in most cases it is found that the 
greater the thickness, the lower is the measured 
dielectric strength. It is not known whether this 
effect is inherent in the material or is dependent 
upon the test method. The observed values of 
dielectric strength vary with the size and shape 
of the electrodes used in applying the voltage, 
and also vary with the frequency and duration of 
the voltage, the temperature, etc. 

There has been very little work done on other 
solid materials having dielectric strengths com- 
parable to those of most glasses. Although much 
attention has been paid to the elimination of elec- 
trode edge discharges, many of the other vari- 
ables mentioned above have been neglected. It is 
probable that the true dielectric strength of glass 
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is even greater than it appears to be. Even so, 
the measured values are very high. Most of the 
work has been done on sections less than 0.50 mm 
thick. The values are usually expressed as peak 
values in the disruptive region with edge effects 
eliminated. Moscicki' obtained a constant value 
of 1360 kv, cm for various thicknesses up to 0.50 
mm. Inge and Walther" derived a constant value 
of 3000 kv cm on certain Russian apparatus 
glasses whose thicknesses varied from 0.05 to 
0.50 mm. Rochow™ found a flint glass to have a 
dielectric strength of 4000 kv, cm at 0.05 mm. 
Moon and Norcross" found values as high as 
3100 kv, cm for a lead glass. 

It is almost impossible to avoid electrode edge 
discharges at the higher voltages required to 
puncture thick sections of glass, and conse- 
quently few results have been published. How- 
ever, there is evidence that puncture strengths 
above 1900 kv/cm can be reached in air on 
1.0-mm sections of “Pyrex” electrical resistant 
glass insulators. The impulse strength of these 
insulators is such that one-half-inch sections have 
repeatedly withstood discharges of a “‘lightning 
generator” at 1500 kv. 


Il. Fiber Glass 


Glass fibers are not new. They have been 
produced for many years on a small scale—first 
being drawn by hand, later on a large wheel: 
hence the name “spun glass.’’ New methods 
developed in this country only a few years ago 
finally made it practicable to produce exceedingly 
fine glass fibers on a commercial scale. 

Two types of fibers are of particular interest : 
staple fibers, 4° to 15” long, and perhaps 
0.00025” in diameter; and continuous fibers, 
which, as the name implies, are drawn as long, 
continuous filaments, usually about 0.0002” in 
diameter. The production of these fibers is in 
pounds per day rather than tons. Even so, one 
glass marble weighing about one-third of an 
ounce could be melted and drawn into a con- 
tinuous fiber 0.0002” in diameter and over one 
hundred miles long. 

Both types of fibers are equally good elec- 
trically. However, the continuous type is 
smoother, slightly better mechanically, and 
available in smaller diameters. Staple fiber 
products are generally better where a heavy type 
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insulation is necessary. The continuous type has 
an advantage in reduced space factor and better 
appearance. 

Fiber glass differs from ordinary lime glass in 
one fundamental way—it is substantially alkali- 
free. This freedom from alkali is essential because 
alkali leaches to the surface and under conditions 
of high humidity forms. ionizable salts which 
decrease the surface resistivity. In a fiber 0.0002” 
in diameter and 1” long, the ratio of surface area 
to volume is 20,000 : 1, whereas in a piece of 
glass 1’ square and }” thick this ratio is only 
12 : 1. The effect of alkali leaching is magnified 
proportionally. 

Glass fibers have a number of astonishing 
properties. In many cases these are substantially 
the same for both types described above. When 
not otherwise specified, the data given below 
may be assumed to apply to either type. 


TENSILE STRENGTH 


In discussing the tensile strength of glass 
fibers, one deals with figures of almost astro- 
nomical magnitude. The tensile strengths for 
continuous fibers are approximately as follows: 


0.0002” diameter 
0.0003” diameter 
0.0004” diameter 


400,000 Ib. /in.2 
250,000 Ib. 
180,000 Ib. 


These figures should be compared with a rough 
figure of 70,000 Ib. in.* for mild steel. The reason 
for the high strength is closely associated with 
the theory of the rupture of brittle bodies, first 
clearly enunciated by Griffith. With extremely 
fine fibers, carefully drawn, it is possible to reach 
specific tensile strengths in excess of 1,000,000 
lb./in.2 No other textile material approaches this 
figure. 

The flexibility of the fibers is due to their small 
diameter, not to any peculiarity of the glass com- 
position. The high tensile strength allows them 
to be bent through a very small radius without 
breaking. The very fact that glass fibers can be 
spun, twisted, and woven on standard textile 
machinery into fabrics as thin as 0.0025” is a 
good illustration of their flexibility. 

When fibers are twisted into yarns, some of 
the effective strength of the fiber is lost, for 
several reasons: (1) lack of parallelism of fibers, 
(2) non-uniformity of stress distribution, (3) slip- 
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page of fibers, and abrasion, and (4) bending or 
crushing stresses, particularly in tightly woven 
cloths. Despite these characteristics, the finer 
varns may have strengths of approximately 
65,000 Ib./in.? and hence, weight for weight, are 
stronger than yarns of most organic fibers, ex- 
cepting manila and linen. 

When yarns are further fabricated into cloths, 
tapes, or braids, there is a further loss of strength. 
In open weaves the strength may be 80 percent 


lance II. Tensile strength in pounds (minimum-average- 


THICKNESS STAPLE FIBER CONTINUOUS FIBER 
INCHES) 4” Wipe 1” 4” Wipe 1” Wipe 


0.005 133-147-172 236-258-280 
010 72— 91-110 153-178-200 229-257-285 312-452-510 
O15 110-134-158 177-230-263 346 595-634-674 
020 134-158-183 213-259-305 — _ 


that of the yarn, in close weaves only 50 percent. 
Table II] shows the breaking strength of some 
of the staple fiber and continuous fiber glass 
tapes used for electrical insulation. In Table III 


TABLE IV." Flexing endurance. Samples: 
Cloth, 0.010" X17". 


YELLOW VARNISHED YELLOW VARNISHED 


TYPE COTTON GLASS 
\s received 2100 2700 
\fter 280°F, 24 hr. 12,000 72,000 
\fter 350°F, 2 hr. 4500 80,000 


the values given are the minimum, average, and 
maximum obtained by combining the results 
obtained from several sources. 

The effect of heat on the tensile strength of 
asbestos, cotton, and fiber glass is shown in Fig. 
5. These curves show that fiber glass-tapes are 
from two to twenty times stronger at room tem- 
perature than commonly used fibrous materials. 
Of equal importance is the fact that the strength 
of glass tape decreases only about 20 percent up 
to 600°F and after twenty-four hours at 800°F 
still retains half its original value. 


RESISTANCE TO ABRASION AND FLEXING 


No conclusive tests have been made to deter- 
mine the abrasion resistance of glass tapes and 
fabrics because there is no unanimity of opinion 
as to what constitutes a proper test. In general, 
however, it seems that there are no outstanding 
differences in this respect between either treated 
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Fic. 5. Effect of heat on tensile strength of electrical 
tapes. (Atkinson.) 


or untreated glass tapes and equivalent cotton 
or asbestos tapes. 

In flexing endurance, glass fiber insulation 
seems to possess a considerable advantage, as 
Atkinson’s" data shows (Table IV). These values 
were obtained on a machine designed to create 
a “whipping” motion and show the total number 
of flexes required to wear out each sample. 


HYGROSCOPICITY 


Another interesting characteristic of these 
fibers which makes them particularly suitable 
for electrical insulation is that they are practi- 
cally nonhygroscopic. Their moisture absorption 
is lower than for any other fibrous insulating 
material. The moisture gain on glass fibers in 
passing from 0 percent to 90 percent relative 
humidity at room temperature will be { of 1 
percent by weight, or less. This small amount 
of moisture is probably adsorbed, rather than 
absorbed, and so does not change the volume 
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characteristics of the glass. It will have an effect 
on the surface characteristics, although much less 
than for organic materials. 


VOLUME REsISTIVITY 


Largely because of the absence of soda and 


potash, the volume resistivity of alkali-free glass 
is very high—at 490°C (914°F) about 2X10? 
times greater than for the lime glass discussed 


——MEASURED VALUES 
——-—EXTRAPOLATED VALUES 
v 7 T T 
| 
NN 
| 
\_ | 
| T T 
z \ } 
+ + + 
T + \ 
x 
> \ 
| \ 
a we ] \ 
\ 
| | | 
TT 
| | 
aw 
[| ~ 
0 400 2400 


800 1200 1600 2000 
TEMPERATURE— DEGREES FAHRENHEIT 


Fic. 6. Volume resistivity of alkali-free glass. (Atkinson. ) 


hy Littleton’ and times greater than for 
his borosilicate glass. Fig. 6, from Atkinson," 
shows the volume resistivity of the alkali-free 
glass for temperatures up to 1315°C (2400°F). 


SURFACE RESISTIVITY 


The resistance of glass fiber insulation at 
ordinary temperatures can be considered almost 
entirely due to its surface resistance. This is 
because of the high ratio of surface area to 
volume or weight. When dry and untreated, this 
resistance is very high—many times greater 
than that of ,cotton or asbestos. This surface 
resistance is, of course, affected by the relative 
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humidity of the atmosphere, but to a lesser 
extent than for most other materials. The surface 
resistance decreases with increased humidity, as 
Table V"* shows, and also with increased time of 
exposure. 

The resistance of treated fiber glass varies with 
the nature of the impregnant and the method of 
applying it. When dry, the glass product behaves 
very much like other porous, woven materials 
similarly treated. However, when exposed to 
high humidity, the glass retains a much greater 
part of its original insulating ability as Fig. 7" 
shows. 

DIELECTRIC STRENGTH 


The dielectric strength of fiber glass products 

as is the case with any fibrous material— 
depends largely upon what fills the interstices of 
the weave. For untreated products, the dielectric 
breakdown is nearly the same as that of air. It 
is, therefore, customary to impregnate with var- 
nishes, gums, or resins. The resulting products 
have dielectric strengths very nearly equal to 
that of the impregnant. Glass is no exception to 
this rule; and specific values, therefore, need not 
be given. Glass, however, has the advantage of 
being nonhygroscopic and has a higher dielectric 
strength than materials like asbestos after ex- 
posure to moisture. 


OTHER PROPERTIES 


Fiber glass possesses other properties which, in 
combination with those already described, make 
it outstanding. Because it is glass and inorganic, 
it is noninflammable and heat resistant; to a 
large extent the operable temperature is limited 
only by the impregnant. As Ferris and Moses!® 
point out, higher temperatures are a means to 
an end: (1) reduction in size of a machine for a 
given rating; (2) reduction in weight; and (3) 
reduction in cost per unit rating. The compact- 
Taste V. Insulation resistance in megohms. Standard 

unimpregnated tapes Measurements after 
humidity held constant 72 hours at 100° F. 


As RECEIVED | WASHED 
Tyre 78% R.H. | 90% R.H. | 75% RH. | 90% RH. 
0.015” asbestos 0.57 1.08 1.50 | 5.5 
0.011” cotton 1.03 35. 30.5 | 380. 
30,100 


0.015” glass | 12.4 200. 1400. 
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ness and good thermal conductivity of fiber glass, 
compared to organic materials, assist in attaining 
the characteristics listed above. The textile 
itself, however, is not the complete answer. In 
addition, new, correlated products are required: 
insulating resins and compounds, insulating 
sheets or boards for both rotating machines and 
transformers, high temperature liquids for liquid- 
cooled transformers, and perhaps even new core 
plate vanishes. It is to be hoped that the advent 
of fiber glass will act as an incentive and a 
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N 
COTTON 
10° 


6 10 14 18 22 26 30 
NUMBER DAYS AT 90 PERCENT RELATIVE HUMIDITY 


Fic. 7. Effect of humidity on resistance of black 
varnished cloth. (Atkinson.) - 
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stimulant for the completion of this compre- 
hensive program. 


Summary 


Some of the advantages of glass as an elec- 
trical insulator may be briefly summarized as 
follows. 

(1) Glass is a homogeneous, nonporous ma- 
terial which is smooth, hard, and highly resistant 
to weathering and chemical attack. 

(2) Glass has good tensile strength and is very 
strong in compression. If a strain member can 
be designed so that the glass part is in com- 
pression, very large loads are possible. Certain 
glass suspension insulators will support twenty- 
five thousand pounds. 

(3) Certain kinds of glass have coefficients of 
expansion lower than any other ceramic sub- 
stance except fused quartz. This characteristic 
tends to minimize the effect of severe temperature 
gradients and improves resistance to power arcs. 

(4) Glass used for most electrical purposes has 
very high volume and surface resistivity. 

(5) Glasses are available with a wide range of 
values of the dielectric constant and of power 
loss per unit volume. 

(6) Glass has very high dielectric strength and 
hence great resistance to high voltage surges. 

(7) Glass is available in the form of fiber glass, 
with the numerous advantages offered by that 
material. 
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A Hydrogen Arc for 
Absorption 
Spectroscopy 


A low voltage hydro- 
gen arc, for use in the 
determination of absorp- 
tion spectra in the ultra- 
violet region was described recently by A. J. 
Allen and R. G. Franklin.' This are emits a 
continuous spectrum in the ultraviolet region, 
compares favorably in intensity with an intense 


NOV 4 
A.C. 3 INCHES 


Schematic drawing of low voltage hydrogen arc. V,;— 
filament voltage, 3.5-4.0 volts; V.—arc voltage, 40.0 volts; 
A—ammeter to measure arc current, 1.0—1.5 amperes; and 
R—variable resistance, 100 ohms, —2.0 amperes. 


spark between iron electrodes and is easy and 
reliable to operate. It may be made to take the 
form of an approximate point or slit according to 
the requirements of the spectrograph used, and 
is, Moreover, inexpensive to construct. 


'J. Opt. Soc. Am. 29, 1 (1939), 
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Resumés of Recent Researeh____ 


A schematic drawing of the essential elements 
of the arc is shown in the accompanying figure. 
The cathode consists of a strip of nickel screen, 
made into a corrugated form and dipped into an 
oxide coating material. The cathode is  sur- 
rounded by a nickel cylinder with an orifice 
which can be adapted to any size or shape. The 
anode is a rectangular piece of molybdenum with 
an orifice to correspond to that of the cathode. 
The cathode-anode assembly is placed in a 
Pyrex bulb onto which a quartz window is sealed 
with glyptal lacquer. The bulb is completely 
water-jacketed for cooling. Ordinary tank hydro- 
gen at a pressure of 4 to 6 mm of mercury has 
been found to be satisfactory for filling the arc. 
Kither alternating or direct current may be used. 
The are voltage is approximately 40 volts and the 
are current 1.0 to 1.5 amperes. The filament 
carries a current of 10 to 12 amperes at 3.5 to 4.0 
volts. 

Complete instructions for making this arc are 
given in the original article. Typical spectrograms 
taken with the are are also shown and data are 
given from which comparisons can be made with 
other light sources. 


Identification of Silicosis develops rath- 
Aluminum Hydrate er quickly in rabbits 
Films of Importance’ exposed to air containing 
in Silicosis 
Prevention 


moderate concentrations 
of quartz particles finer 
than about 5 mu, but is 
completely prevented if aluminum powder is also 
present in the air to the extent of about 1 percent 
by weight of the quartz powder.' It has been 
established that aluminum forms in the lungs a 
protective layer upon the surface of silica parti- 
cles, which prevents them from dissolving and 
thus prevents toxic effects.' The nature of this 
layer has been studied by means of electron 
diffraction.? In these experiments electron diffrac- 
tion patterns have been obtained from thin films 


1]. J. Denny, W. D. Robson and D. A. Irwin, Canadian 
Medical Association J. 37, 1-11 (1937); 40, 213-228 
(1939). 

*L. H. Germer and K. H. Storks, Ind. Eng. Chem. 
Anal. 31, 583 (1939). 
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of silica, about 200A thick, which had been 
previously exposed to aluminum and water at 
38°C (approximately body temperature). A beam 
of high speed electrons is sent through such a 
ireated film and the resulting diffraction pattern 
recorded upon a photographic plate. From studies 
of such patterns, and comparisons with x-ray pat- 
terns of known substances, materials composing 
layers upon silica surfaces have been identified. 

After treatment with aluminum and water at 
high and medium pH values (6 to 9) silica films 
vive sharp diffraction patterns which have been 
identified with oriented crystals of that hydrated 
oxide of aluminum known as aluminum alpha- 


Cover Photograph 


The photograph on the cover of this issue shows several 
sections of the 1,400,000-volt x-ray tube, constructed by 
General Electric for the National Bureau of Standards. A 
description of this tube was given in the Innovations in 
Instruments section of the February Journal of Applied 
Phystes. 


* 


Research Preview 


The Westinghouse Electric and Manufacturing Company 
held a preview of its research work in 1940 on Monday, 
January 29, in the Research Laboratories, East Pittsburgh, 
Pennsylvania. The following program was arranged: 


Morning 


Welcome, by Dr. L. W. Chubb 
Carburization of Steel, by Dr. A. A. Bates 
llydrochloric Acid Relay, by Dr. C. F. Hill 
Westinghouse ‘‘Atom Smasher,”’ by Dr. E. U. Condon 
Inspection tours of the ‘‘atom smasher,"’ mass spectrograph, and 
clasticity and damping of single crystals. 


Afternoon 

Greeting, by M. W. Smith 
\rt in modern design, by D. L. Hadley 
Vhotoelasticity, by Dr. Miklos Hetenyi 
Lightning investigations, including electrical wave model and fulchrono- 

graph, by C. F. Wagner 

Inspection trips to see: control devices; microchemical laboratory; 
exhibit room, equipped with working models; magnetic bearing tester; 
field balancing; creep, fatigue and relaxation; high speed tension; 
K-42-B demonstration; bright annealing in a new electric furnace. 


* 


Appointment of Dr. C. H. Kunsman, long active in 
research work in the U. S. Department of Agriculture, as 
Chief of the Physiochemical and Analytical Division of the 
Department's Western Regional Research Laboratory at 
Albany, California, was announced recently by Dr. Henry 
G. Knight, Chief of the Bureau of Agricultural Chemistry 
and Engineering. 
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monohydrate (Boehmite). At pH 4 monohydrate 
is not formed upon silica. In the presence of 
organic acids oriented soap crystals are formed 
at all pH values from 4 to 9, but this material 
appears not to be readily adsorbed upon silica. 
The experiments prove that aluminum hydrate 
is precipitated fairly rapidly upon silica at pH 
values lying within a range in which lie also the 
PH values of body fluids of men and of animals. 
Since in these experiments aluminum hydrate is 
not formed upon silica at pH 4, it seems highly 
probable that aluminum would not afford pro- 
tection from silicosis to a hypothetical animal 
with body fluids of pH 4. 


Fellowship Fund Announced 


A $25,000 trust fund to provide graduate fellowships in 
electrical engineering has been set up by the Westinghouse 
Electric and Manufacturing Company, the American 
Institute of Electrical Engineers announced recently. The 
fund is established as a memorial to Dr. Charles LeGeyt 
Fortescue who was associated with the Westinghouse 
Company for 38 years. Each award provides a minimum 
allowance of $500, and any student of electrical engineering 
who has received a bachelor’s degree from an accredited 
college is elegible for a fellowship, provided he does not 
hold other fellowships which carry a stipend greater than 
the tuition required by the institution at which he proposes 
to undertake his graduate work. Professor D. F. Miner, of 
Carnegie Institute of Technology, is secretary of the 
fellowship committee. 


* 


Professor Ralph A. Sawyer, who has been granted a 
sabbatical leave from his duties at the University of 
Michigan for the present school semester, wili spend most 
of that time in Ann Arbor devoting himself to study and 
to work on a book on experimental spectroscopy. 


* 
Power Conference Plans 


The Annual Midwest Power Conference, to be held in 
Chicago’s Palmer House, April 9 and 10, includes, in its 
tentative program, plans for sessions on the following 
topics: gas turbines, small power plants, electrical trans- 
mission, power-process, fuel problems of power plants, and 
hydro power. At the opening meeting the address of 
welcome will be given by Loran D. Gayton, City Engineer 
of Chicago; and Philip W. Swain, Editor of Power, will 
speak on “Power Fallacies.”” The Conference is sponsored 
by the Armour Institute of Technology. 
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Royal Society Awards 


Among the awards made recently by the Royal Society 
of Great Britain, the following are of particular interest 
to physicists: 

A Royal Medal was presented to Professor Paul Adrien 
Maurice Dirac whose work in modern theoretical physics 
has aided in the promotion and interpretation of observed 
atomic phenomena. Through careful and original mathe- 
matical treatment Professor Dirac was successful in bridg- 
ing the gap between the new quantum mechanics and 
Newtonian, and later in reconciling quantum and rela- 
tivity mechanics by making the fundamental equations 
of quantum mechanics invariant under a Lorentz trans- 
formation. This broadening of the field of quantum me- 
chanics led to the revis‘on of the theory of the hydrogen 
atom, gave added indication of a quantized electron spin, 
and hinted at the existence of the positive electron, since 
discovered by Anderson. 

‘The Hughes Medal, awarded annually by the Rovyal 
Society for original discovery in the physical sciences, 
particularly in electricity or magnetism or their applica- 
tions, was given this year to Professor George Paget 
Thomson, who has been able experimentally to prove de 
Broglie’s concepts of wave mechanics, both from a qualita- 
tive and a quantitative standpoint. Professor Thomson's 
earlier work was concerned with research in positive rays, 
and he has also published a book and a number of papers 
on aeronautics. He collaborated with his father, Dr. J. J. 
Thomson, in the preparation of the third edition of the 
book, Conduction of Electricity through Gases, the 
important work on that subject 


* 
School Merger Planned 


By September, 1940, a merger of the Armour Institute of 
Technology and Lewis Institute will have been completed, 
according to an announcement made recently by Armour 
Institute. The name of the new school will be Armour 


most 


College of Engineering of Illinois Institute of Technology. 
* 
i International Congress of Mathematicians 


Postponement of the International Congress of Mathe- 
maticians until a more favorable time was announced by 
the Organizing Committee recently. The Congress was 
originally scheduled for September, 1940, in Cambridge, 
Massachusetts. The Committee further stated that due 
notice will be given of any arrangements to hold the Con- 
gress at a later date. 


* 


Science reports that Dr. Elmer O. Kraemer, formerly 
connected with the Experimental Station of E, I. duPont 
de Nemovrrs and Company, has recently returned from 
Europe, where he spent a year as fellow of the Lalor 
Foundation, continuing his researches on the sizes and 
shapes of giant molecules with Professor The Svedberg at 
the University of Upsala, Sweden. Dr. Kraemer is now 
associated with the Biochemical Research Foundation of 
the Franklin Institute in Philadelphia. 
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Illinois Engineering Council 


Professor Harry McCormack, Director of Chemical 
Engineering at Armour Institute of Technology, was 
recently elected President of the Illinois Engineering 
Council. His election took place at the annual meeting of 
the Council held on the campus of the Institute. 

“The Illinois Engineering Council,”’ according to Presi- 
dent McCormack, ‘‘was founded to present a united front 
for the engineering professions in the State of Illinois. 
Actually, the Council consists of about forty representa- 
tives from the various engineering societies in this state, 
and the Council's function is to represent the social, 
economic and legal interests of the engineering profession 
as a whole.” 

The organizations thus far represented in the Council 
Illinois Assembly; 
American Society of Civil Central 
Section; Institute of Chemical Engineers, 
Chicago Section; Illinois Association of County Superin- 
tendents of Highway; Illinois Association of Highway 
Engineers; and the Western Society of Engineers. 


* 


Engineers, Illinois 


American 


Mathematical Biophysics 
The Bulletin of Mathematical Biophysics, edited by 


N. Rashevsky of the University of Chicago, announces the 
following table of contents for the March, 1940, issue: 


A Neural Mechanism for Discrimination: II. Discrimination of Weights. 
A. S. HOUSEHOLDER 


An Approach to the Mathematical Biophysics of Self Regulation and 
Cell Polarity. N. RASHEVSKY 

Contributions to the Mathematical Biophysics of Organic Form: 
I. Formation of Cavities in Cellular Aggregates. N. RASHEVSKY 

Discrimination between Temporally Separated Stimuli. D. 

Convective Diffusion in Parallel Flow Fields. Gate YounG 


The editorial office of the Bulletin of Mathematical 
Biophysics is located at 5822 Drexel Avenue, Chicago, 
Illinois. 


* 
Calendar of Meetings 


4-8 American Society for Testing Materials, Detroit, Michigan 


12-14. American Railway Engineering Association, Chicago, Illinois 
April 

7-13 American Ceramic Society, Toronto, Canada 

8-12 American Chemical Society, Cincinnati, Ohio : 
17-19 American Society of Civil Engineers, Kansas City, Missouri 
22-25 Society of Motion Picture Engineers, Atlantic City, New Jersey 
25-27 American Physical Society, Washington, D. C. 
29-30 Acoustical Society of America, Washington, D. C. 
May 

1-3 American Society of Mechanical Engineers, Worcester, Massa- 


chusetts 
1-4 American Geophysical Union, Washington, D. C. 


6-10 American Foundrymen’s Association, Chicago, Ihnois 
13-15 American Institute of Chemical Engineers, Buffalo, New York 
June 


6-8 National Colloid Symposium, Ann Arbor, Michigan 


9-14 Society of Automotive Engineers, White Sulphur Springs, West 
Virginia 

10-14. American Medical Association, New York, New York 

17-21 American Society of Mechanical Engineers, Milwaukee, 
Wisconsin 

19-21 American Physical Society, Seattle, Washington 

24-28 American Society for Testing Materials, Atlantic City, New 
Jersey 

24-28 American Institute of Electrical Engineers, Swampscott, 
Massachusetts 

27-29 Institute of Radio Engineers, Boston, Massachusetts 
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Contributors to This Issue 


Paul L. Hartman is engaged in vacuum tube work with 
the Bell Telephone Laboratories in New York City. He 
graduated in electrical engineering from the University of 
Nevada in 1934, Dr. Hartman received his Ph.D. degree in 
physics at Cornell in 1938 and held an instructorship in 
that department the following year. 


Eugene C. Bingham is Professor of Research at Lafayette 
College, Easton, Pennsylvania. He received the A.B. and 
D.Sc. degrees from Middlebury College and the Ph.D. from 
Johns Hopkins and he has studied under Sir. J. J. Thomson 
and W. Nernst. 


Paul W. Kinney was employed by the National Adhe- 
sives Corporation in 1936-1937 and is at present carrying 
on research concerning adhesive materials for the Arm- 
strong Cork Company of Lancaster, Pennsylvania. He 
was awarded the Edward Hart Fellowship at Lafayette 
College in 1937 and received his M.S. in chemistry there in 
1938. 


C. J. Phillips is the Sales Manager of the Insulation 
Division, Corning Glass Works, Corning, New York, with 
which company he has been associated since 1931. He 
received his training at the University of Pittsburgh (B.S., 
1930), and at Oberlin College (A.M., 1931). Mr. Phillips 
is at present engaged in preparing a book on the technology 
and applications of glass. 


Pictures and biographical sketches of the following 
contributors have appeared in previous issues of the 


Journal of Applied Physics: 


F. Morgan, p. 391, Vol. 9, No. 6, 1938. 
M. Muskat, p. 391, Vol. 9, No. 6, 1938. 
Lloyd P. Smith, p. 415, Vol. 8, No. 6, 1937. 


William Fuller Brown, Jr., received his A.B. in English 
at Cornell. University in 1925 and his Ph.D. in physics at 
Columbia University in 1937. He was Lecturer in Physics 
at Columbia University from 1928 to 1938; at present 
he is Assistant Professor of Physics at Princeton University. 


C. G. Brennecke F. Breazeale 


Francis Breazeale has been employed 
Cornelius G. Brennecke is Assistant 
d with the American Enka Corporation, 
refessor of Electrical Engineering at the 
a has divided his time between rheological 
eceived his A.B. and E.E. degrees in 1926 
8, hives the University of North Carolina in 1935, 
‘niversity, from which he received his A in physics from the same 
h.D. degree in 1936. institution in 1936. 
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Recent Booklets 


Shop and Laboratory, publication of the Mico Instrument 
Company, Cambridge, Massachusetts, describes a ‘Model 
G Graton-Vanderwilt Polishing Machine,” in the feature 
article of the December, 1939, issue. The table of contents 
also includes articles on: “‘Ruling of Non-Uniform Dials,” 
“High-Intensity Illuminators,” ‘‘Microscope Counting 
Stage,” ‘‘Refractive Index Prism,’ ‘New Instrument 
Lubricants,” “‘Precision Scales,’ ‘Carriage Release for 
10-inch Measuring Engine,”’ and ““Toroid Winder News.” 

The Du Mont Oscillographer for December, 1939, 
January, 1940, carries the second and last part of the 
article, “Study of Phase Displacement in Electrical 
Circuits from Linearly Expanded Lissajous Figures,” by 
H. D. Brailsford, of the Underwriters’ Laboratories, Inc. 
The Oscillographer is published by the Allen B. Du Mont 
Laboratories, Iné., 2 Main Avenue, Passaic, New Jersey. 


“A Broadcast Frequency Monitor for the 20-Cycle 
Rule,” and “One Cycle per Second from the Inverse 
Feedback Oscillator,”’ are the two articles published in the 
January issue of the General Radio Experimenter, publica- 
tion of the General Radio Company, 30 State Street, 
Cambridge A, Massachusetts. 


The Frontier, published by the Research Foundation of 
Armour Institute of Technology, Chicago, Illinois, includes 
among the articles in the December issue the following: 
“Sodium in Industry,” “New Safeguards in War on 
Rust,”’ and “The Secret of the Lens.” 


The Laboratory devotes its feature article and its cover 
drawing for a recent issue to the subject of ““Mendeleeff— 
Visionary and Prophet.” The Laboratory is printed by the 
Fisher Scientific Company of Pittsburgh, Pennsylvania, 
“for those interested in keeping informed on the latest 
developments of laboratory apparatus and technique.” 


New Catalog 


A new set of catalog sheets listing products of the 
Shallcross Manufacturing Company, 10 Jackson Avenue, 
Collingdale, Pennsylvania, has been recently issued. It 
is the Company's series No. 146. The items listed have 
been classified as follows: resistance standards and affiliated 
instruments, Wheatstone bridges, test sets and special 
apparatus, and telephone transmission test equipment and 
kilovolt meters. 


Varied Line of Graph Papers 


Eugene Dietzgen Company, manufacturers of drafting 
and surveying supplies, Chicago, Illinois, announces a 
broad and comprehensive line of graph papers for mathe- 
matical, commercial and industrial use. 
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Innovations in Instruments 


Generator for Radio Testing 


For use in testing radio receivers in the ultra-high 
frequency range, the new Type 804-A U-H-F signal 
generator (sce accompanying figure) has been designed by 
the General Radio Company, 30 State Street, Cambridge, 
Massachusetts. The instrument produces carrier frequencies 
over a wide range, extending from 7.5 to 320 megacycles. 
The output voltage range is 10 microvolts to 20 millivolts. 
This signal generator is direct-reading in both frequency 
and output voltage and is provided with internal 400-cycle 
modulation up to 60 percent. External modulation can also 
be used. A voltage-regulated a.c. power supply is included. 


The oscillator circuit employed makes possible operation 
at frequencies as high as 330 megacycles. The tuning 
condenser is small and is driven by a worm and gear 
assembly. The coils are mounted on a Mycalex disk which 
can be rotated by means of contro! knob on the panel. The 
attenuator is capacitive and is provided with a compen- 
sating mechanism in order to present a constant capaci- 
tance to the oscillating circuit. An output cable is furnished. 


Electronic Devices Announced 


The Farnsworth Television and Radio Corporation of 
Fort Wayne, Indiana, announces the company’s extensive 
research activities on electronic devices, including television 
equipment and electron multiplier tubes of various types. 
Many of these are now being produced commercially. 

The action of the electron multiplier is based on the 
principle of secondary emission from a series of caesiated 
surfaces at successively higher potentials. The impact of 
each primary electron in each successive stage causes the 
emission of from 2 to 5 secondaries, thus amplifying the 
current in geometrical proportion to the number of stages. 
In view of the fact that electron multipliers give promise of 
performing many functions now performed by thermionic 
tubes, the producers of these multiplier tubes believe that 
they ultimately will have revolutionary effects on the 
science of radio communication and its many industrial 
and scientific applications. 
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Contributed Original Research 


Remarks on the Study of Structural Viscosity in Liquids, and the Derivation of 
Equations for the Flux of Such Liquids Through Tubes 


FRANCIS BREAZEALE 


The American Enka Corporation, Enka, North Carolina 


I 


HE aim of this paper is briefly to indicate 

how the elementary flow properties of 
certain imperfect viscous fluids may be found; 
and the manner by which may be mathematically 
developed the expression for the flux (through a 
tube) of these liquids, which possess (1) struc- 
tural viscosity alone, or (2) structural viscosity, 
changing to some thixotropic condition at high 
stress. It is believed that the expressions for flux 
of such liquids thus obtained are more accurate 
and more amenable to quantitative evaluation, 
over very considerable ranges of rates of shear, 
than any that have heretofore been presented." * 
Comparison of experimental and theoretical re- 
sults for the flux~pressure relationship is given 
for a typical example (cellulosic ‘‘solution’’). 

Flux~pressure curves obtained by forcing 
these fluids through capillaries (or pipes), fall 
into two general types, namely Bingham* and 
Ostwald.* Fig. 1 shows these two types. 

In so-called Ostwald flow, the curve is sigmoid 
in shape, first approaching an asymtote directed 
towards the pressure axis, and then rising to one 
directed towards the origin of the coordinates. 
In Bingham flow, the curve approaches an 
asymtote directed towards the P axis, and does 
not arise from it at the most elevated pressure 
considered. It thus appears that Ostwald flow is 
hut a special extension of Bingham flow. 


II 


The possible causes for the shape of either of 
‘hese two types of curves may be several; 
namely, 


Farrow, Lowe, and Neal, J. Textile Institute 19, 18 (1928). 
* RK. K. Schofield, Physics 4, 122 (1933). 


on C. Bingham and H. Green, J. Am. Soc. Testing Mat. 


* W. Ostwald, Kolloid Zeits. 47, 176 (1929). 
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1. Existence of a yield value. 

2. Creep, and not true flow. (This is also an indication 
of rigidity (elasticity).) 

3. (a) Structural viscosity—that is, lowering of viscosity 
with increase of shearing stress with very rapid recovery 
upon release of stress. 

(b) Thixotropy (slowly reversible change in liquid 
structure, manifested as a lowering of viscosity, from 
application of shear), Although the term is usually applied 
to sol-gel systems which show this property, we shall use 
it to describe the possible analogous condition in cellulose 
“solution.” 

(c) A hypethetical property which we shall call “pro- 
gressive thixotropy,”’ and define as a lowering of viscosity 
with shear; a process in all cases reversible, but whose 
time of reversibility is dependent on and increases with 
the magnitude of the disturbing shear. That is, at low 
shears, it acts as structural viscosity; at high shears as if 
a permanent liquid destruction had been effected. 

4. Slip at the tube walls, which factor might increase 
faster than does the stress. 


Items 1 and 2 are easily answered through use 
of a simple, chronophotographic method we have 
developed. A pigmented specimen of the liquid 
under study is forced through an untreated 
sample of the same liquid. From the. displace- 
ment ~time relationship of the resulting profile; 


Flux 


Pressure 


Pressure 
Bingham 


Ostwald 


Fic. 1, Flux-pressure curves obtained by forcing certain 


imperfect fluids through capillaries (pipes). 
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and from the form of the profile itself, it is 
readily determined whether yield value or creep 
occurs. 
Figure 2 is a sketch of the apparatus. The 
long tube with side arm, attached to glass jar 
with thick-walled tubing, is held in vertical 
position and the pigmented liquid from the jar 
allowed to ascend to the side arm. The main 
stopcock is now closed, the clear liquid poured 
down the tube and allowed to settle on the 
treated liquid. The side arm cock is now closed 
and the tube placed horizontally in the glass 


bath container, as shown. The tube has the same - 


(or nearly same) index of refraction as the liquid 
understudy, and the bath likewise is of this 
index. Hence distortion due to refraction is not 
troublesome. In our case, soft glass and xvlol 
were used for tube and bath. 

The clamp is used only when disassembling, 
to prevent the liquid from the jar from pouring 
out upon removal of the main glass observation 
tube. 

The main cock is now opened, and the profile 
between pigmented and unpigmented liquid pho- 
tographed at regular intervals. Three typical 
1, 25-sec. exposures at 30-sec. intervals, the first 
a few seconds after the cock was opened, (to allow 
for kinetic energy effects) are shown in Fig. 5. 

This record shows no appreciable yield value ; 
no creep (no rigidity or elasticity), as the profile 
is not blunt, and the displacement ~ time rela- 
tionship is linear. 

To investigate Item 3 a special Couette vis- 
cometer was constructed in our instrument shop, 
a photograph of which is shown (Fig. 6). It is 
‘featured by its ruggedness; speed range (10-300 


| 


Fic. 2. Sketch of apparatus used in the chronophoto- 
graphic method of determining existence of yield value and 
creep. 
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Rate @& Shear 


| Viscosity 


No 3.08 


Experimental 
formula 


33 

~& 


Rate of Shear (-sec~) 


Fic. 3. Experimental and calculated curves for the relation- 
ship between viscosity and rate of shear. 


r.p.m.) and small distance between cylinders in 
respect to their radii. This latter circumstance 
results in little change in rate of shear through 
the sample being investigated. Dimensions: inner 
cylinder-radius, 93.5 mm; height, 121 mm; outer 
cylinder-radius 99.5 mm; depth 106 mm. Clear- 
ance usually employed at bottom, 12 mm. 
Deflection~rotation speed curves ob- 
tained are transformed to viscosity~rate of 
shear curves by (1) assuming a perfect viscous 
fluid, and calculating rate of shear at inner 
cylinder from the familiar formula connecting 
this quantity and rotational speed,* and (2) 
dividing various points along the deflection ~ rate 
of shear curve by the corresponding rate of shear, 
and plotting the dividends as relative viscosity, 
in place of deflection. The determination of ab- 
solute values of viscosity, to replace relative ones 
is made from oil calibration curves. The assump- 
tion in the first step of a perfect viscous fluid does 
not defeat the purpose of the operation, inasmuch 
as the error thus involved can hardly be more 
than the change in rate of shear between inner 
and outer cylinders for a perfect viscous fluid, in 
our case, 11 percent. It will be seen that an 
error of this magnitude near the origin of the 
viscosity ~ rate of shear curve is allowable. It was 
found that there was no change of liquid struc- 
ture under stress, for all ranges of speeds; that 
is, the fluid was not thixotropic in this test. 


*Rate of shear at inner cylinder =—— 


X28. 
1/R,i—1/R, 


where: 
@=angular velocity of outer cylinder 


R,, R.=radii of inner and outer cylinders, respectively. 
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Fic. 4. Flux-pressure curves for capillaries of the same 
21 ratio, with the cellulosic “solution’’ used for 
Print 3. 


Bottom corrections were made by filling the 
instrument to two different heights; top cor- 
rections (piling up) are not important under 70 
r.p.m., as indicated by oil calibration curves. 
For this reason the speed range employed was 
from 10 to 60 r.p.m. 

The fourth point, concerning slip at walls, is 
probably best answered through measurement of 
the Ostwald asymptotes of two capillaries of 
equal length but of different radii. It has been 
found for solutions of the type under discussion, 
that at any given pressure on these asymptotes, 
flux (F)~r‘, within small experimental error. 
This is an indication of parabolic velocity profile ; 
and no slip.* We therefore say that in our case 
there is no slip in the lower region. 


Ill 


In this solution we are therefore led to believe 
that structural viscosity alone accounts for the 
lower region of both curves. Writing a general 
Viscosity ~rate of shear relationship from an 
examination of the experimental curve obtained 
‘s outlined above (Fig. 3) 


vhere 


= Viscosity at infinite rate of shear, 
A, B, constants 
d\ dr=rate of shear (negative in a capillary). 


* H. Lamb, Hydrodynamics (Cambridge University Press, 
th Edition, 1932), p. 586. 
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This condition, of course, says that true flow 
may never be reached, as the asymptote cannot 
pass through the origin (Bingham flow). Then 
how may be expressed the circumstance that in 
Ostwald flow, a condition of F~P is finally 
reached? Since we know from experiment at low 
rates of shear that the »~dV/dr function is 
immediately reversible, it seems necessary to 
postulate that above a certain critical stress, So, 
corresponding to a critical rate of shear, dV /dr., 
the action is not reversible. Although a gradual 
transition from reversibility to irreversibility (in 
respect to time—therefore ‘progressive thixo- 
tropy’’) would be much more reasonable, the 
discontinuity is postulated for mathematical 
simplicity. 

For Ostwald flow, then, 


A dr/dr=0 
n= [n+ 
B-d\ dr dV ‘dr=dV dre 


and above Sy» 
dV/dr=dV/dre 


The problem of determining Sy is now to be 
solved. Although this must be done to some 
degree intuitively, as will be apparent, it leads 
to results which are in fair accord with experi- 
ment. We find that in the Bingham assumption 
(no irreversibility) a term independent of pres- 
sure appears in the expression for the flux, and 
this term can easily be shown to be the ordinate 
distance between the Bingham asymptote and a 
parallel line through the origin (Ostwald asymp- 
tote for n=7n.). We may therefore suppose that 
when, in the Ostwald system, the curve leaves the 
Bingham axis, dV /dr, has been exceeded at the 
wall (where rate of shear is, of course, a maxi- 
mum), and that at this point, 8, So=pa/2l 
(where p=pressure at 8). Here we say that n has 
dropped to its lowest reversible value, 7.. We 
assume for simplicity, moreover, that the profile 
is parabolic here, although it is, of course, 
Bingham, (this assumption may be shown in 
most cases, a posteriori, to be a good one), and 
calculate ». by Poiseuille’s formula. 7, is, of 
course, solved exactly through Poiseuille’s for- 
mula, at any point on the Ostwald asymptote. 
At this point (assuming parabolic flow), rate of 
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shear at the wall=dV dr.=4F/ra*, (where made to some extent. Fig. 3 shows the experi 
F=flux, and a=radius of tube). A and B may mental »~dV'/dr curve to be about 50 sec~' 


: be solved through the simultaneous relations: (given by the viscometer) ; as well as »,, and n,. 
| ' The equation which appears to fit these values 
best is shown in dotted line. 

? ne Nx, AV dr., no being known from above. 

The values of A and B may have to be modified 

to some extent if the liquid exhibits an extremely 

rapid drop of viscosity with increased shear, at 

; low rates-of shear. In the case shown, which is 

ft an extreme one, this modification had to be 


Fic. 6. Photograph of viscometer. 


1. Driving mechanism powered by a }-hp Reeves varispeed motor. 
Two changes of large pulley and wheel, along with two changes of drive 
gears gives a continuous speed variation to the rotating cylinder from 
about 0.2 to 5 r.p.s. This mechanism is mounted upon a separate table 
from that which supports the instrument, in order to prevent vibration 


. from entering the viscometer. 

L 2. Housing for V belt which drives viscometer.’This belt connects 
. horizontal pulley on vertical shaft of driving mechz inism and one on 
7 vertical shaft which supports cylinder. 

: 3. Shaft which supports and drives outer cylinder. Cylinder is 
i removable. 

4, Outer cylinder. 

i 5. Inne~ cylinder, closed at bottom and top and weighted to sink in 
5 cellulose ‘‘solution”’ (sp. gr.: 1.1). 

my 6. Ball-bearing to keep inner cylinder cente red and/or supported. 
% This bearing is engaged or disengaged (inner cylinder hanging free) at 
a will. Its frictional resistance compared to the torque on the inner 
i. cylinder is quite ne gligible, however. 
‘a 7. Torsion wire suspending inner cylinder. 

ee 8. Bayonet locks fastening wire to shaft of inner cylinder and 
head at top. 


Centering head. 
10. Support for torsion wire, vertically adjustable only. 
- 180° scale graduated in 4° 
- Pointer attached to inner cylinder to measure deflection. 
is Support for inner cylinder assembly. It is removable from 
support (14) into which it slides, so that the whole assembly of inner 
cylinder, scale and suspension wire may be taken off together. 
Ruler is graduated in inches. 
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From static considerations we obtain the two 
expressions for flux (Bingham and Ostwald). 
For Bingham flow we have: 


A 
net —rr’P, 
B—dV/dr/ dr 


l=length of tube, P=pressure differ- 


where 


20a—L)(Q? 
4P 


f (= 


where D=A/n,+B; L=2D-—4B; A=P/2ln,z. 


ence between ends, r=differential radius. From 
whence we derive a first-order, second-degree 
differential equation, which is solved by means 
of the quadratic formula, the roots of which are 
found to be real and unequal. This gives V=f(r). 
By multiplying V by 2zr, and integrating be- 
tween r=0 and r=a, the Bingham expression for 
flux is obtained: 


2 r=a 
(2Q0r —L)(Q°r? —QLr+D?)!rdr 


4D?— 
160 


1 
(Q2r? —OLr+D?)}+0r—L/2 


—Da'*/6 is the constant term (independent of 


pressure) mentioned previously as the distance between Ostwald and Bingham asymptotes. 
It is of interest to note that the first term is exactly } of Poiseuille’s value for a constant viscosity 


=Na; 


In the case of Ostwald flow, 


condition that at r= 21S) 


P S,2 
= — 
4Aln, P? 


The results of these two are combined in one expression for Ostwald flow 


‘So? “|= 


rl’ So! So? 


P*n, Ps 


r=21S0/P 


40 r=0 


rPfa* 414S,* 
r= 
2in,L 4 | 


7(4D? 


The curves (Fig. 4) show F~P for capillaries 
of the same a 2/ ratio, (a=0.5 mm, and 0.75 mm; 
=12 mm and 18 mm) with the cellulosic 
“solution” used for Print 3. B is found to be at 
about the same pressure for both, as the theory 
vould demand (25 Ib./in.*). The only correction 
of any importance is the entrance end one. 
his was effected by first obtaining F~ P curves 
vith an orifice of capillary diameter, and then 
repeating the operation after attaching the 
capillary. The difference between the curves is 
robably as near the true F~P relationship as 
an be obtained with reasonable effort. The 
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(20r—L)(Q?r? — 
40°?  2OLSol 


— r=21S9/P 
In 
80) 


and that the third and fourth terms approach this value as P— ~. The last term—0 as P> ~. 


we perform a regular Poiseuille integration between r=a, and 
r=21S)/P; and then a Bingham one, as above, between r= 2/S)P 


and r=0, with the boundary 


) from Poiseuille development. 


(P>p) 


P 
OLr+D?)'rdr 


201S» 


2 


rdr. (2) 


of Ne Ne are lower 
curve (smaller capillary) as indicated previously. 
nx» is found to be 3.85 Poises; 7., 11.6 Poises. 
The Bingham (below 8) and Ostwald (above 8) 
curves are given by substitution of the values 
of A and B as found (and modified) in the 
preceding formulas Nos. 1 and 2. A=31,250; 
B=1010; dV /dr.= —3000 sec.-'. The curves 
plotted for both the large and small capillaries 
from these data agree with the experimental one 
rather well in the Bingham region; in the Ostwald 
region only in character. This is undoubtedly 
because there is no discontinuity between rever- 
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sibility and irreversibility in the viscosity ~rate 
of shear function at 8, but only the quickening of 
a transition (in terms of time) at this point. It 
will be noticed that the calculated Bingham axis 
for the larger capillary is in fair agreement with 
the drawn one (parallel to the calculated Ostwald 
asymptote, and tangent to the curve). This offers 
opportunity for further mathematical work. 

The graphical evaluation of the integrals in the 
expression for Flux (Bingham and Ostwald) is not 
difficult, for the F=f(r)-dr curves are single 
valued and fairly symmetrical from r=0 to r=a, 
and hence a few points, say at r=}, 3, 3, 4 a, 
together with the known points of f(r)dr=0, at 
r=0(0, and at r=a, are all that are necessary to 
draw an accurate curve. 


SUMMARY 


It is shown how the flux (through a tube) of 
certain fluids which show no yield value but do 
show structural viscosity, may be calculated with 
fair accuracy in the Bingham region. If we 


assume that the viscosity ~rate of shear relation 
is irreversible above S» (corresponding to 8 on 
the curve) the expression we derive for flux in 
this region agrees in character with the experi- 
mental curve obtained with the “‘solution”’ under 
discussion. The actual discrepancies between the 
two curves will, according to our theory, de- 
crease as the transition from reversibility to 
irreversibility becomes more rapid (magnitude of 
the “progressive thixotropy’’). In our case, the 
theory describes the transition as comparatively 
gradual. 


ACKNOWLEDGMENTS 


It is indeed a pleasure to acknowledge thanks 
to: Dr. A. J. L. Moritz, Vice President and 
Technical Manager, for permission to publish this 
paper; Mr. J. J. Schilthuis, Chief Chemist, and 
Dr. G. I. Thurmond, Research Director, for 
sympathetic advice and assistance in the course 
of the work; Mr. E. C. Jolliff, Microscopist, for 
assistance in photography. 


The Relation of Fluidity to Volume in Organic Liquids* 
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The Batschinski expression relating the fluidity of a compound to the specific volume at 


various temperatures is inadequate in the cases of associated compounds and those of high 
molecular weight. A modification, herein given, consists of the introduction of a correction 
term. The new expression, which relates linearly the fluidity and the corrected molecular free 
i volume, has universal applicability and it holds with a high degree of precision. The constants 
of the new relationship are related to the complexity of the compounds in homologous series. 
From the constants given, the fluidity of undetermined members of the homologous series 


INTRODUCTION 


HIS paper is based on a series of investi- 

gations of the relationship existing between 

fluidity and molecular volume in pure organic 

liquids. 

Using the valuable data obtained by Thorpe 

_* Presented at the Tenth Annual Meeting of the Society 

of Rheology, held at Pittsburgh, Pennsylvania, De- 
cember, 1938. 


** Present address: Armstrong Cork Company, Lan- 
caster, Pennsylvania. 
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studied can be calculated if the density is known at that temperature. 


and Rodger in their experimental investigations, 
combined with experimental data of later workers 
including the unpublished data of Bayer,' and the 
theoretical work of Batschinski, Coombs? in 1934 
under the direction of Professor E. C. Bingham 
began a thorough investigation of the relation 
between fluidity and volume. The work of 


1A. R. Bayer, Unpublished thesis, Lafayette College 
(1933). 

?C. E. Coombs, Unpublished thesis, Lafayette College 
(1935). 
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Coombs was continued by Grant,® Holderith,* 
Stookey,® Kinney,® and Adams.’ 

Thorpe and Rodger,* * measured the vis- 
cosities and densities of 87 pure liquids at 
temperatures, ranging from 0°C to their boiling 
points. Their studies included a number of com- 
pounds of each of several homologous series 
including the aliphatic and aromatic hydro- 
carbons, alkyl ethers, bromides, chlorides, iodides 
and sulphides and aliphatic acids, acid anhydrides, 
alcohols, esters, ketones and aldehydes. 

It was known that the fluidity increased with 
the temperature. Thorpe and Rodger found that 
the fluidities at a given temperature decreased 
successively as the molecular weight increased in 
a homologous series. This rule does not hold true 
in the case of the lower members of certain 
series, notably the aliphatic acids. Acetic acid is 
more fluid than formic acid; propionic acid is 
more fluid than acetic acid, and then butyric acid 
is less fluid than propionic acid, and the higher 
members behave in a normal fashion, that is, the 
fluidity decreases as the molecular weight 
increases. 

This anomaly is attributed to association, 
since in other physical chemical measurements 
the low members of the acid series behave not as 
if they are simple, single molecules but as if they 
are associated into multi-molecular groups. 

Batschinski' " indicated that the viscosity of 
a liquid is related to the specific volume by 
Eq. (1), 

(1) 


where 7 is the viscosity, v is the specific volume, 
i.e., the reciprocal of the density and c is the 
proportionality constant. The term w is defined 
as the limiting volume or the value which v 
approaches as the viscosity approaches infinity. 
—w) was named the “‘free volume.” Batschinski 


n=c/(v—w), 


E. Grant, Unpublished thesis, Lafayette College 
1936), 
‘A. E. Holderith, Unpublished thesis, Lafayette College 
1936). 


1938), 
* H. Adams, Unpublished thesis, Lafayette College (1939). 
*T. E. Thorpe and J. W. Rodger, Trans. Roy. Soc 
London) A185, II (1894). 
'T. E. Thorpe and J. W. Rodger, Trans. Roy. Soc. 
London) A189 (1897). 
'° A. Batschinski, Zeits. f. physik. Chemie 84, 643 (1913). 
"A, Batschinski, Physik. Zeits. 13, 1157 (1913). 
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himself, however, found that the relationship 
expressed by Eq. (1) did not hold in the case of 
associated substances, namely the lower acids 
and alcohols, water, etc., or even in nonassoci- 
ated liquids of high molecular weight. 

Equation (1) may be stated, 


o=K(V—-Q), (2) 


where ¢ is the fluidity, V is the molecular volume, 
K is the proportionality constant and Q is the 
molecular limiting volume. Eq. (2) states that 
the fluidity is proportional to the molecular free 
volume. 

When the fluidity at various temperatures is 
plotted against the molecular volume at the 
same temperature for the compounds of a 
homologous series, curves are obtained which are 
not linear at low fluidities but which instead 
appear to be hyperbolic. 

Figure 1 shows a series of curves obtained 
when the fluidity ¢ is thus plotted against the 
molecular volume V in the case of the homologous 
series of n-alkyl bromides. The bromide series is 
typical and it is here given as an example because 
it is the most complete with respect to the 
fluidity data. It will be noted that the curves are 
apparently linear in the lower members of the 
series but as the complexity increases the curves 
lose their linearity. Batschinski’s equation is 
apparently applicable to the lower members but 
not to the upper members. The same result has 
been obtained in the other series investigated. 


8 


Fiuidity 


Fic. 1. Fluidity vs. molecular volume; bromide series. 
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NeW RELATIONSHIP 


Bingham and Coombs have developed an 
equation which fits the data much more ac- 
curately than Batschinski’s equation. The new 
relationship is, 


B 
V=Ago——+C (3) 
where | and @ are, respectively, the molecular 
volume and fluidity, and A, B and C are con- 
stants of this hyperbolic equation. 
The constants A, B and C of Eq. (3) have been 
calculated by use of the three-point method 
using the following equations: 


Ve—V, B 
A= (4) 
1 1 
= )/( ). (5) 
C=—A¢,+B (6) 


Each compound has a value of A, B and C. 
These values are given in Table I. 

The hyperbolic equation (3) has been found to 
fit every curve. It holds for associated compounds 
and those of high molecular weight. The average 
percentage deviation between observed V and 
V’ calculated by use of the constants A, B and C 
for 50 compounds investigated is 0.05 percent, 
including associated compounds and those with 
chains of as high as 18 carbon atoms. As an 
example the observed molecular volumes and the 
molecular volumes calculated by use of Eq. (3) 
are given in Table II, for four members of the 
n-alkyl bromide series, representing 2, 6, 12 and 
18 carbon atoms, respectively. The percentage 
deviation between observed and calculated V for 
the various compounds and the average per- 
centage deviation for the various series are given 
in column 5 of Table I. 

The fact that the Bingham-Coombs hyperbolic 
equation (3) requires three constants for each 
compound where the Batschinski equation re- 
quired only two constants may at first appear 
cumbersome. However, it is shown later on in 
this paper how the three constants can be 
related to molecular complexity by rather simple 
relationships and thereby cut down the total 
number of constants needed in any homologous 
series of compounds. 
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The equation, 


V—2=A¢, (7) 


where A =1 K, is obtained from Eq. (2). Also by 
rearranging Eq. (3) it can be seen that this 
equation is very closely related to the rearranged 
Batschinski Eq. (7) and differs essentially only in 
the addition of a correction term to Eq. (7). 
Thus Eq. (8) is obtained, 


V—C=Aqo—B/¢. (8) 


Batschinski found his limiting volume 2 to be 
additive and he assigned values to several 
elements and to various molecular configurations. 
Bingham and Grant found the C term of Eq. (3) 
also to be additive and nearly identical to the 
Batschinski limiting volume 2. The C term can 
be considered to be a limiting volume. The 
elemental values of limiting volume as calculated 
by Batschinski from Eq. (7) and those as 
calculated from Eq. (3) are given in Table ITI. 
For purposes of comparing the two equations, 
the limiting volume 2 as obtained by Batschinski 
can be applied without loss of validity. Hence 2 
may be substituted for C in Eq. (8) to give, 


V—-2=Ao—B  ¢. (9) 


Eq. (9), which is the rearranged Bingham- 
Coombs equation, therefore differs from the 
Batschinski equation, as written in Eq. (7), only 
by the term —B ¢. 

In calculating C from additive 2, no values are 
needed for the double bond or for the iso grouping, 
in the case of the new hyperbolic relationship. 

The substitution of additive 2 for C in Eq. (9) 
is called restriction I. The average percentage 
deviation between observed and calculated values 
of V after restriction I is 0.11 percent for 98 
compounds. Since there is a slight difference 
between C and &, there is a corresponding differ- 
ence between the other two constants A and B 
of the hyperbolic equation. The values of A and 
B which were calculated using 2 in place of C are 
designated by the subscript 2. Thus A» and Be 
are the constants of Eq. (3) after restriction I. 
k-q. (3) becomes Eq. (10). 


V=Aod— B2/o4+2. (10) 


The substitution of additive @ for C has cut 
down the total number of constants needed in 
the calculation of molecular volume quite con- 
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siderably, 


since once the additive values of & 


have been evaluated, we do not need to calculate 


PABLE 


I. The values of the constants 


A, B, C, As, Bo, Q, As, Bz and By, for various compounds studies and the percentage 


C for every compound. The value of of Eq. (10) 


for any compound can be easily computed by 
use of the atomic limiting volumes in Table ITI. 


The values of As, Be and © for the 


various 


deviation between the observed molecular volumes and the molecular volumes calculated from the fluidity data by 


using the values of the constants in Eqs. (3), 


(10), (15) and (17), respectively. 


COMPOUND 


n- Hydrecarton Sevies 


Pentane 
llexane 
lleptane 
Octane 
Nonane 
Decane 
Undecane 
Dodecane 


n-Ether Series 
Diethyl ether 
Ethyl propyl ether 
Dipropyl ether 
Dibutyl ether 
Diamyl ether 


n-Mercaptan Series 
Ethyl mercaptan 
Propy! mercaptan 
Butyl mercaptan 

Amyl mercaptan 

Hlexyl mercaptan 
Heptyl mercaptan 
Octyl mercaptan 

Nonyl mercaptan 


Thiol Series 
Propane thiol- 
Pentane thiol-2 
Hexane thiol-2 
Hleptane thiol-2 
Octane thiol-2 
Nonane thiol-2 


n-Bromide Series 
Ethyl bromide 
Propyl bromide 
Butyl bromide 
\myl bromide 
llexyl bromide 
Hleptyl bromide 
Oetyl bromide 
Nonyl bromide 
Decyl bromide 
Dodecyl bromide 
lridecyl bromide 
letradecyl bromide 
Ventadecyl bromide 
\lexadecyl bromide 
lleptadecyl bromide 
Octadecyl bromide 


n-Acid Series 

rmic acid 
\cetic acid 
‘ropionic acid 
Butyric acid 
iproic ac id 
leptoic acid 
aprylic acid 
Laurie acid 
\ivristic acid 
Palmitic acid 
tearic acid 


n-Alcohol Series 

Methyl alcohol 

thy! alcohol 

Vropyl alcohol 

sutyl alcohol 

leptyl alcohol 

alcohol 
Oectadecyl alcohol 


0. 04420 | 
| 0.05505 | 
| 0.06067 | 
| 0.07236 | 
| 0.09275 
| 0.10130 
0.12337 | 
| 0.13184 


| 


| 
| 0.03 368 | 


0557 3| 
0.08399 | 
0.11357 | 


0.04134 
0.04196 | 
0.06100 
0.06390 
0.07625 | 
0.08875 
0.10763 
0.12745 


0.05607 
0.07378 
0.09228 | 
0.10342 | 


0.03753 

0.05090 
0.07032 

0.07490 | 
0.09713 

0.10631 

0.12319 

| 0.14584 

| 0.17309 

0.22520 

0.21233 
| 0.26662 
0.260360 
0.31434 
| 0.33739 
| 0.41086 


0.02011 
| 0.03710 
0.04428 | 
0.06701 | 
| 0.09049 | 
13560 | 


| 0.2 2001 


| 


| 0.01278 
| 0.02407 
10.0 2970 | 
| 0.04. 348 | 
0.08214 | 
0.13105 
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0.00 
93.65 
445.71) 
372.23 
28.18) 
204.53 
124.56 
234.14 


— 1,236.40) 
57 78.23 
200.16 
170.03) 


— 1,358.65) 
147.28 
—404.12 
109.24 
119.69 
137.31 
115.63 

97 44) 


431.04 
77.06 
46.66 
82.05 


340.67 
187.18 


35.89) 
90.32 
213.23 
56.26 
115.24 
51 02) 


121 00 


50.98) 


58} 158. 58 


C 


96.07 
113.67 
133.53 
151.09 
166.00 


185.94 | 


202.29 
222.43 
Aver. 
92.43 


126.43 


158.95 | 


193.74 
Aver. 


56.46 
80.74 


93.01 | 


114.60 
132.03 
149.92 
167.04 


184.00 | 


Aver 


117.12 


131.58 
148.42 
166.63 


Aver. 


67.74 


82.21 | 
96.71 | 


118.52 


130.99 | 


150.51 
168.15 
184.70 
201.45 
236.34 
256.91 
272.45 


293.05 | 


309.91 
327.56 
341.35 

Aver. 


37.27 
$5.27 
72.94 


88.37 | 


126.09 
140.91 


234.58 


382.44 
Aver. 


39.46 


57.40 


75.59 
92.15 
142.98 


Aver. 


0.01 
0.01 
0.02 
0.04 
0.02 
| 0.05 
| 0.05 
=0.03 


0.03 


0.05 
0.03 
| 0.02 
=0.03 


0.02 
0.06 
0.10 
0.03 
0.02 
0.04 
0.02 
0.03 


. =0.04 


0.06 
0.03 
0.03 
0.02 


0.01 
0.01 
0.12 
O11 
0.05 
0.04 
0.08 
| 0.05 
| 0.02 
0.06 
0.09 
0.09 
0.08 
0.08 
0.05 
0.18 
=0.07 
| 

0.02 
| 0.04 
| 0.07 
0.06 
| 0.03 


0.11 


=0.06 | 


0.01 
0.01 
0.13 


| 0.14 
| 0.14 | 


0.08 


0.04255 


0.05257 | 


0.06122 
0.07241 


| 0.08619 


| 


0.10034 
0.11557 


| 0.13392 


| 0.03672 


| 0.10575 
| 


0.04619 | 
0.05822 | 


0.08087 


| 0.03521 


0.04782 
0.05605 
0.06731 
0.07918 
0.09429 
0.11955 
0.12523 


0.05150 


| 0.10350 | 


0.063 18 


0.07538 | 


0.09000 


0.12176 


| 0.05077 


| 0.19861 | 


0.05804 
0.07047 


0.08492 | 
0.09161 | 


0.10747 
0.12425 
0.13913 
0.16154 


0.21950 


0.25358 | 


| 0.26514 


0.29491 
0.33968 


| 0.36586 | 


| 0.02104 


| 


| 
| 


| 


| 


0.03836 | 
| 0.04560 | 


0.05769 
0.09254 
0.10980 
0.11135 
0.24127 
0.25569 
| 0. 33531 
| 0.41927 


0.02107 
0.03662 
0.04798 
0.05815 
0.08641 | 
0.10248 
(0. 23811 


| 402.89] 
413.14] 
352.50 
286.75 
259.79 
231.70 
209.67 | 


| 


| 
= 


—81.92 


11.73 
102.32 


| 
| 
| 


—498.76 
} —6.31 
39.80 
64.80 
86.45 
89.82 
| 
| 
| —426.38 
— 128.02) 
21.81 
—54.60) 
102.52) 
102.39 
109.28) 
113.92 
127.46 
116.71| 
126.48 
139.47 
134.05 
126.85 
146.47) 
150.68! 
| 
25.32) 
53.36] 
186.37) 
140.91) 
106.37 
100.18 
110.48) 
79.41) 
79.26) 


87.25) 303.46 | 


49.68) 


Q 


97.72 
115.49 
133.24 
151.00 | 
168.76 
186.52 


204.28 | 


222.04 
Aver. 


89.61 | 


107.37 
125.13 
160.65 
196.17 
Aver. 


60.07 | 


77.83 
95.59 
113.35 
131.11 
148.87 
166.63 
184.39 
Aver. 


77.83 
113.35 
131.11 
148.87 
166.63 
184.39 

Aver. 


61.33 
79.09 
96.85 
114.61 
132.37 
150.13 
167.89 
185.65 
203.41 


238.93 | 
256.69 | 


274.45 
292.21 
309.97 
327.73 
345.49 

Aver. 


37.06 | 


54.82 
72.58 
90.34 
125.86 


143.62 | 


161.38 


232.42 | 


267.94 


338.98 


Aver. 


—88.78 
—84.20 
—50.79) 


36.33 | 
54.09 | 
71.85 | 


—5.22| 89.61 


56.85} 142.89 


43.78| 160.65 
78. 79| 302.73 | 0.07 
Aver. =0.30 | 


% 
| Dev 
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0.02 
| 0.02 
0.04 
0.06 
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| 0.05 
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0.06 
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0.06 
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0.03 
0.06 
0.16 
0.05 
=0.06 
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0.02 
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=0.04 
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0.11 
0.04 
| 0.16 

0.14 

0.22 
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| 
0.04 
| 0.04 


0.39 


as | 


04541 
0.05257 
0.06158 | 
| 0.07243 
0.085 12 
0.09966 
| 0.11604 
0.13427 | 
| 0.03672 
| 0.04666 | 
0.05727 | 
0.08049 | 
| 0.10637 


| 
| 0.04113 | 
0.04782 | 
| 0.05637 | 
0.06677 | 
0.07903 | 
0.09315 | 
0.10911 | 
0.12694 | 
| | 
| 0.04191 | 
| 0.06318 
0.07567 
0.08941 
0.10438 | 
0.12060 | 


0.05077 | 
0.05939 | 
0.06942 | 
0.08087 | 
0.09375 
0.10804 
0.12376 | 
0.14089 | 
(0.15944 | 
| 0.20081 | 
0.2 22363 | 
| 0.24786 | 
| 0.27351 
10. 30059 | 
| 0.32908 
| 0.35900 


| 
0.02104 
| 0.03154 
0.04437 | 
10.057 20 | 
08909 
0.10736 
0.12720 
| 0.22208 
| 0.27886 


7 | 0.34186 


0.41109 


0.02641 
| 0.03662 

0.04710 
0.05784 
| 0.09163 
| 0.10342 
| 0.23563 


—238.4 | 
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78.3 | 
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= 
= 


> 
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~= 


= 
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103.4 |(1.09)| 
—84.2 | 0.26 
0.53 | 


By Dev 
719.6| 0.33 
543.1} 0.31 
428.4) 0.05 
348.5| 0.02 
200.8) 0.14 
247.1| 0.05 
213.3) 0.14 
186.6} 0.13 

| 0.15 
709.5 | 0.18 
586.8 | 0.13 
502.5} 0.42 
393.4| 0.09 
325.4} 0.12 


—305.6| 0.48 
—145.5) 0.32 
—54.2|} 0.05 

5.7 | 0.10 
48.9} 0.03 
81.7 | 0.05 
107.5 | (1.35) 
128.6| 0.10 

0.16 

—1,329.0| 0.84 

—20.9) 0.08 

38.8 0.05 
72.6 0.09 
93.5| 0.04 
107.2) 0.17 
0.21 
—383.9| 0.23 
—118.7| 0.22 
—16.1} 0.12 
35.0| 0.99 
64.3! 0.42 
82.9! 0.15 
95.4| 0.07 
104.3 | 0.14 
110.9| 0.09 
119.8! 0.18 
122.9| 0.11 
125.5| O.11 
127.6| 0.21 
129.3 | 0.17 
130.7 | 0.16 
132.0| 0.22 
| 0.22 
373.1 |(10.12) 
220.2 | (3.72) 
161.8| 0.15 
130.0| 0.11 
95.5| 0.22 
84.9) 0.25 
76.7) 0.99 
56.4) 0.44 
50.1| 0.76 
45.3| 0.62 
41.4) 0.26 
| 0.42 

| 
—236.0 | (4.50) 
—103.4| 0.33 
—48.2| 0.56 
—16.4| 0.67 
31.6| 0.78 
40.8 | 0.39 
83.3) O.11 
0.47 


a 
A | B | @ | Dev] As Bs Dev.| 
| | ‘3 
308 57 0.20 
0.03 
0.04 | 
233.1 | 0.09 
236.8 | 0.06 
— 
594.92 
| 396.66 
395.98 
0.19 
| | | 
| | 
— 199.18 0.03 
ue 
9.39 
| 58.65 
91.63 
| 
| 
—1,254.6 | 0.79 
—29.9 | 0.08 
55.5 | 0.05 
47.3 | 0.06 
99.8 | 0.06 = 
— 70.5 | 0.08 
aa 
} | | 
—71.9 | 0.13 i 
285.74 —198.1 | 0.60 
5.00 155.7 | 
131.22 112.0 
129.16 88.4 | 0.11 
89.02 129.1 | 0.11 
57.37 108.8 | 0.10 
73.79 125.9 | 0.17 
119.39 144.0 | 0.12 : 
90.05 123.7 
157.20 173.6 || 
114.90 141.1 | 0.19 
123.59) 110.5 | 0.15 
79.97 124.6 | 0.19 
| 0.16 
| 
25 
— 105 92 
0.06 160 ( 
0.11 134 
0.03 88 
i— 0.17 157 64 
| — 22 6 
- | - 0.04 88 
190.20) 0.16 
109.67, 0.26 
92.50 0.46 
65 36 0.59 —8.6 | 0.6 
77.6 | 0.09 


compounds studied are given in Table I, together 
with the percentage deviation between observed 
molecular volumes and those calculated by use . 
of Eq. (10). In Table II are given four typical where (l’—2)+ Bz ¢ is the corrected free volume 


foregoing statement. 


F=A9, (11) 


members of the n-alkyl bromide series, showing ¥. The proportionality factor Az of Eq. (11) isa \ 
the observed and calculated molecular volume, Very important constant. It has been designated ‘ 
calculated from the fluidity at various tempera- by Stookey as the “rheolar free volume” and c 
tures, after restriction I. is theoretically defined as the free volume 
necessary to give unit fluidity to a liquid. Eq. 
Stupy or CONSTANTS A AND B (12) may further clarify the rheolar free volume ) 
The constants Az and By of Eq. (10) were 
investigated to determine how they varied from Aue 8/4 (12) | 
compound to compound in a homologous series. A» has a different value for each compound. 
k-q. (7) and Eq. (9) both show that the fluidity is The value of Az» for the various members of a | 
proportional to the free volume. In the first case homologous series becomes larger as we increase ! 
it has been shown that the Batschinski expression — the molecular complexity of the compounds of the ‘ 
is only an approximation and that in reality the — series. That is, in the larger molecules more free { 
fluidity is proportional to a “corrected free volume is needed to give the same fluidity. Fig. 2 ( 
volume.”’ This corrected free volume expression, gives the curve obtained when the values of A, \ 
since it holds true in every case, may be con- are plotted against the number of carbon atoms i 


sidered a fundamental property and may be in the case of the n-alkyl bromides. Similar 
designated by the term F. Eq. (11) expresses the curves are obtained in the cases of the other 


TABLE I (Continued). 


| & % ; A 
COMPOUND A B ( Dev. A: By 2 Dev. A; By Dev.| By Dev. 
Ester Series | | 
Methyl formate 0.04400 —916.35 45.57 0.02802 355.6 54.82 | 0.06 | 0.02802 355.7 0.06 1,068.0 | (4.30) 
Ethyl! formate 0.04552 —256.56 67.93 0.03697 356.5 72.58 | 0.05 | 0.03663 332.0 | 0.04 787.3; (2.20) 
n-propyl! formate 0.05037 74.13) 87.81 0.04544 365.7 90.34 0.06 0.04584 391.9 0.09 634.0 | (1.09) 
n-butyl formate 0.06117 86.62) 107.51 0.05963 131.9 108.10 | 0.01 | 0.05564 20.7 | 0.26 535.9 | (2.86) 
n-butyl acetate 0.06987 167.73) 123.58 0.06448 363.1 125.86 | 0.10 | 0.06602 | 428.1) 0.19 | 467.2; 0.28 
n-butyl n-propionate 0.08367 174.26 140.06 ~ | 0.07404 439.6 143.62 0.13 0.07700 535.0 | 0.31 416.0 0.38 
Ethyl n-valerate 0.07461 430.05 143.62 | 0.15 | 0.07700 507.1 | 0.23 4160) 0.34 
n-butyl n-butyrate 0.09595 186.88 157.45 0.08343 431.22, 161.38 | 0.17 | 0.08857 552.1 0.40 376.3 0.64 
n-butyl »-valerate 0.11272 137.38) 174.97 0.09699 353.8 | 179.14) 0.15 | 0.10073 | 417.7 | 0.25 344.3 0.32 
Ethyl n-heptylate 0.09578 361.99) 179.14 | 0.18 | 0.10073 444.6 | 0.33 344.3. OA3 
n-hepty! acetate 0.10733 286.98 179.14 0.22 | 0.10073) 207.2 | 0.42 344.3 0.74 
n-amyl n-butyrate 0.10878 178.21) 176.06 0.09800 343.79 179.14 0.14 0.10073 392.2) 0.20 344.3 0.23 
n-hexadecyl acetate 0.36388 90.94) 332.12 0.26169 174.09 338.98 | 0.18 | 0.23674 135.9 | 0.38 204.4 0.63 


Aver. =0.12 0.24 0.44 


Miscellaneous Hydrocarbons 
‘ Benzene 0.03260 404.88) 86.47 0.05016 —200.70| 79.80 | 0.03 


Toluene 0.05337 —59.18) 96.90 0.05182 9.88) 97.56 | 0.05 | - 

o-xylene 0.06010 | 101.03) 114.40 - 0.05767 175.71; 115.32 | 0.06 | - 
m-xylene 0.05989 —32.63) 112.79 | - | 0.05456 233.07) 115.32 | 0.07 | - . - 
p-xylene 0.05615 142.06; 115.11 | 0.05586 176.43) 115.32 | 0.03 | = 
isoprene 0.04514 | —1,537.37| 77.17 | | 0.04247 —855.12) 79.88 | 0.02 - 
Dially! 0.04261 161.41, 104.09 - |0.05108  —1,029.14) 97.64| 0.07 - - 
Isopentane 0.04831 —661.14) 93.64 — | 0.04308 278.5 97.72 | 0.04 - ~ 
Isohexane 0.04715 | 954.48 118.01 0.05106 286.42) 115.49 | 0.52 — 
Isoheptane 0.06014 | 432.39 133.08 0.05724 498.74 133.24 | 0.65 - - - 


Miscellaneous Halides } | 
Allyl chloride 0.03765 114.75) 70.80 0.04297 —351.69| 67.61 | 0.03 | _ 
Propyl chloride 0.04500 | —33.10 75.59 - 0.04349 | 110.22 76.53 | 0.02 - — 


Isopropyl! chloride 0.05188 —714.39, 73.18 0.04618 76.53 | 0.02 
Isobuty! chloride 0.05430 —64.15) 93.69 | 0.05307 28.25) 94.29 | 0.08 -- 
Carbon tetrachloride 0.06702 115.17) 90.76} 0.09171 —305.03) 84.00 | 0.27 — 
Isopropyl bromide 0.04988 245.62, 85.02 0.06327 —391.0 79.09 | 0.06 

Isobutyi bromide 0.06484 171.43) 103.19 0.08046 —423.04 96.85 | 0.24 - - 
Methyl iodide 0.04224 —84.67| 53.20 0.04037 —9.40) 53.96 | 0.01 | — 


Ethyl iodide | 0.05123 112.37| 72.49 | 0.05319 41.18) 71.72 | 0.02 

Propyl iodide 0.06408 | 102.66, 89.45 - | 0.06460 | 115.58, 89.48 | 0.07 — ‘ om 

Isopropyl iodide 0.06726 35.54, 90.17 — | 0.06902 —22.64) 89.48 | 0.04 - 

Isobuty! iodide 9.07142 230.53, 108.62 0.07502 | 108.50, 107.24 | 0.06 | - | - 

Acetic anhydride | 0.05799 124.67; 89.75 — |0.05418 218.11, 91.07 | 0.13 - 

Propionic anhydride 0.07654 183.92) 123.84 | — | 0.06767 358.81) 126.59 | 0.20 — 
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series investigated. The curve of A vs. 1 is found 
io be parabolic, obeying the equation, 


(13) 


where » represents the number of methylene 
eroups in the molecule. The validity of Eq. (13) 
can be shown by taking the first derivative of 
Iq. (13) to give, 


dA»s dn=b+2cn. (14) 


In Fig. 2, AAe/An is plotted against n. The 
resulting straight line proves the validity of 
leq. 14. 

The various terms of the parabolic equation all 
have theoretical justification. The ‘‘a’’ term is 
not dependent upon the number of carbon atoms 
and is the free volume which must be supplied by 
the residue, i.e., the He, HBr, HS, etc.; the ‘‘d” 
of the ‘‘bn’”’ term must be the rheolar free 
volume value of a single CH. group and therefore 
in a chain, the value of a single CH» group is 


multiplied by the number of carbon atoms. 
These two terms alone would satisfy the straight 
line relationship A2=a+dn, which is what would 
be expected when a number of CH: groups and a 
residue are involved. There is still another term, 
the ‘‘cn*”’ term for which we must account. The 
fact that the length of the chain is squared 
suggests that an area is involved. If the chain 
may be considered to be rotating in propellor 
fashion, then the effective volume, or volume of 
influence, of the molecule is not alone that of the 
actual molecule itself. The influence of the area 
of rotation must be taken into account. 

We find an illustrative parallel to this n- 
squared relationship in the work of Staudinger™ 
on the viscosity of solutions of long chain 
molecules. He finds that the range of action in 
solution of the thread molecules does not increase 
Staudinger, Die 


hochmolekiilaren organischen 


Verbindungen Kautschuk und Cellulose (J. Springer, Berlin, 
1932). 


137.78 


234.80 


Hexyl Bromide 


Dodecyl Bromide 


Octadecyl Bromide 


| V (oss.) | | 
Wr.) V (CALC.) V (cALc.) V (CALC.) V (CALC.) 
- — ZO. % AFTER % AFTER % AFTER % 
| SP. GR (3) Dev. Rest. I DEv. Rest. Il DEv. Rest. III Dev. 
| 
Ethyl Bromide 
0 | 209.2 73.96 73.96 | 0.00 73.99 0.04 73.99 0.04 73.79 0.23 
10 | 231.7 74.96 74.97 | 0.01 | 74.93 0.04 74.93 0.04 74.75 0.28 
20 255.1 | 76.00 75.97 0.04 | 75.95 0.07 75.95 0.07 75.78 0.29 
280.1 77.07 77.07 76.92 


3.4 | 0.01 137.69 0.07 137.13 0.48 138.37 0.42 

20 | 99.1 140.43 140.37 0.04 140.42 0.01 140.07 0.26 141.01 - 0.41 
40 | 128.4 143.27 143.27 0.00 143.33 0.04 143.20 0.05 143.91 0.45 
60 161.0 146.45 146.44 0.01 146.43 0.02 146.49 0.03 147.06 0.42 
80 | 195.5 149.85 149.85 0.00 149.76 0.06 149.90 0.03 150.37 0.35 
100 | = 231.5 153.01 153.37 - 0.24 153.14 0.08 153.40 0.25 153.79 0.51 
Aver. =0.05 0.43 


0.23 234.13 0.51 233.55 0.76 233.95 0.59 

20 26.2 239.35 239.42 0.03 239.68 0.14 239.38 0.01 239.62 0.11 

40 40.4 243.61 243.61 0.00 244.06 0.18 243.92 0.13 244.07 0.19 

60 | 57.1 247.97 247.91 0.02 248.23 0.10 248.20 0.10 248.30 0.13 

80 | 76.4 252.58 252.58 0.00 252.57 0.00 252.62 0.02 252.70 0.05 

100 | 97.0 257.14 257.42 0.11 257.00 0.05 257.11 0.01 257.17 0.01 
Aver. =0.06 0.16 : 


40 | 13.0 342.72 340.54 0.64 338.66 1.18 340.58 0.62 340.01 0.79 
60 | 26.6 349.28 349.27 0.00 349.55 0.08 350.36 0.31 350.08 0.23 
80 37.5 355.15 354.63 0.15 355.19 0.01 355.63 0.14 355.43 0.08 
100 50.7 361.21 360.60 0.17 361.25 0.01 361.23 0.01 361.09 0.03 
120 67.2 367.77 367.77 0.00 367.84 0.02 367.76 0.00 367.65 0.03 
0. 374.50 : 
Jl 
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| 
lance Il, Example of calculation of molecular volumes at various temperatures from the fluidities at the same temperatures S| 
by use of the Bingham-Coombs equation and after restrictions I, II and III, respectively. : 
4 
| | | Aver. =0.01 | | 0.05 | | 0.23 
197 


proportionally to the length of the molecule but 
to the square of the length. In the case of very 
long molecules the effect of the “a” and “‘bn”’ 
terms of our parabolic equation would be negli- 
gible with respect to the effect of the ‘‘cn®” 
term. At the limiting viscosity the volume of the 
solvent available is just equal to the effective 
volume of the dissolved thread molecules, so 
that the molecules no longer possess freedom of 
movement. In our considerations a rough parallel 
might be drawn between Staudinger’s volume of 
solvent and our free volume ; between his volume 
of limiting viscosity, i.e., the volume of solvent 
plus volume of molecules and our molecular 
volume; between his volume of molecules and 
our limiting volume; and between his range of 
action in solution and our effective volume. The 
comparison cannot be carried beyond the quali- 
tative The fact that the n 
squared relationship occurs in both realms is the 


stage, however. 
significant relationship. 

-The values of a, 6 and ¢ of parabolic equation 
(13) for the various series studied are given in 
Table IV. The values of A calculated by use of 
Iq. (13) were given the subscript 3. Restriction 


TABLE IIL. Values of limiting volume from Batschinski and 
from Bingham-Coombs equation. 


ELEMENT 
OR BATSCHINSKI 
GROUP 


GRANT'S 
BINGHAM-C OOMBS 
VALUES 


8.84 
4.46 
17.76 
9.65 
15.63 
18.79 
21.35 
31.74 
0 


ge 


CH, 
S 
Cl 
Br 
I 
Double 
Bond 
iso- 0 


we — 


TA 


SERIES 


Bromides 
Hydrocarbons 
Mercaptans 
Thiols 
Alcohols 
Esters 

Acids 


Ethers 
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0.03780 
0.03727 
0.03331 
0.01930 
0.01645 
0.01256 
0.01210 
0.00361 


b 


0.00506 


— 0.00298 
0.00205 
0.00567 
0.00983 
0.00714 
0.00816 
0.00694 


0.0 


Fic. 2. Parabolic relationship between A» and » for 
bromide series. Also linear relationship between AA2/An 
and n. 


II comprises the substitution of Ag, in. the 
hyperbolic equation (3) to give, 


(15) 


The average percentage deviation between ob- 
served values of volume and those calculated 
using Kq. (15) is 0.20 percent for 72 compounds. 
The values of A; and By; for the various com- 
pounds studied are given in Table I, together 
with the percentage deviation between observed 
molecular volumes and those calculated by use 
of Eq. (15). In Table I] are given four typical 
members of the n-alkyl bromide series, showing 
the observed and calculated molecular volume, 
calculated from the fluidity at various tempera- 
tures, after restriction IT. 


BLE IV. Values of constants of Eq. (13) and of Eq. (16). 


0.00071 
0.00092 
0.00093 
0.00062 
0.60013 
0.00030 
0.00078 
0.00033 


3.9350 
— 3.0399 
— 3.8789 
— 2.5866 
3.2554 
2.5718 
3.3630 


— 2.3551 
—0.6071 
—0.7529 
— 0.7607 
—0.8505 
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Fic. 3. Interrelations of the functions of # and V. ‘‘a’— 
Batschinski region. ‘‘b’’—Hyperbolic region. 


Consideration of the significance of the B/¢ 
term and the B; constant now follows. The B/¢@ 
term is the difference between the Batschinski 
equation and the Bingham-Coombs equation and 
is therefore the most significant part of these 
researches involving the new relationship. The 
difference between the two equations is better 
understood by the graphical illustration given in 
lig. 3. The straight line corresponding to the 
equation 1 —Q=A¢@ is the Batschinski relation- 
ship. The hyperbolic equation is defined by the 
equation V—Q=A@—B/¢. In other words the 
B® term indicates, and is a measure of, the 
departure of the hyperbola from the straight 
line. It is possible to make the degree of deviation 
more clear by plotting ¢@ against B/¢. It will be 
noted that as @ increases B/@ decreases and the 
hyperbola approaches the straight line. In this 
fact lies the reason why Batschinski made his 
relationship one of a straight line. He worked in 
the region of high fluidity. The fluidity in the 
case of the lower members of a given series is, at 
the temperature of measurement (0°C to the 
boiling point), rather high and therefore in such 
compounds as pentane, hexane, diethyl ether, 
etc. Batschinski did obtain straight lines. These 
straight lines are in reality the asymptotes of the 
‘rue hyperbolas. 

As shown by Volkman and Petrokubi," who 
studied the fluidity of pentane at temperatures 


"R. G. Volkman and J. L. Petrokubi, Unpublished 
hesis, Lafayette College (1936). 
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as low as —70°C, the straight line obtained by 
Batschinski becomes the characteristic hyperbola 
at low fluidities. This observation emphasizes the 
importance of the magnitude of the fluidity at 
which the comparisons are made. 

The need for the correction term B/¢ has been 
attributed to the effect of intermolecular attrac- 
tion. Nearly all the values of B/@ are positive. 
Certain compounds, especially those of high 
fluidity, have negative values of B. Inaccuracies 
in application of the three-point method for the 
calculation of B prebably accounts for these 
negative values. 

When the values of B3, calculated after 
restriction II, are plotted against the number of 
carbon atoms, curves are obtained similar to 
those shown in Fig. 4. These curves appear to be 
hyperbolic and obey roughly the equation 


B;=E+dn’/, (16) 


where E is the horizontal asymptote and d and f 
are constants. The sign of the dn’ term is 
positive for curves of which that of the hydro- 
carbons is typical and negative for curves of 
which the curve of the bromides is typical. 

The curves obtained above appear, at first 
consideration, to be poor but it must be pointed 
out that the B;/¢ term of the hyperbolic equation 
is never more than 3 percent of the molecular 
volume V and is usually much less than 3 
percent, so that quite a variation in the value of 
B; will not greatly effect the value of V. 


Fic. 4. Hyperbolic relationship between B; and n. The 
open-circle curve is for hydrocarbons and the closed-circle 
curve for bromides. 
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prises the substitution of B, into the Bingham- 
Coombs hyperbolic equation to give equation 


E+dn/‘ 


V=(a+bn+cn’)o-—- +2. (17) 
\ 


The average percentage deviation between calcu- 
lated values of V after restriction III is 0.2 
percent for unassociated series. Association ap- 
pears to affect the calculated values of B; 
although due to the possibilities for accumulated 
error in B; the degree of association cannot be 
calculated. For the “associated” series, if we do 
not consider the highly associated formic acid, 
acetic acid, methanol and the lower esters, we 
have an average percentage deviation of 0.4 
percent. 

The values of B, for the various compounds 
studied are given in Table I, together with the 
percentage deviation between observed molecular 


is 18 volumes and those calculated by use of Eq. (17). 

rn In Table II are given four typical members of the 

Fic. 5. Log F vs, n-hydrocarbons, n-alkyl bromide series, showing the observed and 

calculated molecular volume, calculated from the 

The values of B; calculated by use of Eq. (16) fluidity at various temperatures, after restric- 
are given the subscript 4. Restriction III com- tion III. 


nw 


NUMBER ie 20°C 40°C | 60°C | 80°C 100°C 


C ATOMS Oss. Case. Oss. CaALc. Oss. | Ops. CALc. Ons. Cate. | OBs. CALc. 


5 15.13 15.17 | 18.32 18.32 | — — | 

6 13.32 13.32 83 | 

7 11.83 11.69 | 14.93 14.86 | 18.35 18.32 22.14 22.20 
8 

9 


a 
a 
w 


26.37 26.233 | — 
| 24.95 25.08 | 29.57 29.67 


10.16 10.26 | 13.44 13.38 | 16.91 16.92 | 20.74 20.91 


9.09 9.00 12.15 12.06 | 15.58 15.63 | 19.68 19.69 | 23.96 23.96 | 28.86 28.76 
10 7.95 7.90 | 10.88 10.85 | 1442 14.44 | 1846 18.55 | 22.81 22.91 | 27.89 27.87 
11 | 6.91 6.93 | 9.73 9.78 | 13.36 13.34 | 16.90 17.47 | 21.96 21.90 | 27.00 26.99 
12 | 6.01 6.09 8.76 8.81 | 12.32 12.33 | 16.54 16.45 | 20.96 20.93 | 26.14 26.17 
| 


Tance VI. Corrected free volume F at various temperatures, observed and calculated for aliphatic n-acids, along with calculated 
apparent associations of lower members. 


60°C | 100°C 
No. OF RaT1o0 RATIO Ratio | AVERAGE 
C ATOMS OBs. Calc. Carc./OBs Ops. CALc., Catc./ Os. Oss. Catc./OBS., ASSOCIATION 
1 2.39 7.99 3.34 | 3.09 10.17 3.19 3.85 12.02 2.12. | 3.22 
2 5.50 7.52 1.37 | 6.86 9.74 1.42 8.31 11.74 1.41 1.40 
3 6.85 7.08 1.03 8.42 9.32 1.11 9.99 11.47 1.15 1.10 
4 6.74 6.66 0.99 | 8.61 8.92 1.04 | 10.66 11.20 1.05 1.03 
6 6.18 5.90 0.96 8.31 8.18 0.98 10.61 10.67 1.01 | 0.98 
7 5.55 S30 1.00 7.87 7.83 1.00 10.73 10.43 0.97 0.99 
8 4.99 5.23 1.05 7.35 7.50 1.02 9.95 10.18 1.02 1.03 
12 4.16 4.10 0.99 6.62 6.31 0.95 10.22 9.26 0.91 0.95 
14 | 3.58 3.63 1.01 5.36 5.78 1.08 8.38 8.83 1.05 1.05 
16 — — — 5.44 5.30 0.97 8.04 8.42 1.05 | 1.01 
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CALCULATION OF FLUIDITY FROM 
VoLUME DATA 


If it is desired to calculate the fluidity when 
the volume is known, this can be done by use of 
the equation 


V-2\?2\3 
—)). as 
2A; A; 2A; 
A value of this research, aside from theoretical 
interest, is that we can calculate very accurately 
the fluidity of any member of the series inves- 
tigated at any temperature if we know the 
density of the material at that temperature. The 
normal straight-chain series of aliphatic hydro- 
carbons, ethers, bromides, acids, alcohols, esters, 
mercaptans and thiols have been investigated. 
The values of A; and By, for any member of 
these series can be calculated by use of Eqs. (13) 
and (16), using constants given in Table IV. 
The value of Q for the member can be obtained 
by using the additive values of limiting volume 
given in Table ITT. 


FREE VOLUME VERSUS CHAIN LENGTH 


The results of our investigations lead us to 
believe that the corrected free volume is a fun- 
damental property and therefore should have a 
definite relation to the complexity of the mem- 
hers of a homologous series. The corrected free 
volume F has been found to be related exponen- 
tially to the number of carbon atoms in the 
various series at any chosen temperature. This 
relationship is expressed 


F=De-*", (18) 


where D and K are constants. 

Figure 5 shows the linear curves obtained 
vhen the logarithm of the corrected free volume 
is plotted against the number of carbon atoms 
it the various temperatures for the hydrocarbon 
series. Similar sets of curves are obtained in the 
ases of bromides, ethers, thiols and mercaptans 
eries. Table V gives the values of F observed 
nd calculated by Eq. (18). The “observed” 
alues are (V—)+B/¢. 

In the case of series in which the lower mem- 
ers are associated, the results shown in Fig. 6 
re obtained. This acid series may be regarded 
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Fic. 6. Log F vs. n-acids. 


as typical. The deviation of the lower members 
from the straight line drawn through the upper 
members of the series gives us a method for the 
calculation of the degree of association. Thus it 
has been found that the ratio of the calculated 
free volume to the observed free volume is equal 
to the apparent molecular association. The values 
of observed and calculated F are given in Table 
VI, along with the degree of association. The 
associations calculated in this manner compare 
favorably with those calculated by other workers. 
Similar curves were obtained for the esters series 
but no good curves could be obtained for the 
alcohols. 

The authors wish to acknowledge indebtedness 
to Coombs, Grant, Holderith and Stookey for 
liberal use of their calculations and results. 


SUMMARY 


1. A new relationship between fluidity and 
volume is given. In it a correction term is applied 
to the Batschinski equation which states that 
the fluidity is proportional to the free volume. 
The Batschinski equation is limited in its appli- 
cation. The new equation is universal. 
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2. The three constants necessary in the new 
equation are studied. C is found to be additive 
and closely related to the Batschinski limiting 
volume. The limiting volume is regarded as the 
volume of the atoms themselves. 

3. The other two constants of the new equa- 
tion are related to the number of carbon atoms 
of the members of homologous series. 

4. The need for the B’@ correction term is 


attributed to intermolecular attraction at low 
temperatures (and low fluidities). 

5. The corrected free volume is a fundamental 
property. The fluidity is proportional to the 
corrected free volume. 

6. In a given series and at any chosen tem- 
perature, the corrected free volume is related 
exponentially to the number of carbon atoms in 
the chain. 


Conductivity Measurements on Potassium Halides 
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The electrical conductivities of flat crystal plates of KBr and KCI were measured over a 


temperature range of from 25°C to 600°C, at stresses of from 0.05 to 0.90 kv/cm. Ohm's law 
was found to hold, when cold conductivities were measured after preliminary heating. Current- 
time data indicated an exponential decay of current with time at the lower temperatures, and 
a current constant with time at the higher temperatures. An exponential relation between con- 
. ductivity and temperature was found to hold for both KBr and KCl, but in the latter case 
an abrupt change in the slope of the curve was found to take place at a temperature between 


200°C and 250°C. 


INTRODUCTION 


HILE much work has been done on the 

electrical characteristics of dielectric ma- 
terials, the complexity of the subject seems to 
justify additional research, since only sufficient 
data will make it possible finally to correlate all 
the phenomena involved with a relatively simple 
theory. 

In simple ionic crystals such as the alkali 
halides, at normal temperatures and low field 
strengths, metallic conduction is absent because 
of the closed configuration of the ions.' The 
current observed under such conditions is now 
ascribed by many!“ to a “‘lattice conductivity” ; 
that is, to the presence of points in the crystal 
lattice from which the ions are missing, and to a 
resultant ionic migration under the influence of 
the field. The current furthermore seems to be a 
linear function of the field i.e., to obey Ohm's 


1N. F. Mott, Trans. Faraday Soc. 34, 822 (1938). 

?W. Schottky, Zeits. f. physik. Chemie 29, 335 (1935). 

* A. Smekal, Physik. Zeits. 37, 150 (1936). 

*E. Koch and C. Wagrer, Zeits. f. physik. Chemie 38, 
295 (1937). 

®’ A. Wenderowitsch and R. Drisina, Zeits. f. Physik 98, 
108 (1935). 
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law. When, however, still at normal tempera- 
tures, potential differences approaching the order 
of magnitude of the breakdown potential are 
applied across the crystals, Ohm’s law no longer 
holds,® and the relation between current and 
voltage becomes exponential, somewhat in the 
manner first predicted by Poole.’ The equation 
proposed by Poole is 


k=ae"* 


(1) 
where «x is the conductivity, x the field strength, 
and a and b arbitrary constants, but it gives only 
a first approximation to experimental data. On 
the basis of much work’ in this field, von Hippel 
concludes that the phenomena in the region of 
breakdown are chiefly electronic in nature: 
Electrons locally bound to the lattice points 
become first excited and then conducting, moving 
in preferential lattice directions. If sufficiently 
accelerated by the field, these electrons ionize 
by impact. The ionization becomes cumulative, 

®F. Quittner, Physik. Zeits. Sowjetunion 11, 359 (1937). 

7H.H. Poole, Phil. Mag. 32, 112 (1916); 34, 195 (1917); 
42, 488 (1921). 


* A. von Hippel, Zeits. f. Physik 75, 145 (1932); J. App. 
Phys. 8, 815 (1937); Phys. Rev. 54, 1096 (1938). 
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the resulting space charges distort the field, the 
material is weakened, and breakdown occurs. 
This theory, of course, denies the older one of 
Joffé® and his co-workers, who maintained that 
ions loosened from the lattice produce the impact 
ionization which leads to breakdown. 

Again, when at low field strengths the tem- 
perature of the crystal is raised, the conductivity 
rises sharply. The relation 


(2) 


where «x is the conductivity, T the Kelvin tem- 
perature, and A and B constants, was proposed 
by Rasch and Hinrichsen'® to describe this 


*A.F. Joffé, Tae Physics of Crystals (McGraw-Hill Book 
Co., New York, 1928). 


1” A. Rasch and H. Hinrichsen, Zeits. f. Elektrochem. 
14, 41 (1908). 
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phenomenon, but it has been found to hold 


exactly, over the entire temperature range, for 


very few of the solid salts which have been inves- 
tigated. The whole question is complicated 
further by the presencé of the so-called anomalous 
charging current first noted by the Curies'! in 
1888, and variously explained as polarization, a 
space-charge effect, or the like. The current, in 
most cases observed,"-" appears to follow an 
empirical law of the form 


I=kt-™ (3) 


"J. and P. Curie, Ann. chim. phys. 17, 385; 18, 203 
(1889). 

2A, D. Goldhammer, Physik. Zeits. Sowjetunion 5, 
523 (1934). 

13 E. Hippauf and R. Stein, Physik. Zeits. 39, 90 (1938). 
4 F, Seidl, Physik. Zeits. 39, 714 (1938). 
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where ¢ is the time in seconds and k and m are 
constants. 

In the present work the current passed by 
pure crystal plates of KBr and KC] was measured 
over a wide range of temperatures. At each tem- 
perature a number of rather low voltages (well 
below breakdown potential) were applied, and 
current-time data were taken at each voltage. 


DESCRIPTION OF APPARATUS 


Figure 1 shows the specimen-holder, the fur- 
nace in which it was placed, and the electrometer 
circuit used for measurement. Upper electrode A, 
lower electrode B, guard ring R, and connecting 
rod E were all made of Chromel, to withstand 
the temperatures used. Rod E slides freely in the 
Chromel bushing D, which is fastened by two 
Chromel nuts K to the upper porcelain plate P. 
The entire upper assembly, including electrode 
A, is removable for insertion of the specimen Sp. 
With the specimen in place, the upper assembly 
may be set on the three guide rods U (two only 
are shown) and fastened in place with the leveling 
nuts S, whereupon rod E is moved through bush- 
ing D until contact is established with the 
specimen. The ball and socket joint J ensures 
even contact between electrode A and_ the 
specimen. To provide electrical connection, 
Chromel wires were hammered into holes drilled 
in electrode B and guard-ring R, and a third 
wire was fastened by collar /] to rod E. The 
Chromel-Alumel couple C, by means of which 
the specimen temperature was measured, was 
fastened in contact with guard-ring R. Electrode 
B and guard-ring R were accurately centered and 
permanently fastened with Chromel screws to 
the lower porcelain plate P. The separation thus 
achieved was considered sufficient insulation for 
the guard-ring, since the compensating resistor 
R. was used to maintain it at the same potential 
as the electrode B. The conductivity of the entire 
specimen-holder was measured with no specimen 
in place and the electrodes separated by about 
the thickness of a specimen, and was found 
negligible in comparison with the specimen con- 
ductivities later measured. 

The specimen-holder was mounted in an 
electric furnace with Alundum walls, roughly 
cubical in shape, with the inside edge measuring 
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eight inches. Approximately a foot of mineral 
wool was used to separate the furnace from its 
outer container. This container was lined with a 
grounded electrostatic shield ES, with only the 
specimen connections, furnace power wires, and 
thermocouple leads emerging, through thin 
Alundum tubes 7. 

The specimen leads were carried directly to the 
electrometer circuit, the source of potential, and 
the calibrating circuit (not shown in Fig. 1). A 
Westinghouse RH-507 electrometer tube” was 
used in the simple arrangement shown, with 
Westinghouse WL-756 special resistors (sput- 
tered and mounted in evacuated glass tubes) in 
the input circuit. A sufficient selection of input 
resistors to cover the range of currents measured 
(10-'* to 10-° ampere) was provided, although it 
was found that at the higher currents the gal- 
vanometer G, having a sensitivity of 10° mm per 
ampere, could be used directly with almost as 
great an accuracy. Battery Eo» and resistors Re 
and R; were used to balance out the normal plate 
current of the RH-507, so that the galvanometer 
read only increments due to currents passing 
through the specimens. The universal shunt S$ 
provided a number of changes of range with the 
same input resistor. Linear amplification was 
achieved on each range by adjustment of the 
tube voltages. 


RH-507 


VV 
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Fic. 2. Calibrating circuit. 


In using the RH-507 for the measurement of 
very low currents, it is important that the grid 
bias voltage be so adjusted that no grid current 
flows. A check on this was furnished by the 
calibrating circuit shown in Fig. 2. Resistor Reg 
was less than 10* ohms: of much lower value than 
the WL-756 input resistor, which for the smallest 
current readings was 10" ohms. Thus, if grid 


% Bulletin TD 83 (1937), Westinghouse Electric and 
Manufacturing Company. 
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Fic. 3. Detail of specimen. 


current was flowing and switch S, was closed (.S2 
remaining open), the grid potential would 
change, producing a deflection of the galva- 
nometer in the plate circuit. Accordingly, the 
potentiometer across battery E,. was adjusted 
until no deflection was noted upon the closing 
of Sj. 

For calibration, switch S: was closed, S; 
remaining closed after the grid bias adjustment. 
R, and Rx were decade resistance boxes, so used 
that the total resistance Ry+R,4+Rz, in series 
with battery E was 2X10* ohms. Thus the cali- 
brating potential applied to the RH-507 grid was 
0.5X10-*K Re XE. The voltage of battery E was 
checked, with the circuit in operation, by means 
of a thermionic voltmeter, also used to check the 
potential applied to the specimens. 

As a source of adjustable potential, a heavily 
filtered transformer-rectifier circuit was used, 
with a General Radio Company ‘‘Variac’’ con- 
trolling the input voltage. Batteries were sub- 
stituted for several check readings, but no 
appreciable difference in the readings was noted. 
The potential source, the calibrating circuit, the 
electrometer circuit, and the furnace control 
circuits were all thoroughly shielded and bonded 
to the furnace shield, such precautions proving 
quite necessary for proper operation of the elec- 
trometer tube. 


EXPERIMENTAL TECHNIQUE 


Four specimens were tested, two of KBr and 
two of KCl. All were flat plates two cm square, 
with thicknesses as follows: KBr I, 3.64+0.02 
mm; KBr II, 2.78+0.02 mm; KCI I, 4.62+0.02 
mm; and KCI II, 4.14+0.02 mm. After grinding 


VOLUME 11, MARCH, 1940 


to the tolerance indicated, the specimens were 
uniformly sputtered with platinum. As shown in 
Fig. 3, the sputtering was then removed from 
the ring R separating the contact areas of lower 
electrode and guard-ring, and from the four 
faces of width f¢. 

By taking some care to reach thermal equi- 
librium, it was found easily possible to hold the 
temperature constant to within 0.5°C (=0.02 mv 
on the precision potentiometer used for measure- 
ment) during a run. After constant temperature 
was achieved, one operator adjusted, measured, 
and applied the potentiai, and signaled time 
intervals. A second operator noted the electrom- 
eter readings at intervals of a few seconds. This 


. procedure, while somewhat crude, was found 


sufficiently accurate for the determination of the 
first-throw current, and for current-time plots 
from which the value of the constant m in Eq. 
(3) could be deduced. After taking readings over 
a period of several minutes, the potential was 
removed and the specimen short-circuited. 
Readings at the next voltage were not taken 
until the residual (reversed) charging current 
was reduced to a negligible proportion of the 
recorded currents. These precautions became 
less necessary at the higher temperatures, where 
the current was essentially constant with time. 

Readings were taken on each specimen at a 
number of temperatures between 25°C and 600°C. 
At the lower temperatures potentials between 20 
and 250 volts were applied, but at the higher 
temperatures the maximum voltage was 100, to 
avoid the possibility of breakdown. 


RESULTS 
a. Reduction in cold conductivity 


Measurements in the range 25°C-50°C were 
taken on all crystals, before heating to the 
higher temperatures, after heating to the highest 
temperature, and after several reheatings. In two 
cases the cold conductivity seemed to rise 
sharply after heating, but this was found due tc 
poor original contact between specimen and 
electrodes. With this difficulty remedied, it was 
found in all cases that the first heating (taking 
from 36 to 48 hours) reduced the cold conduc- 
tivity by 15 to 25 percent, room temperature 
conductivity being reduced in greatest degree. 
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No reduction in the conductivity at temperatures 
above 50°C was noted, and a second heating 
produced no further change in the cold con- 
ductivity. This effect seems similar to that 
reported by Smekal* for rocksalt, and by 
Rahimi'® for both rocksalt and quartz. The 
reduced conductivities are those used in the 
curves discussed below. 


b. Current-voltage relation 


No appreciable departure from Ohm's law 
was noted for any of the specimens tested. In 
Fig. 4 first-throw galvanometer deflection (= con- 
stant Xcurrent) is plotted against voltage for 
five typical cases. Only in the case of specimen 
KBr II, measured at 29°C, are deviations from 
the linear relation at all prominent, and even 
there they are small enough to be ascribed to 
experimental error. 


Current-time relation 


In Figs. 5 and 6, where the logarithm of the 
resistivity in ohm-cm is plotted against the 
reciprocal of the absolute temperature for KBr 
and KCl, respectively, small curves are super- 
imposed upon the main graphs to show the 
approximate current-time relations for the 
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Voltage Applied 


Fic. 4. Current-voltage relation for all specimens. The 
ordinates must be multiplied by a different constant in 
each case, to obtain the conduction current in amperes. 


© NM. G. Rahimi, J. de phys. et rad. 9, 291 (1938). 
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various regions. Eq. (3) describes the dependence 
of current on time rather closely. In fact, it 
appears from the data that the effect of increas- 
ing temperature is to reduce the value of the 
exponent m in Eq. (3) from about unity at the 
lowest temperatures to successively lower values, 
until (at about 200°C in all cases) m is so small 
that current decay with time was no longer 
measureable. The current-time curves can only 
be considered a first approximation, however, 
and it is planned to use a more accurate experi- 
mental method for a thorough investigation of 
the current-time relationship for these and other 
crystals at various temperatures. 


d. Resistivity-temperature relation 


Since Ohm’s law seemed to hold, a resistivity- 
temperature curve for one voltage only was 
plotted for each of the four specimens tested. 
The voltage chosen was 100, for convenience in 
calculation. A great variation in resistivity with 
temperature was found, as was to be expected: 
from nearly 10" ohm-cm at room temperatures 
it fell to less than 107 ohm-cm for KBr and less 
than 10° ohm-cm for KCI at the highest tem- 
peratures used. 


~ 
te 


8 


of Resistivity ta Cm 


Rec Degrees Kelvin 


Fic. 5. Resistivity-temperature relation for KBr, with ap- 
proximate Current-time Curves. 
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Figure 5 shows the results for KBr. Both 
specimens obey the exponential law of Eq. (2), 
and the curves have, within experimental error, 
the same slope, although the intercepts are quite 
different, constant A of Eq. (2) having the value 
20.0 for KBr I and 2.24 for KBr II. Constant B, 
on the other hand, is 12,200 for KBr I and 12,300 
for KBr II, a difference of only 0.82 percent. The 
activation energy deducible from these constants 
is 24.5+0.2 kcal. per mole. 

The curves for KCl, shown in Fig. 6, are not 
so simple. Eq. (2) might be said to hold on either 
side of a critical temperature (235°C for KCI I 
and 208°C for KCl ID), at which an abrupt 
change of slope occurs. Similar, though smaller, 
deviations from the Rasch and Hinrichsen equa- 
tion have been reported for sodium," potassium," 
and lithium!® halides by Phipps and his co- 
workers, and may be explainable by the assump- 
tion that at the lower temperatures only the 
potassium cations take part in the conduction, 
while at the higher temperatures the chlorine 
anions become active. On this assumption, 
Smekal®® has proposed as a modification of Eq. 
(2) the form 


+ A og 82/7 (4) 


where B, and By are constants whose respective 
values depend on the activation energies of the 
two kinds of conducting ions present, but this 
form seems inadequate to describe the conduc- 
tivity-temperature relation found for KCI in the 
present work. Furthermore, while at the higher 
temperatures; i.e., to the left of the critical 
point in Fig. 6, the slopes for KCI I and KCI II 
agree quite well, this is not the case at the lower 
temperatures. This may indicate that in the low 
temperature region the conductivity is more a 
function of the particular specimen tested than 
of the material of which it consists. In the higher 
temperature region the slope of the resistivity 


7 Phipps, Lansing and Cooke, J. Am. Chem. Soc. 48, 112 
(1926). 

18 T. E. Phipps and E. G. Partridge, J. Am. Chem. Soc. 
51, 1331 (1929). 

' —D. C. Ginnings and T. E. Phipps, J. Am. Chem. Soc. 
52, 1340 (1930). 

20 A. Smekal, Zeits. f. tech. Physik 8, 561 (1927); Zeits. 
f. Elektrochem. 34, 472 (1928). 
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Fic. 6. Resistivity-temperature relation for KCI, with ap- 
proximate current-time curves. 


curve is 10,050 for KCI I and 9900 for KCI IT, 
corresponding to an over-all activation energy 
of 45.8+0.2 kcal. per mole. 


e. Summary of results 


From the data taken, it may be concluded 
that: (1) Ohm’s law holds for KBr and KCl, at 
temperatures between 25°C and 600°C, and at 
fields of from 0.05 to 0.90 kv,/cm. (2) There is 
an exponential relation, most prominent at lower 
temperatures, between current and time, for KBr 


and KCl. (3) An exponential relation holds for - 


KBr between conductivity and temperature in 
the range tested. (4) There is an exponential 
relation between conductivity and temperature 
for KCl, but the value of the exponential con- 
stant changes at a critical temperature between 
200°C and 250°C. 

The author wishes to acknowledge the valuable 
assistance of Professor M. F. Manning, who 
collaborated in planning and beginning this 
work, and of Mr. Fred Willard, who assisted in 
construction of apparatus and the taking of data. 
The specimens used were obtained through the 
courtesy of Professor H. M. Randall. 
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Studies in Lubrication 


VII. The Lubrication of Plane Sliders of Finite Width 


Muskat, F. MorGan anp M. W. MERES 


Gulf Research & Development Company, Pittsburgh, Pennsylvania 


(Received August 30, 1939) 


The theory of the lubrication of plane sliders of finite width—thrust block bearings—has 


been developed by a method that is more convenient for numerical calculations than that of 
Michell and Boswall. Whereas previously all the derived properties of such systems have had 
to be obtained by direct graphical and numerical procedures, explicit analytical expressions 
have been obtained in the present theory for the friction coefficients, the inclination of the 
wedges, the load-carrying capacities, the minimum film thickness, and the lubricant flow. The 
theory has been applied numerically to a much wider range of conditions than has been done 
before. In particular, with respect to self-aligning thrust bearings or sliders, numerical results 
have been obtained so as to cover both the range where the pivot is very near the center of the 
slider and where it is set within short distances of the trailing edge. With such extended ranges 
it has been possible to construct complete curves for the variation of the important derived 


INTRODUCTION 


HI: theory of the viscous lubrication of plane 

sliders of finite width, as based upon the 
Reynolds equation, was first given by A. G. M. 
Michell in 1905.' Except for some additional 
numerical calculations carried out by means of 
Michell’s formulas by Martin,? and Boswall,’ no 
extensions or improvements have been made in 
the theory. Duffing* has developed a solution for 
the mathematically simpler case where the film 
viscosity is proportional to the film thickness, 
and has pointed out that in some cases, at least, 
the first derivative term in the Reynolds equa- 
tion does not greatly affect the final solution. 
However, no advance was made in the solution 
of. Michell’s original problem. 

The main disadvantage of Michell’s theory is 
that it not only reduces to an inherently nu- 
merical method, but moreover requires inde- 
pendent numerical or graphical integrations for 
the determination of every feature derivable from 
the solutions for the pressure distribution, such 
as the load-carrying capacity, friction coeff- 
cients, the center of pressure, and the lubricant 
flow. As a consequence, the oniy calculations 


'A.G. M. Michell, Zeits. f. Math. Physik 52, 123 (1905). 

2H. M. Martin, Eng. 109, 233 (1920). 

?R. O. Boswall, The Theory of Film Lubrication (Long- 
mans, Green and Co., 1928). 

4G. Duffing, Handbuch der Physik Tech. Mech., Vol. 5 
(1931), p. 839. 
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properties of the slider system as functions of the pivot position. 


heretofore made on sliders of finite width have 
referred to the two cases where the ratios of 
inlet to outlet film thicknesses have been either 
3.0 or 2.0. The results so obtained, for various 
length to width ratios, have been summarized 
by Boswall. 

Because the plane slider of finite width is 
essentially nothing else than the thrust block 
bearing so widely used in industry, it is desirable 
to have a more comprehensive treatment of its 
lubrication properties. In the following we shall 
therefore develop solutions of the Reynolds equa- 
tion for such sliders which are more convenient 
for numerical calculation than those of Michell. 
In principle they do not differ greatly from 
Michell’s solution. However, they give explicit 
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ries representations of the resultant pressure 
distribution in terms of Bessel functions which 
have been completely listed in standard tabula- 
‘ions. The one troublesome step in the solution 
is the determination of the roots of certain tran- 
scendental Bessel function equations which must 
be calculated separately for each ratio of the 
inlet to outlet film thickness. However, once 
these have been found, the subsequent calcula- 
tions for variable slider widths can be made 
very simply. 


THE SOLUTION OF THE REYNOLDS EQUATION 


The Reynolds equation appropriate to the 
lubricant film between plane sliders may be 
written as: 


0 dap dp oh 
( )+ ( 
Ox\ Ox Oz\ Ox 


where (cf. Fig. 1) p is the pressure distribution in 
the film, x is the coordinate in the direction of 
motion of the slider, h is the film thickness, 
u its viscosity, z the transverse coordinate, and 
l’ is the velocity of the sliding plate, which is 
assumed to be moving in the direction of the 
negative X axis. 

Assuming a constant viscosity, and a film 
thickness given by: 


h=mx, (2) 
(1) reduces to: 


ep 3dp ap 
+ = (3) 
Ox? xdx dz? 


Michell’s method consisted essentially in ex- 
panding both p and the right side of Eq. (3) as 
lourier series in z, thereby obtaining a series of 
ordinary inhomogeneous equations in x. While 
the homogeneous solutions of these equations 
were the first-order Bessel functions of the third 
kind, the particular integrals were expressed 
simply as power or Laurent series. The coeff- 
cients of the various terms of the solutions were 
then determined numerically so as to give 
vanishing pressures at the leading and trailing 
edges of the film, the Fourier series terms being 
so adjusted initially that the pressure distribu- 
tion automatically vanished at the lateral edges 
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of the film. In this way solutions were obtained 
in the form of a series of functions known only 
numerically, and which had to be computed 
separately for each different set of physical 
parameters specifying the lubrication system. 

Here we shall formally eliminate the inhomo- 
geneity of Eq. (3) by first subtracting off from 
the solution that corresponding to the slider of 
infinite width, that is by making the substi- 
tution: 


(x; —x)(x—Xp9) 


m?(x;+2Xo) 


where x; represents the inflow or leading edge of 
the film and xo the outflow or trailing edge. 
The equation for p is clearly the homogeneous 
equation corresponding to Eq. (3). We take its 
solution to be of the form: 


Ou . C,U,(x) cosh a,z (s 
) 
m? x cosh a,w/2 


where 
U,(x) J \(anXo) Yila aX) — Yila J ait), (6) 


J, and Y, being first-order Bessel functions of 
the first and second kind, respectively.’ The 
slider width is taken as w and the x axis is 
assumed to lie in the plane of symmetry. 

The constants a, and the coefficients C, must 
now be chosen so that the resultant series of 
Eq. (5) will satisfy the boundary conditions of 
the problem. Since p itself must vanish on all 
four sides of the film, j must vanish on the 
leading and trailing edges, and must cancel off 
the distribution in the infinitely wide slider at 
the lateral edges. These conditions will be 
satisfied if the a, are chosen so that: 


=0, (7) 
and if the C, are determined so that: 


(x;—x)(xX9—x) 
> C,U,(x)= (8) 
(xi: +X0)x 


5 The Bessel functions appearing here are those of the 
first and second kind, rather than of the third kind ap- 
pearing in Michell’s analysis, because of the use of hyper- 
bolic functions in the terms involving z rather than 
trigonometric. This difference is introduced here to ac- 
celerate the convergence of the series from the numerical 
point of view. 
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By carrying out the standard manipulations for the determination of the coefficients of 
orthogonal expansions, it is finally found that the pressure distribution can be expressed as: 


m? (x;-+xX0)x? 2 x cosha,w/2 
where 
J ila f Yola »x)dx — Yila f Jola nX)dx. (10) 
x0 xo 


By applying Eq. (8) to the first term in Eq. (9) the pressure distribution can also be expressed as: 


¥ U(x) | cosh | 
lz | (11) 


m* J — J 17( cosh a 2 | 


THE DERIVED PROPERTIES OF THE PLANE SLIDER 


Once the pressure distribution has been determined, the other properties of the system can be 
obtained by integration processes. First, however, it is convenient to introduce the following addi- 
tional notation :° 


=6,; (Brp), (12) 
2 
where L is the length of the slider in the direction of motion and W is the total load supported by it. 
From these definitions it will be clear that w represents the dimensionless width of the slider, d= ~ 
corresponding to the infinitely wide slider, and w= 1 to the square slider. The dimensionless parameter 
S characterizes the physical conditions to which the lubricant film is subjected, since it includes the 
relative speed of the slider, the load which the film must support, and the viscosity of the film, all 
of which are to be considered as being externally imposed on the system. 

With this notation the load supported by the film will be given by: 


ari Uw 2(p—1) 4 tanh 2 


S 1 6 


j m* 2(p—1) 4 6,” tanh B,(p—1)& 2 (14) 
log p— 


so that 


It will be noted that here as well as in Eqs. (9) integrating over the bounding surfaces the stress 
and (13) the series terms represent the effect of | components 
the finiteness of the slider width. Because of the 


factor 1 w, its contribution will vanish as the a oa P (15) 
slider becomes infinitely wide, as indeed it ; h 20x 

should. 


The frictional forces may be obtained by Where the upper sign refers to the moving plate 
and the lower to the stationary slider. These 
_ © The parameter p, which is also equal to the ratio of the components are essentially nothing more than 
inlet to outlet film thickness, is equivalent to 1 plus the 
parameter 6 introduced by Boswall. From the manner in the surface tractions in the lubricant film at the 
which itenters the subsequent analysis it will be clear that — hounding surfaces, with the signs adjusted so as 
it is more convenient for our purposes than would be the : : : 
parameter b. to take account of the direction of motion of 
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he moving slider and to give directly the 
\umerical magnitude of the individual external 
‘orces required to keep the moving slider moving 
and the fixed plate stationary. Carrying out 
‘hese integrations formally, evaluating the 
second terms by integrations by parts, and 
introducing the above notation, it is found that 
the friction coefficients for the moving slider will 
be given by 


— =~—+— log p, (16) 
m 2 m? 


and that on the stationary slider by: 


f 
= ——-+-— log p. (17) 
m 2 m? 


It should be noted that the latter refers to the 
frictional force as measured parallel to the sta- 
tionary slider. Defined in this way its numerical 
value will become negative for those values of 
S m?® where f/m is less than 1, since the differ- 
ence between f/m and f/m will always equal 1. 
Such negative values of the friction coefficient 


9 


Ax 1 6571 
xo p—1 m?*L2 


4 


Wp(p—1)* 


imply physically that the stationary slider tends 
to be pushed back by the viscous film in the 
direction opposite to the moving slider. This 
apparently paradoxical behavior may be readily 
understood by a closer analysis of the velocity 
distribution in the film under such conditions 
where the negative values of f are obtained. 
At the same time, however, it should be pointed 
out that because the normal stresses on the 
slider, which provide the load-supporting power 
equal to W, are also inclined to the vertical by 
the angle between the sliders, they provide a 
horizontal component which counterbalances the 
negative contribution due to f. Moreover they 
leave a resultant friction coefficient parallel to 
the moving slider which is not only positive, but is, 
in addition, exactly equal to that on the moving 
slider. 

Another property of interest of the lubricant 
film between the plane sliders is the center of 
pressure. By symmetry it obviously must lie 
midway between the lateral ends. Its position 
along the central traverse in the direction of 
motion is given by: 


log p 


Ji(Bn) — pJi(Bnp)} tanh 
—Ji?(Bnp)} 


where again the series term’ represents the effect of the finiteness of the width of the slider and 
vanishes as the slider becomes infinitely wide (w— ~ ). 

Finally it is of interest to derive the expressions for the lubricant flow between the sliders. From 
the basic hydrodynamic expressions for the average flow between parallel plates, and from which 
I-q. (1) is derived, it is readily found that the total lubricant inflow at the leading edge will be 


given by: 


h® ap Uh UmwLp 
o- | ( — —-+- ) dz= [ 
—w/2 Ox 2 x=x —1 


2(p+1) 


tanh B,(p—1)wv/2 
Brit —J17(Brp)} 


w(p—1) 


While the outflow at the trailing edge will be given by evaluating the integral of Eq. (19) with 
the integrand taken at x=Xo, it is more convenient to consider directly the side leakage, which may 


be calculated from the equations: 


h® ap 2UmwL 
xo \12p O27 wW(p—1)? 


— pJi(B,)} tanh B,(p—1)v/2 
—J:2(Bnp)} 


(20) 


7 The terms outside the series, here, and in Eqs. (14), (19) and (20), which represent the values for the infinitely wide 
slider, could be included in the series merely by replacing tanh 8,(p—1)@/2 by: tanh 8,(p—1)@/2—8,(p—1)@/2. 
However, this change would slow up the convergence of the series and would mask the effect of the side leakage as a 


correction factor for the corresponding systems of infinite width. 


VOLUME 11, MARCH, 1940 


211 


p 
| 
| 
+_ 
a 
= 


TABLE I. Values of 8,. 


= 


e. 1 2 3 4 5 6 7 S 
1.2 15.7277 31.4259 47.1305 | 62.8368 | 78.5438 94.2511 
15 | 6.3218 12.5861 | 18.8628 | 25.1427 | 31.4239 | 37.7057 
2.0 3.1917 6.3116 | 94446 | 12.5812 15.7199 18.8595 
3.0 1.6356 3.1789 | 4.7381 | 6.3028 | 7.8696 9.4375 
5.0 0.84715 | 1.61105 2.38532 | 3.1642 3.9454 | 4.7279 | §.$112 
10.0 0.394094 | 0.733057 1.074838 1.418864 1.764331 | 2.110731 2.457757 2.80522 


As it should, the whole of this term vanishes as 
the slider width is made _ indefinitely large. 
Moreover, if Q. be subtracted from the total 
inflow, as given by Eq. (19), the flow out of the 
trailing edge will be obtained. 

These complete the formal derivations of the 
expressions for the significant physical properties 
of plane slider lubrication systems. We shall now 
proceed with the problem of evaluating these 
expressions numerically and with the interpreta- 
tion of the final results. 


PLANE SLIDERS WITH FIXED WEDGE ANGLE 


As already indicated, the main problem in- 
volved in the evaluation of the analytical ex- 
pressions derived in the last section lies in the 
determination of the roots 8, of Eq. (7). In 
terms of the notation of Eq. (12) this equation 
may be rewritten as: 


Ji(Bn) Vi(Bnp) — (21) 


As these roots are not given in the standard 
texts on Bessel functions, and as this type of 
equation appears in other physical problems, 
those that were used in the present investigation 
are reproduced in Table I. The values for p= 1.2, 
1.5, and 2.0 were taken from a paper by A. 
Kalahne,* after correcting Kalahne’s first three 
roots for p= 2.0. The roots for the other values of 
p were calculated explicitly for the purpose of 
the present problem. It should be emphasized 
that these roots, especially the lower ones and 
those for large values of p, must be known with 
a high degree of accuracy. In fact, for p=10, 
even the seventh decimal places in the lower 8,, 
as determined by graphical interpolation, had 
some noticeable effect upon the left side of 
Eq. (21). Moreover all the subsequent calcu- 
lations, and in particular those of the functions 


8 A. Kalahne, Zeits. f. Math, Physik 44, 55 (1907). 
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5,, must be carried out to a high degree of 
accuracy or else the final results will be quite 
meaningless. While this is a rather tedious 
problem, it can fortunately be carried out with 
satisfactory precision by making use of the 
Bessel function tables given in Watson's treatise® 
and in the British Association Tables.'® 

The functions y, defined by Eq. (10), and 
appearing in 6,, were obtained by expressing 
the integrals of Eq. (16) in terms of Struve 
functions, i.e., as: 


2x0 
(Bnp) 


11\(B,p) 
2 


| 
| (22) 
2 


the J/,; being the Struve functions of the first 
order. Account was taken of the fact that the 
8, are-roots of Eq. (21). These Struve functions, 
up to arguments equal to 16, were taken from 
the tables of Watson. For higher values asymp- 
totic expansions were used. For such values 
where both of the Struve functions in Eq. (22) 
permitted the use of asymptotic expansions the 
resultant values of 6, took the form: 


1 


3 
B 2 


n 


Ji(Bn) 
-| 


45 Ji(Bn) 1 
+— | (23) 
p® 


n 


Once the 8, and 6, are calculated, it is a rela- 


*G. N. Watson, Theory of Bessel Functions (Cambridge 
University Press, 1922). 

10 British Association Mathematical Tables, Vol. 6, 
Bessel Functions, Part 1, 1937. 
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Fic. 2. The variation of the friction coefficient, f, with S, 
for sliders of finite width. S=uUw/W; w=film viscosity; 
('=slider speed; w=slider width; L =slider length, in the 
direction of motion; W=total load; m=slope of slider 
inclination; @=w/L. 


tively simple matter to carry through the evalua- 
tion of the final physical quantities of interest 
such as the load carried by the film, the friction 
coefficients, center of pressure, etc. Moreover the 
major elements in the terms of the various 
series remain the same for the calculation of 
different slider widths, as these can be taken 
care of merely by changing the arguments in the 
hyperbolic tangent terms. 

For the slider systems of fixed wedge angle the 
final parameter used for plotting the various 
physical properties of the system was S/m? as 
ziven by Eq. (14), although this in itself is 
determined by the dimensionless slider width W, 
and the ratio of the maximum to minimum film 
thicknesses, p. 

In Fig. 2 is plotted a series of curves showing 
the variation of the friction coefficient vs. S/m? 
as calculated by Eqs. (14) and (16). The curve 
for W=% corresponds to the infinitely -wide 
slider as given by the Reynolds theory. It will 
be seen that for fixed values of S/m? the friction 
coefficient (divided by m) increases with de- 
creasing slider width. The magnitude of this 
increase becomes more pronounced for larger 
values of S/m*. This effect is caused, of course, 
by the side leakage, which in turn tends to lower 
the average film pressure and thus necessitates 
that the surfaces come closer together in order 
\o support the external load. 

As in the case of the infinitely wide slider, the 
friction coefficient increases with S/m* more 
slowly than linearly. In fact, except at very 
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small values of S/m?, the various curves are 
essentially parallel on double logarithmic paper, 
so that the magnitudes of the friction coefficient 
may be approximately expressed as a function of 
the width times a function of S/m? corresponding 
to the infinitely wide slider. Thus taking the 
values for the infinitely wide slider as unity, the 
values of f/m at S/m?=50 are 1.14, 1.33, 1.77, 
2.16, and 2.56 for #=2, 1, 3, 3, and }, respec- 
tively. That is, a square slider will have a friction 
coefficient roughly 4/3 as great as that of the 
infinitely wide slider, and one whose width is 
only } as great as its length—but with the same 
value of m—will have coefficients about 23 
times as great as the slider of infinite width. 
These ratios break down at very low values of 
S/m*, and in the limit of vanishing S/m?, f/m 
approaches the value 3 regardless of the slider 
width. 

The effect of the slider width on the minimum 
film thickness is shown in Fig. 3. Here, as would 
be expected, the minimum film thickness uni- 
formly decreases with decreasing slider width, 
this effect being the consequence of the side 
leakage. Moreover, as in the case of the friction 
coefficients, the curves are approximately parallel 
on double logarithmic paper, so that the values 
of ho/mL may also be roughly expressed as a 
function of the width times a function of S/m? 
corresponding to the infinitely wide slider. At 
Sm? = 50 these width factors are 0.86, 0.70, 0.48, 
0.37, and 0.29 for #=2, 1, 3, 3, and 4, respec- 
tively. For vanishing values of S/m* the minimum 
film thickness will in all cases go to zero. 

The variation of the center of pressure with 


P ri |_| 
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Fic. 3. The variation of the minimum film thickness, 
ho, with S, for sliders of finite width. Notation same as in 
Fig. 2. 
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S/m®* is plotted in Fig. 4. The ordinates are the 
distances of the center of pressure from the 
trailing edge divided by the slider length in the 
direction of motion. It will be seen that the 
center of pressure is shifted closer and closer 
toward the trailing edge as the slider width 
decreases. Here the curves are roughly parallel on 
the Cartesian scale, so that the various curves 
are related to each other by differences which are 
approximately constant, rather than uniform 
factors. For very small values of S m?*, however, 
the steepness of the curves completely vitiates 
the parallelism obtaining for large values of S/m?. 

For fixed large values of S the total 
lubricant flow into the slider decreases with 
decreasing slider width (cfs Fig. 5). For small! 
values of S m®* there is a crossing of the curves 
and a reversal of the variation with the slider 
width. That is, for large values of S m?*, Q 
decreases monotonically with decreasing slider 
width. For S m?=10, Q increases first as the 
slider width decreases, reaches a maximum, and 
then decreases as the width decreases still 
further. Moreover the maximum becomes shifted 
to decreasing widths as S m®* is made to ap- 
proach 0. 

As would be expected, the side leakage in- 
creases with decreasing slider width for all values 
of S/m*. Moreover within the errors of the 
calculations the magnitude of the side leakage is 
essentially independent of S m?*, for large values 
of the latter. However, as S' m* approaches 0 the 
curves develop maxima and ultimately fall off 
again before the Q, axis is reached. 


» 40 50 i= 8) 70 60 30 100 
S/mé2 


Fic. 4. The variation of the center of pressure with S, 
for sliders of finite width. Ax=distance of center of 


pressure from trailing edge. Rest of notation same as in 
Fig. 2. 
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Fic. 5. The variation of the lubricant flow with S for 
sliders of finite width. Q=Q/UwL; Q.=Q./UwL; Q=total 
lubricant inflow at leading edge; Q.=total side leakage. 
Rest of notation same as in Fig. 2. 


SELF-ALIGNING SLIDERS OF FINITE 
TRANSVERSE WIDTH 
It will be clear that the results just presented 


may be given a direct physical interpretation 
only when the slope of the inclined surfaces, as 


represented by m, is considered to be fixed or 


known. As under the conditions of practical 
operation the value of m will be extremely small, 
it is not convenient to achieve fixed and pre- 
determined wedge angles in practice. A much 
more practical procedure is to permit the sliders 
to be self-adjusting and to take up whatever 
mutual inclination is necessary to support the 
externally applied load. As originally pointed out 
by Michell and Kingsbury, this can be attained 
by pivoting the load supporting slider and per- 
mitting it to rotate about the pivot line. Since the 
external load must act through the pivot line, the 
supporting film pressures will also have to have 
their center of pressure along this line, and, 
depending upon the magnitudes of the speed, 
load, and lubricant viscosity, will force an 
inclination of the pivoted sliders so that the 
center of pressure will fall under the pivot line." 

This is the principle of the pivoted slider block. 
By its definition the system is essentially de- 
termined by the position of the pivot, which 
must coincide with the center of pressure, the 
inclination of the surfaces being merely a second- 

" To obtain exact equilibrium of the moments about the 
pivot line the actual center of pressure would have to be 
displaced somewhat from the pivot line because of the 
moments arising from the friction forces on the slider. In 


practice, however, where the height of the pivot above the 
sliding’surface is not large, this effect will be quite negligible. 
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Fic. 6. The variation of the slider 


inclination, m, with the pivot posi- . m o6 LX N : \ 
‘ion for pivoted sliders of finite » 

Zz 
width. Ax/L=(distance of pivot 0s » % \ N 
from trailing edge)/(slider length, \ 


L.). Rest of notation same asin Fig. 2. 


fundamental independent variable is the quantity 


Ax/L, which gives the position of the center of 
pressure or the distance of the pivot line from the 


12(p—1) 24 
D*(p) =6 log p—- 


ary parameter. For such systems therefore the 


p+1 w(p—1) 
we will have directly from Eqs. (14), (16), (17), (19), and (20); 


trailing edge, although from an analytical point 
of view the fundamental parameter is p, which 
determines Ax/L through Eqs. (14) and (18). In 
fact, if we define the function D(p) by: 


6,2 tanh B,(p—1)w/2 Ax 
=p*( ). (24) 


(25) 


p D log p 
s=(- + S; j=(- (26) 
2 D 2 D 
Q Dp 2(p+1) 6,J1(B,) tanh B,(p—1)w,2 
Q. 2D — tanh B,(p—1)0/2 
ULw wW(p—1)* nt J17(Bn) —J1?(Bnp)} 


These relations show that all the fundamental 
characteristics of the pivoted sliders including 
the slope m, the friction coefficients f and f, the 
minimum film thickness fo, and both the total 
lubricant flow Q and side leakage Q., are pro- 
portional to \/S for a fixed-pivot position. The 
effect of the pivot position and the slider width 
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enter only in the coefficients of Eqs. (25)—(28). 
These are plotted in Figs. 6-11. 

It will be seen from Fig. 6 that for all values of 
the slider width the sliders will adjust themselves 
to greater wedge angles as the pivot is moved 
closer to the trailing edge.’* Thus for a square 


The ending of the various curves in Figs. 6-11 at 
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Fic. 7. The variation of the friction 


coefficient, f, with the pivot position for 
# | pivoted sliders of finite width. Ax/L 
5 at at 4 = (distance of pivot from trailing edge) 
ba | /(slider length, L). Rest of notation 
| same as in Fig. 2. 
iy 
| } 
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slider a shift of the pivot position from 0.45 to 
0.40—-measured as fractions of the total length 
from the trailing edge—will increase the wedge 
slope, for the same values of S, by a factor of 3. 

The wedge angle will also increase for a fixed- 
pivot position with increasing slider width. In 
particular, for the fixed-pivot position of 0.40, 
cutting the width of the square slider by two will 
decrease the slope of the wedge by a factor of 2.1 
for the same load or speed, that is, S. This 
behavior is clearly due to the fact that the higher 
film pressures obtaining in the sliders of greater 
width permit the load to be carried with the 
greater wedge angles. As indicated directly by 
Kq. (13), if the pivot position (p) and width are 
kept fixed the load-carrying power will vary 
inversely as the square of the slope of the wedge 
angle. 

The variation of f/\/S is shown in Fig. 7, as a 
function of the pivot position. The general 
behavior of all the curves indicates a sharp fall in 
f/\/S as the pivot position is moved away from 
the center of the slider, very flat minima in the 
curves for pivot positions less than 0.39 of the 
slider length, and rather moderate rises to higher 
values of f \/S as the pivot is moved still further 


different values of Ax/L is due to the fact that the values of 
Ax/L for the individual values of p—and in particular for 
the largest value, p= 10—decrease with increasing #. As 
previously indicated, the calculations were made by 
assuming p and computing the Ax/Z corresponding to it by 
Eq. (18). 
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toward the trailing edge. In the case of the slider 
of infinite width, the minimum in the curve lies at 
Ax/L=0.39. For the sliders of finite width the 
minima are shifted to smaller values of Ax/L, 
and for #<1, the curves are still falling over the 
range of Ax /L covered by the calculations. The 
absolute values of f/\/.S uniformly increase with 
decreasing slider width for large values of Ax/L. 
Thus for the square slider the friction coefficient 
with fixed S for a pivot position of 0.40 will be 41 
percent greater than for the infinitely wide 
slider, and for a slider width equal to } of the 
length the coefficient will be 3.48 times as great 
as for the infinitely wide slider. However, for 
small values of Ax/L the curves tend to merge 
and ultimately cross. In fact, the infinitely wide 
slider curve is crossed by those for #=2 and 1 at 
Ax L>0.26, and even that for #=} should cross 
after a small extrapolation. 

Figure 8 gives the corresponding variations of 
the minimum film thickness which is, of course, 
that at the trailing edge. Here, regardless of the 
slider width, the minimum film thickness first 
rises as the pivot position is moved away from the 
center of the slider, reaches a maximum, and then 


falls rather slowly as the pivot position ap- 
proaches the trailing edge. For large values of 
Ax L the minimum film thickness decreases with 
decreasing slider width, but as Ax/L is decreased 
the curves for finite width ultimately cross that 
for w= ~. Moreover, the maxima of the curves 
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Fic. 8. The variation of the mini- 


mum film thickness, Ao, with the 


pivot position for pivoted sliders of T 
finite width. Ax/L=(distance of UWS 55, | | 

pivot from trailing edge) /(slider 

length, L). Rest of notation same as 016 4 t = 
in Fig. 2. } 


shift to smaller values of Ax/L with decreasing 
slider width, and the rise to and fall from the 
maxima become increasingly more gradual as the 
slider width decreases. Hence both from the 
point of view of obtaining minimum friction and 
maximum values of the trailing edge film thick- 
ness the pivot positions in the finite width sliders 
should be shifted closer to the trailing edge than 
is indicated by the theory of the infinitely wide 
slider. Even so, the absolute values of the 
minimum film thickness will be materially less in 
the sliders of finite width near their maximum 
values. Thus the maximum value of the trailing 
edge thickness in the square slider will be only 67 


percent of that for the infinitely wide slider with 
the same value of S, while sliders twice as long as 
they are wide will have maximum trailing edge 
thicknesses only 44 percent of those of the 
infinitely wide slider. 

At first thought, the result that a slider of 
finite width could operate with higher trailing 
edge thickness and lower friction coefficient than 
the infinitely wide slider with the same pivot 
position and value of S is indeed surprising. For 
the pressures in a film open at the sides must 
certainly be lower than in one of the same 
geometry except for its infinite width. And 
indeed this is also directly supported by Fig. 3, 


Fic. 9. Relative load-carrying ca- 


pacities of pivoted sliders of finite 

width. W/W...= (load per unit width 

0.7 /(load per unit width carrie 

= slider of infinite width); Ax/L 

=(distance of pivot from trailing 

x edge) /(slider length, L); w=slider 

05 Ae, width; #=w/L; solid curves: load- 

ts carrying capacities for sliders of 

04 N equal minimum film thickness; 

| = dotted curves: load-carrying capac- 

03 = NY PINS ~ ities for sliders of equal frictional 
force per unit width. 
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Fic. 10. The variation of the total 
lubricant flow with the pivot posi- 
tion for pivoted sliders of finite 
width. Q=Q/UwL; Q=total lubri- 
cant inflow at leading edge; Ax/L 
=(distance of pivot from trailing 


edge)/(slider length, L). Rest of 
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notation same as in Fig. 2. 


which shows that for fixed values of m and hy the 
load-carrying powers of the finite width films will 
be less than those of the infinite width film. 
However, it must be noted that in Figs. 7 and 8 
‘the films compared are for sliders with the same 
pivot position, or center of pressure. The actual 
geometry of the film will then be materially 
different for the sliders of different width. In 
particular, as Fig. 6 shows, the slopes m will be 
much smaller for the finite width sliders than for 
the infinitely wide slider, and this drop in m will 
be especially large for small values of Ax L. 
Since the smaller values of m inherently lead to 
larger load-carrying capacities, this factor will 
ultimately counterbalance the direct end effect of 
the open sides of the finite width slider and enable 
the latter to carry ‘the applied load with a 
minimum film thickness that is even higher than 
that: required for the infinitely wide slider. And 
associated with the higher values of the minimum 
film thicknesses, is the tendency of the friction 
coefficients to fall below that of the infinitely 
wide slider. However, the actual crossing of the 
friction coefficient curves does not develop as 
soon as in the case of the minimum film thickness 
because the friction forces depend on the average 
film thickness as well as on the minimum 
thickness. 
A somewhat more practical way of expressing 
the results shown in Figs. 7 and 8 may be had by 
.considering the relative load-carrying capacities 
of the sliders of finite width as compared to those 
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of infinite width. Thus supposing that the sliders 
are of the same length, have the same pivot 
position, and operate at the same speed with 
lubricants of the same viscosity, the ratios of the 
loads per unit width carried by the finite and 
infinite width sliders for equal trailing edge film 
thicknesses will be given by the reciprocals of the 
squares of the ordinates of Fig. 8. These ratios are 
plotted as the solid curves in Fig. 9. While these 
curves are subject to some error because of the 
graphical interpolations required for their calcu- 
lation, there can be little question with regard to 
their main features. Thus, except for the shallow 
minima in the region of large Ax /L, the relative 
load-carrying capacities will be seen to increase 
monotonically as the pivot line approaches the 
trailing edge. And for the smaller values of w this 
increase is so marked that the value unity is 
exceeded even within the range of the calcu- 
lations. This corresponds, of course, to the 
crossing of the curves of Fig. 8, previously 
discussed. In particular, if Ax L =0.40 the square 
slider will be able to carry a load, per unit 
width, only 46 percent as great as the infinitely 
wide slider, for the same speed, film viscosity, and 
minimum film thickness. And if the slider width 
is cut to } its length the load-carrying capacity 
will be reduced to only 5 percent of the infinitely 
wide slider. However, if Ax/L=0.29 or less, the 
load-carrying capacities of all sliders whose 
widths are not less than } of their length will 
even exceed that of the infinitely wide slider. 
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Fic. 11. The variation of the 


total side leakage with the pivot 


position_for pivoted sliders of finite 


width. Q.=Q,./UwL; Q,=total side 


leakage; Ax/L=(distance of pivot 


from trailing edge)/(slider length, 


L). Rest of notation same as in 


Fig. 2. 


Considering sliders having the same total 
frictional force per unit width, the relative load- 
carrying capacities will be as indicated by the 
dotted curves of Fig. 9. These were obtained from 
the ratios of the reciprocals of the squares of the 
ordinates of Fig. 7 rather than Fig. 8. It will be 
seen that for practical purposes these relative 
load-carrying capacities are pretty much the 
same as those for which the minimum film 
thickness is kept fixed, although the rise for the 
sliders of small widths at small values of Ax/L 
is not at all as sharp as in the case of the solid 
curves. 

In Fig. 10 are shown the ratios of Q/\/S for 
different slider widths as functions of Ax/L. The 
quantities Q are the dimensionless total lubricant 
flows, as defined by Eq. (27). Here it will be seen 
that the details of the curves depend markedly 
upon the geometry of the slider. Thus the value 
of Q for the infinitely wide slider increases to a 
maximum at Ax/L=0.39 as the pivot position is 
moved away from the center of the slider and 
then decreases as the pivot approaches the 
trailing edge. For sliders of relatively great width 
the general behavior is rather similar, although 
the maxima in the curves are shifted toward 
considerably lower values of Ax/L. For narrow 
sliders, however, the curves show inflections and 
continue to rise with decreasing values of Ax/L 
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even up to values as low as 0.20. As functions of 
the slider width, for fixed values of Ax/L, the 
quantities Q/./S will therefore show monotonic 
increases with increasing width when Ax/L is 
high, but will reach maxima at finite values of the 
width and fall back again to lower values as the 
width is made indefinitely large, when the pivot 
position is appreciably shifted from the center of 
the slider. 

The corresponding results for the dimensionless 
side leakage, Q., are shown in Fig. 11. Here there 
appears to be a monotonic increase in side 
leakage as the pivot position is moved over 
toward the trailing edge, the side leakage always 
being zero when the pivot is at the center of the 
slider. In its variation with the slider width, 
however, it again increases first as the slider 
width is diminished, reaches a maximum, and . 
ultimately falls again as the sliders are made still 
narrower. With respect to the absolute magnitude 
of the side leakage, it may be noted that while it 
is necessarily smaller than the total lubricant 
inflow Q, it becomes of the same order of magni- 
tude as Q except when the pivot position is very 
close to the centers of the sliders. 

The authors are indebted to Dr. Paul D. 
Foote, Executive Vice President, Gulf Research 
& Development Company, for permission to 
publish this paper. 
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The Formation and Maintenance of Electron and Ion Beams 


Litoyp P. SmitH AND PauL L. HARTMAN 


Department of Physics, Cornell University, Ithaca, New York 


The radial potential distribution in an ionic beam of 
circular cross section is calculated and the maximum beam 
current which can be obtained in a beam of given radius 
and boundary conditions is computed. In the case of ions 
the maximum beam currents may be quite small. The ions 
in such beams have a considerable velocity distribution 
which in turn leads to a greater beam divergence than 
former calculations indicate. In Section 3 a method is 
worked out for holding the beam together during the initial 


INTRODUCTION 


a of electrons and positive ions have 
assumed a place of great importance in 
modern electronic devices and in researches in 
physics and biology..In these fields it has become 
more and niore important to increase the current 
density in the beam. Considerable difficulty has 
been experienced in obtaining positive ion beams 
of high intensity and in some cases even with 
electron beams. In fact, the following analysis 
will show that there is a definite upper limit to 
the current that can be obtained in a beam of 
given cross section.under specified voltage condi- 
tions. For positive ions this current is quite small. 

The problem of beam formation divides itself 
naturally into five parts: (1) that of keeping the 
ions together during the period of initial accelera- 
tion; (2) the space charge divergence throughout 
the remainder of the beam length; (3) that of 
preventing divergence over the high velocity 


‘ 


| 


Fic. 1. Key for the designation of the various 
beam potentials, 


(Received September 1, 1939) 


accelerations. A disadvantage of cylindrical lenses is 
pointed out in this connection. The use of electrostatic 
lenses for holding together a beam composed of high energy 
particles is discussed in some detail. Calculations show the 
special arrangements of the lenses necessary for producing 
convergence of the beam. A few other arrangements to 
prevent divergence which were tried experimentally are 
described. 


region; (4) the radial potential or velocity dis- 
tribution in the beam and (5) the maximum 
beam current that can be realized. These aspects 
of the problem will be considered in detail in 
the following sections and various practical 
methods of holding the beam together will be 
discussed. 


1. RADIAL POTENTIAL DISTRIBUTION AND 
MAXIMUM CURRENT IN BEAM* 


In a beam of electrons whose center lies on 
the axis of a cylindrical surface of definite poten- 
tial, the potential distribution will in general 
depend on the distance r from the axis and the 
distance z along the beam. This general case will 
not be treated here but it will be assumed that 
the potential is a function of the radius only. In 
practice this would correspond to a case where 
the spreading of the beam was small or to a 
situation in which the electrons or ions were made 
to travel parallel to the axis of the beam by the 
application of an intense magnetic field parallel 
to the beam axis. Under these circumstances 
Haeff ' has determined the potential distribution 
through a two-dimensionally infinite lamina of 
finite thickness. 

Consider a beam of electrons of radius p 
coaxial with a cylindrical conductor of radius R. 
The potential at the center of the beam will be 
taken as zero and the potential of the cylindrical 


* The calculations in the first two sections were com- 
pleted by the first author while temporarily associated with 
the Research and Engineering Department, RCA Manu- 
facturing Company, Harrison, New Jersey. 


' A. V. Haeff, Proc. I. R. E. 27, 586 (1939). 
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surface as Ver while the potential difference 
between cylinder and source of electrons will be 
denoted by Vo as shown in Fig. 1. The velocity 
of an electron or ion of mass M at a distance r 
from the axis is 


VetV(r) 


and the charge density is = —i/v where 7 is the 
constant current density. In cylindrical coor- 
dinates, Poisson’s equation for the resulting 
potential distribution is 


IdsdV M\! 
rdrX\ dr m 


x- 
([2e m Vo—Vet+Vin) 


The above equation may conveniently be ex- 
pressed in terms of the dimensionless variables 
s and ¢ where 


r=as and g=V(r)/(Vo— Vr) (1) 
and 
(m/ M)*(2e m)} 


Vo— (2) 


In these variables the differential equation 


becomes 
ld/fde 1 
(: )- - (3) 
sdsX\ ds (1+ ¢)} 


Before discussing the actual solution of this 
equation, it is advantageous to discuss the con- 
ditions at the boundary of the beam and their 
relation to the total beam current. In the region 
p<r=R the charge density is zero and it is 
readily found that the potential in this region is 
given by 

(Vr—V,) log p+V, log R/p 
log R/p 
. 


where VV, is the potential at the beam boundary. 
This potential must smoothly join the potential 
for r<p at r=p which requires that 


dV’ /dr},.,=dV/dr} 
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Fic. 2. Curve a shows gas a function of u. Curve 6 shows 
¢ as a function of s. Heavy straight line is a plot of Eq. (5). 
Dashed line indicates procedure for determining Sp. 


In terms of the variables ¢ and s, this relation 


becomes 
(=) Vre-V,1 
a ds ,, log Rip p 


R 
¢o= ¢r—log (5) 
p ds/ , 


or 


where 
¢,=V,/(Vo—Vr) and s,=p/a. (6) 


Eq. (5) may be regarded as an equation for the 
determination of the value of s corresponding to 
the beam boundary for proper matching. This 
together with Eq. (2) then determines the total 
beam current in terms of Vo and Vr. The form 
of Eq. (5) suggests the following procedure for 
determining s,. Assuming that the solution of 
Eq. (3) is known and letting sdg/ds=u, ¢(s) is 
plotted as a function of u as shown by curve a 
in Fig. 2. Eq. (5) may be written 


¢=¢r—log R/p-u 


and is represented by the straight line in Fig. 2 
whose slope is —log Rp. The intersection of this 
line and curve a determines the values of ¢ and 
u at the beam boundary, namely ¢, and u,. By 
plotting ¢ as a function of s in the same figure, 
curve 6b, the value of s corresponding to ¢, can 
be found as indicated by the dotted line in Fig. 2. 

With the value of s, thus determined, the 
value of the total beam current J as a function 
of Vr/Vo can be found, since a=p/s and the 
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Fic. 3. Ve/ Vo for maximum beam current as a function 


of R/p. 


relation (2) give 


6.07 X10' om)*T(amp.) ( 
} s,”. (7) 


'V(volts)! 

The solution of Eq. (3) cannot be given in 
terms of elementary functions. For values of s 
not too large a useful series solution satisfying 
the conditions g=dg ds=0 at s=0 may be 
found in the usual way. The first three terms in 
this series give 


s* 7s! 
g= ~+ (8) 
4 2:4" 4°-6° 


This gives values of ¢ accurate to less than 1 
percent even at s=2. For large values of s the 
equation was integrated numerically to provide 
the curves plotted in Fig. 2. 

In order to calculate s, analytically, (when s, 
is not too large) Eq. (8) is substituted in Eq. 
(5), which after some rearrangement gives 


1 4y 1+4 log R/ p 
1+2 log R/ 2-4" 


] R 
p 


The value of s, can be obtained by solving this 
equation by the method of successive approxima- 
tions starting with 


1 
142 log R/p 1—x 


as the first approximation, where x= Vpr/ Vo. 
Three approximations are warranted here and 
the third approximation yields 


x | x {x . 
$,°=A ,1+B +¢( (9) 
1-—x 


where 


1 1 
A=--: 
414+2loge R p 
1+4 log R p 
B= (10) 


(142 log R p)? 
7 14+6log R p 


16L2(14+2 log R/p)* 9 (142 log R/p)! 


From Eq. (7) the total beam current is propor- 
tional to (1—x)!s,2=w. When Cx (1—x)<1 the 
value of x which makes w a maximum can be 
determined in the usual way and its value is 
readily found to be 


2 1 14+4 log R o 
(x= Vp/ Ve) mex [i+ | (11) 
3 4(14+2 loge R p)* 


The slowly varying function of R/ p on the right 
is plotted in Fig. 3. Eq. (11) furnishes the rather 
surprising result that the maximum beam current 
does not occur when Ve= Vo but Vr~jVo. By 
comparing the analytical values of x for maximum 
I with those found by the graphical method, it 
has been established that the values computed 
from Eq. (11) are very accurate for the entire 
range 1=R/p< a. 

In Fig. 4, w has been plotted as a function of 
Ve/ Vo for several values of R/p including the 
interesting case when R/p=1. The beam current 
for any given value of ‘9 can then be calculated 
by means of relation (7). It is interesting to note 
that for R/p=10 and Vo= 1000 volts, the maxi- 
mum beam current of molecular hydrogen ions 
is only 2.6 ma. 

The radial voltage distribution in the beam is 
of very great practical importance. On account 
of it, the electrons or ions have a considerable 
range of velocities. It is more useful to express 
the potential at any point in the beam with 
respect to the source of electrons or ions. This 
potential will be denoted by V, and may be 
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expressed in terms of the potentials already used 
by writing 


Vo- Vat V(r) = (12) 


or Vs/Vo=U—Ver/Vo)(A+¢). ¢ may be ex- 
pressed in terms of r/p by replacing s in (8) by 
s=r/a=(r/p)s,, which gives 


Ve 1 
Vo 4 8 4 


7 
+——(s,?4)?(r, r/p=1. (13) 
4-62 


The function of r contained in the brackets | } 
is normally a rather slowly varying function of 
r especially for R/p>>1, so that the voltage dis- 
tribution in the beam is rather closely represented 
by 


Ver 5," 
V, ). r p=1. (14) 
Vo 4 


The voltage distribution outside of the beam 


a | 


| | | 
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Fic. 4. w, which is proportional to beam current J, 
plotted as a function of Ve/Vo. Arrows indicate value cf 
Vr/Vo for maximum w as computed by Eq. (11). Crosses 
indicate values of w computed from Eqs. (7) and (9). 
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Fic. 5. Potential distribution from center of beam to 
boundary cylinder for R/p=4 and maximum beam 
current. 


1=r/p=R/p is readily found from (4) so that 
V, Vo=(1— Vr/ Vo) 


(Vr/ Vo— V,/ Vo) log r/p+(V,/ Vo) log R/p 
——,, (15) 
log R/p 


1=7r/p=R/p. 


The voltage distribution for maximum beam 
current for R/p=4 as computed from expressions 
(13) and (15) is shown in Fig. 5. This clearly 
brings out the fact that all of the electrons in 
the beam are traveling with velocities very 
much different than that corresponding to the 
applied voltage Vo. This is an extremely im- 
portant matter in certain uses of electron and 
ion beams. In particular when the confining 
magnetic fieid is removed the beam diverges 
faster than one would expect. This is shown in. 
the next section. 


2. DIVERGENCE OF THE BEAM 


The divergence of an electron or ion beam 
produced by the mutual repulsion of the charged 
particles in the beam has been computed by 
several people.2 In these computations it has 
been assumed that all of the electrons or ions in 


the beam are traveling with the same velocity. 
The analysis of the preceding section shows that 
this will not be the case except for extremely low 
beam currents. The velocity distribution of the 


2 E. E. Watson, Phil. Mag. 3, 849 (1924). M. Knoll and 
E. Ruska, Ann. d. Physik 12, 604 (1932). F. D. Fowler and 
G. F. Gibson, Phys. Rev. 46, 1075 (1934). 
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electrons in the beam will cause the beam to 
diverge more rapidly than the above-mentioned 
computations would indicate. In calculating the 
divergence, it will be assumed that the electrons 
retain the velocity distribution characteristic of 
the initially confined beam. This is justified 
provided the divergence is not large. 

The radial force acting on an electron initially 
at the outside edge of the beam may be found by 
differentiating expression (15) with respect to r. 
This gives 

Vre-V,1 
E,= = 
log pr r 


(Vo-—Ver)u, 
(16) 


The term on the right was obtained with the help 
of Eq. (5). The equation of motion in the radial 
direction is 


(lo— 
_ (17) 
dt? r 


The equation of the path may be obtained by 
writing v,=dz dt where z is the coordinate along 
the beam axis and 2, is the velocity of the electron 
initially at the beam edge. Then 


d*r (Vo—Ver)u, 


Integrating this equation twice gives 


s=2p{ —-—- (18) 
0 


The expression for z obtained by Fowler and 


Gibson? may be written 


(log r p)} 
=2p(2, wf 


0 


(19) 


where w is defined by (7). The ratio of the 
distances required for the beam to diverge the 
same amount in the two cases is 


(1+ %,)w)! 


2u, 


¢ 
~ 


(20) 


This ratio is plotted as a function of w in Fig. 6 
for the case R/ p= 10. The ratio is less than unity 
for all values of w up to its maximum value 
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except at w=0 where the ratio is unity. Thus the 
two computations are in agreement only for 
extremely small beam currents. At maximum 
beam current (w=20.3 in this case) the beam will 
diverge a given amount in about one-half the 
distance computed by Fowler and Gibson. 


3. PREVENTION OF BEAM DIVERGENCE 
DURING ACCELERATION 


In the early stages of beam formation the 
electrons or ions start with a relatively low 
velocity and are accelerated. The natural diver- 
gence computed in the last section is especially 
marked at low velocities and must be prevented. 
It will be shown in this section that it is always 
possible in principle to hold the beam together 
during acceleration. It is only necessary to 
produce a radial force toward the center of the 
beam which is somewhat greater than that 
produced by the charge in the beam. Such a 
radial field arises whenever the gradient of the 
electric field in the beam direction increases with 
distance along the beam. 
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Fic. 6. Ratio of distances for beam to diverge a given 
amount taking into account and not taking into account 
the voltage distribution through the beam. R/p=10. 
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When Gauss’ theorem is applied to a cylin- 
drical surface of length Az coaxial with the beam 
of not too large a radius, it is readily found that 
the radial component of electric field is given by 


— 0E,/02)r, 


where ¢=i/v is the charge density and v is the 
velocity of the particles in the beam and £, is 
the component of electrical field along the axis 
of the beam. Expressing E, and v in terms of the 
potential at a given point we arrive at the fol- 
lowing differential equation for the way in which 
|’ must vary in order that E,=0: 


—=0,. 
dz* (2Ve/m)} 


The signs are to be so chosen that Ve>0O and 
id V/dz<0. The solution of this equation subject 
to the boundary conditions V=dV/dz=0 at 


z=0 is 
34/8 2/3 
(2e 


This may conveniently be written 
V(volts) = 5680! i(amp.) | (21) 


where W/m is the ratio of the mass of the par- 
ticle considered to that of an electron. It is of 
interest to note that this is the potential vari- 
ation given by the Child-Langmuir equation for 
the case of a diode with infinite parallel plane 
electrodes. 

A convenient practical way of obtaining this 
variation of potential along the axis is to pass 
the beam through holes in a succession of closely 
spaced disk electrodes the potential -of which 
increases so as to produce the proper variation 
along the axis. If the potential increases faster 
than the } power of z then the beam will be 
confined to a smaller cross section. It should be 
noticed that for beams of not too large a radius 
the space charge force is compensated throughout 
the entire cross section of the beam. 

It is of considerable practical importance to 
realize that the use of cylindrical electrodes for 
rapidly increasing the velocity of relatively slow 
electrons or ions may result in the loss of a large 
number of particles from the beam by over- 
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focusing. In order to obtain some idea of the 
order of magnitude of this effect the path of an 
electron or ion was computed for a particular 
case neglecting space charge but using the other- 
wise exact equation® 


2U d’*r dV dr AV 
- — —+—-=0, (22) 
1+(dr/dz)?dz* dz dz Or 


where Ue>0 is the energy of the particle and V 
is the potential within two adjacent coaxial 
cylinders of radius R separated by a plane. 
Since this equation is invariant under a trans- 
formation of the unit of length, the unit will be 
taken as the radius R of the cylinder. For con- 
venience the initial energy Ue in this case will 
be taken to be Ue= Vo—V(ro,20) where 2Vo is 
the potential difference between the two cylinders 
and ro=0.9 and z9= —1, z being measured from 
the plane of separation of the electrodes. 

The potential within the cylinders in the 
neighborhood of the plane of separation for long 
cylinders and for electron acceleration is given by 

1 Jo(iur) sin ws 
V/Vo= 
Jo(iu) 


‘ 
n=1 Jilyn)¥n 


10V 
= —— sin 
J (tp) 


Vo or Tv 


(23) 
1 Jo(iur) 
=— —— cos 
Vo TY Jo(iu) 
=+2 


where J(r) denotes the Bessel’s function of the 
first kind and y, denotes the zeros of Jo(y). The 
convergence of the above series is sufficiently 
rapid to use them in the numerical integration 
of Eq. (22). This was done for the initial con- 
ditions dr/dz=0 at ro=9.9 and and 


3E. g., see I. G. Maloff and D. W. Epstein, Electron 
Optics in Television (McGraw-Hill, 1938), p. 80. 
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Fic. 7. Trajectory showing overfocusing when a relatively 
slow charged particle is appreciably accelerated through 
evlindrical electrodes. 


initial energy 


Vi ro,20) 
U =1— 
Vo Vo 


V(0.9, —1) 
20.02. 


This means that for 100,000 volts across the 
cylinders the initial energy corresponds to 
U=1000 volts. The resulting path is shown in 
Fig. 7. It is readily seen that an electron or ion 
attains a considerable radial component of 
velocity which it tends to retain. As a result the 
ion strikes the cylinder wall about two and one- 
half radii from the gap. This shows the over- 
focusing characteristic of cylindrical lenses for 
‘slow particles rather far from the center of the 
beam and indicates that for large initial accelera- 
tions such cylindrical 
avoided. 


electrodes should’ be 


4. Tue Use or LENS SYSTEMS 
FOR SUPPRESSING DIVERGENCE IN HIGH 
ENERGY BEAMS 


In this section the effect of a succession of 
electrostatic lenses with cylindrical symmetry 
will be investigated. It is found that the beam 
will diverge to somie extent except for certain 
special cases. The velocity of the particles will 
be assumed large enough so that space charge 
may be neglected. The focusing effect of a suc- 
cession of two-dimensional lenses has _ been 
treated by Roset who also took into account the 


*M. E. Rose, Phys. Rev. 53, 392 (1938). 
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effect of a time variation of the potential dif- 
ference across the electrodes. In the following 
analysis of electric lenses with cylindrical sym- 
metry no time variation of the fields will be 
considered. 

The differential equation of the path of a 
positive ion through the mth lens for (dr/dz)?*<«1 
is 


d dr aV 
2(U,— V')- |-- -, (24) 
dz dz 


where V is the potential at the point 7, 2, such 
that V(r,0)=0 and U, is the voltage the ion 
has fallen through by the time it has reached 
the median plane z=0 as illustrated in Fig. 8. 
In the following calculation only those ion paths 
will be considered for which the potential at any 
point in the lens may be represented by V(r,z 
= V(0,2) — 4). When this expression 
is used Eq. (24) becomes 


d*r dVidz dr 1 
—— (25) 
dz? 2(1U,-—V) dz 4(U,-—J) ds? 
It may be readily be shown? that 
r=(U,-—V)-'y, (26) 
where y satisfies the differential equation 
dy 3,dV,dz)\? 
161U,—-—V 


Since this discussion pertains to particles having 
a high energy, powers of V>/U, higher than the 
first will be neglected where 2 Vy is the potential 
difference across the lens. In this approximation 
k:q. (27) reduces to d*y/dz*=0 so that 


(28) 
and 
dr 
—=B}(U,—Virz)}- 
dz 
dV 
(29) 
dz 


This shows that the path of the particle suf- 
ficiently far from the median plane to either side 
is a straight line as shown in Fig. 8 because V 


5E.g., see E. T. Whittaker and G. N. Watson, Modern 
Analysts (Cambridge University Press, 1920), p. 194. 
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becomes constant and dV/dz=0 here. If these 
straight portions of the path are continued to the 
median plane their points of intersection with 
this plane will be a distance Ar, apart. The slopes 
of the straight line portion of the paths to the 
left and right of z=0 will be designated as 
and respectively, and the difference 
in slope will be Aa,=a@n4;—a@n_;. The effect of 
any one lens will be determined by the values 
of Ar, and Aa,. These quantities may now be 
calculated from Eq. (28). For large negative 
values of z, V(r,z)= Vo and for large positive z 
\"(r,s)= —Vo. The constants A and B in Eq. 
(28) may be determined by noting that r=r, 
for s=0 and dr/dz},--.=an_, and dr/dz} 
=«,,, and have the values 


A=U,\,, (30) 


B=ay_\| Un— (31) 


Consequently the paths of the particle on the 
right and left of the median plane are 


r1—Vir.z) 
ne=(1 -- On+12, 
U, L 14+V0/U, 


U, 
Jon 


The value of Ar,, computed from these expres- 
sions as mentioned above, is 


Vo 
Ar, = 
U, 


It is easily found that 


dr 
| 
1 Vo 


(1+V./U,)\- 
\1—Vo Val 


(33) 


n 


The change in 7, per lens depends on the slope 
of the particles’ path and distance /, between 
the (n—1)st and mth lens and the displacement 
Ar,. Thus we may write 


dr ,/dn = tArn. 
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Fic. 8. The effect of a cylindrical electric lens on the path 
of a charged particle near the axis. 


By differentiating this equation and eliminating 
a», one obtains the differential equation for the 


particles’ path through a succession of lenses 
namely 


\idl, d 
-— Ar, 


(34) 
il, dn dn 


When it is assumed that the particle receives the 
same increase in energy at each lens, U,=2Von 
and this together with expressions (32), (33) and 
(34) gives 


1 dl, 


-—1,=0, 
2n/ dn Anl, dn 


idl, 1dr, 


l, dn 


(35) 
dn? 

where terms of the order 1/n? have been neg- 
lected. For generality, let the lenses be spaced 
according to the law 


(36) 


in which case the solution of Eq. (35) neglecting 
terms in 1, * and higher is 


B, (37) 


This result is interesting in that when the law 
of spacing of the lenses is such that m> —}, the 
ion beam will always diverge but when m< —} 
the beam will converge. In the case of a linear 
resonance accelerator® the lengths are given by 


1, 


so that here the beam diverges according to the 
law 


V,=An+B. 


Another lens of interest may be handled in 
much the same manner, namely, the lens ob- 
tained by placing a plane grid over the face of 
the high velocity electrode of the previously 
considered cylindrical lens. In this case the 


6 P. L. Hartman and Lloyd P. Smith, Rev. Sci. Inst. 10, 
223 (1939). 
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n+ 

| 

| 

dn? \dn | £ 

| 

| 

| 

| 

U,, 

1 Vo | : 

rn. (32) 

2U, 

= 
q 


potential of the grid at z=0 and at greater values 
of z is assumed to be zero, and at values of large 
negative z is assumed to be 2V». Instead, how- 
ever, of taking (dr/dz).., =a,,4, to evaluate con- 
stants, it is necessary to take (dV /dz),0=an, 
and evaluate constants with this and_ the 
previous boundary conditions, the reason being 
that the derivatives of the potential are discon- 
tinuous at this point. Carrying through the 
same procedure used previously in evaluating 
the lens effects one finds the same Ar, : 


Ar,=- Vola (38) 


For Aa one obtains the somewhat different result 


lo dV 
Aa=— ( -) (39) 
2U. 4U,\ dz7 .-0 


Here (d1'/dz).~» is interpreted as being the limit 
of the electric field as z approaches zero from the 
negative side. Since (dV, dz).» is negative, this 
term gives the lens stronger focusing action than 
the ordinary cylindrical lens. 

For a succession of such lenses such as previ- 
ously considered for the symmetrical cylindrical 
lens one obtains the differential equation for 7, 


d*r,, dl, 1 
dn* l, du dn 
1 dl, l, sdV 
+—(—) =0. (40) 
4nl, dn S8nVo\ dz 


Again taking /], =/)n”" one finds the equation to 
have the solution 


‘ 


7, = 
where 
d*S 


S=0. 
dn* 81, dz 1) 


Solving this by the W.K.B. method and sub- 
stituting above gives for a succession of lenses 
of this type 


ly dV ‘ 
81, dz z=0 


lo dV 
xexp| +i/-—-(—) —| (41) 
| 8Vo\dz7 m+1 
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It is seen that this is an oscillating function the 
amplitude and wave-length increasing or decreas- 
ing depending on m. When m<0O the solution 
converges and for m>0O the solution diverges, 
the results being markedly different in this 
respect from the ordinary cylindrical lens. For 
the linear resonance accelerator referred to 
one has 


above; where m= 4 
ly fdv ‘ 
8V> dz z=0 
fdV 
| 8Vo\dz } 3 


The beam diverges very much less in this case 
than for the symmetrical cylindrical lens dis- 
cussed above. 

Instead of merely taking /, to be some par- 
ticular function of m one could also take U, to 
be some other function of than 2nVo» and 
determine requirements for a focused beam. 

All the results here obtained are independent 
of the e m value of the accelerated particle. 


5. EXPERIMENTAL ATTEMPTS AT Low 
VELOCITY FOCUSING 


In connection with the confining of low energy 
ion beams of considerable intensity, it may be 
of interest to consider some experiments made 
in this connection in work on an acceleration 
device of the linear resonance type. In the par- 
ticular device it was desired to pass the beam of 
ions from the source through a delivery tube 
about 2.5 centimeters in diameter to the accel- 
erator, a distance of about 30 centimeters. The 
source was a low voltage type’ and it was desir- 
able to refrain from the use of high voltages to 
confine the ion beam. But with low voltages of a 
couple of kilovolts or less, a beam of only about 
a microampere was obtainable at the accelerator 
when perhaps a milliampere or so was delivered 
»y the source. 

The first attempt at improvement consisted of 
mounting inside, and coaxial with the delivery 
tube, a helical grid about 2 centimeters in 
diameter and biased negatively with respect to 
the delivery tube around it. It was thought that 


7G. W. Scott, Phys. Rev. 56, 954 (1939). 
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Fic. 9. The electric field in a cylindrical grid type of 
focusing device. 


there might be sufficient field penetration from 
the outside cylinder to keep the beam confined 
within the grid. As shown in the sketch of Fig. 9 
the grid wires do act as lenses with non-image- 
forming properties. In region 1 of each wire in 
the grid there is convergent action and in region 
2 divergent action acting for a shorter period of 
time due to the particle’s higher velocity. Again, 
region 3 is diverging but also acting for a shorter 
time than the convergent action of region 4 due 
to the particle’s slower velocity in region 4. One 
cannot say easily what the net lens action will be, 
due to the stronger field in near the grid wire 
than away from it. At any rate not much im- 
provement over the ordinary tube was noted. 

A second attempt was made with another 
non-image-forming device. A 5-mil tungsten wire 
was stretched along the axis of the tube and 
biased negatively with respect to the tube itself, 
giving rise to a purely radial field driving par- 
ticles toward the axis. With about 700 volts 
potential difference between the wire and the 
tube, a beam of about 50 microamperes of 2000 
volt hydrogen ions was obtained at the accel- 
erator end. With the wire heated to copious 
electron emission temperature a current of the 
same order of magnitude was obtained’ with no 
radial field other than that due to the voltage 
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drop in the wire itself, which was less than 50 
volts. The greater part of the focusing here was 
due primarily to some ion space charge neutral- 
ization. 

A third improvement was made with the use 
of cylindrical lenses. It is, of course, well known 
that a cylindrical lens produces convergence 
whether or not it accelerates or decelerates the 
particles passing through it. With this in mind, 
a series of about twenty cylindrical lenses were 
mounted within the delivery tube, which alter- 
nately accelerated and decelerated the particles 
giving rise to a considerable net radial deflection. 
With about 300 volts between successive 
cylinders, a current of about 150 microamperes 
was obtained at the acceleration end of the 
delivery tube. 

The most satisfactory arrangement so far was 
had by surrounding the delivery tube with a 
solenoid producing a magnetic field of a couple 
hundred gauss. Electrons were fired into the 
delivery tube from the accelerator end and ions 
in at the other end. The electrons produce a 
“potential trough” through which the ions may 
pass. If both beams of particles have the same 
energy, for space charge neutralization the 
currents must be in the ratio of the square roots 
of their masses. This arrangement has allowed 
currents as high as 1.5 milliamperes of ions to be 
discharged into the accelerator system. There is 
sometimes considerable ion production in the 
tube itself, as may be evidenced by the fact that 
it performs quite well when only electrons are 
fired in at both ends. 

Work is under way at the present time toward 
the maintenance of a beam by increasing the 
potential along the axis as the 4 power of the 
distance as indicated by Eq. (21), for balancing 
out the space charge radial field. 


The Mechanism of Static Friction 


Amontons’ law of friction,! which places tangential 
frictional force directly proportional to load, has long been 
an accepted relationship, but there has been much differ- 
ence of opinion concerning its theoretical explanation. 
Coulomb? proposed that the force of friction is due to the 
opposition which the normal load produces when inter- 
locking asperities of opposing solid surfaces must be lifted 
over each other. As a second possible mechanism he 
suggested that opposition to slip might be caused by the 
interaction of the force fields of opposing solid surfaces. 
However, he rejected this latter possibility for cases where 
Amontons’ law held. 

During the past three years an extensive investigation 
of friction and related phenomena has been carried on 
under the joint cooperation of the Cincinnati Milling 
Machine Company and the University of Cincinnati. It is 
the purpose of the present publication to show, on the basis 
of the results of this investigation, that interaction of 
surface force fields is as much a part of the mechanism of 
static friction as is interlocking of surface irregularities. 

First, to harmonize the two views, let us consider 
contact between two tdealiy smooth solids of different 
curvature. There can be no interlocking of irregularities 
here. However, let it be assumed that the surface fields do 
interact at the area of actual contact, A, and thus a finite 
tangential force, F, will be required to initiate relative 
sliding. Then as a consequence, F must be proportional to 
A,or F=SA, S thus being the force per unit area necessary 
to cause slip or shear of the contact area. It follows, then, 
that if the frictional force, F, is to be also proportional to 
normal force, N, i.e., if Amontons’ law is to hold, A must 
be proportional to V, or N=HA. Thus H, the compressive 
force per unit area of actual contact, is a constant of the 
material, and so independent of load. That this last 
conclusion is generally true has been mentioned by 
Tomlinson (as having been observed by Stott)* and has 
been more fully confirmed by Bowden and Tabor.‘ It is 
readily observed ,on a macroscopic scale in ordinary 
engineering hardness testing (although the constancy of [/ 
is there somewhat obscured by the phenomenon of work 
hardening). Hence H is what is generally known as the 
mean pressure hardness of the softer of two contacting 
solids. Now, by definition, ~=F/N. From the above, 
therefore, u4=SA/HA =S/H, where ua is the coefficient of 
static friction for cases where no interlocking of asperities 
takes place. 

Let us consider contact similar to the above 
between a pair of rough flat surfaces for which the value of 
S is sero. It is apparent that the angle of repose for these 
ideal solids will be approximately equal to the average 
angle, 6, which the planes of slip (in this case the tiny 
surface ‘‘hillsides’’) make with the general plane of the 
surfaces. Thus the coefficient of friction, ug, will be equal 
to tan @. Bikerman and Rideal® have recently suggested 


now 
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this same conclusion, but have assumed that this term 
alone is the source of all friction. 

In the foregoing we have shown, then, that when 6=0, 
u=ua=S/H, and when S=0, p=ug=tan 6. It follows 
therefore that, in general, h«=yat+us=S/H-+tan 6, which 
we now propose as the general equation of static friction. 

That any interlocking of surface irregularities must have 
some effect on static friction is self evident, but the exist- 
ence of the term S//H/ in this expression for ~« demands 
thorough verification. This has been furnished in the 
course of the previously mentioned investigation of friction. 
Coefficients of static friction were determined for metal 
surfaces dry turned in a lathe and quickly transferred toa 
vacuum chamber. Strict precautions were taken to prevent 
contamination of the freshly exposed metal. The surfaces 
were therefore clean except for a certain amount of 
adsorbed gas, and were probably quite fully work-hardened 
by the machining process. Details of the experimental 
work and apparatus will be given in a more complete 
discussion to be published later. 

The results of the investigation are in agreement with 
the theory outlined above. They may be summarized as 
follows: 

(1) The coefficient of friction was independent of load 
(except for very minute loads) thus verifying Amontons’ 
law for the case of these clean surfaces, and for the range of 
loads used in obtaining the succeeding facts. 

(2) For like-metal pairs, the points of contact were 
severely torn and scored after a single slip, indicating 
actual shearing of the metal itself, and thus a very high 
value for S. 

(3) Carefully polished like-metal pairs were prepared by 
a dry rouge technique which left the surfaces sufficiently 
clean to be freely wet by water. These gave virtually the 
same coefficients of friction as the rough machined surfaces, 
even for the very first slip on an unscratched finish. This 
indicates that the effect of surface roughness is ordinarily 
negligible when the metal itself shears, (i.e., the planes of 
shear are practically parallel to the surface). 

(4) For like-metal pairs, values of » were independent 
of temperature and approximately equal to one for metals 
of cubic crystal structure tested below their annealing 
range. When the metals were heated above their annealing 
temperature, however, the coefficient of friction increased 
to much higher values, and this increase was maintained 
after cooling and even after an exposure to air equivalent 
to that experienced by a newly prepared specimen. 

Quantitative confirmation of the correctness of the 
term S/H will be published later. 


M. E. MERCHANT 
Research Department, 
Cincinnati Milling Machine Company, 
Cincinnati, Ohio, 
January 18, 1940. 


! Amontons, Mem. Acad. Roy. Sci. 206 (1699); 96 (1704). 

2 Coulomb, Mem. Acad, Roy. Sci. 161 (1785). 

3G. A. Tomlinson, Phil. Mag. 7, 905 (1929), see p. 920. 

4F. P. Bowden and D. Tabor, Proc. Roy. Soc. A169, 391 (1939). 
5 J. J. Bikerman and E. K. Rideal, Phil. Mag. 27, 687 (1939). 
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Wratten Light Filters 


for Photography. Photomicrography. 
Spectroscopy. Photometry. and Astronomy 


Wrhrasrren LIGHT FILTERS are universally used in widely sep- 
arated fields of scientific research. 

The Wratten list includes filters for use with all types of 
color-sensitive films and plates. Over one hundred different 
filters are available. 

Complete spectrophotometric data are given in the Eastman 
publication, Wratten Light Filters; price, 50 cents. A price list 
of filters will be furnished free on request. 


EASTMAN KODAK COMPANY 
Research Laboratories ROCHESTER, N. Y. 


GEIGER-MUELLER TUBES wits. WIDE PLATEAU 


by Professor Francis R. Shonka, De Paul University, Chicago 


Typical Shonka-Welch G-M Tube Curve 
Showing extremely broad plateau SPECIAL FEATURES 
ai Extremely Wide Plateau. 

Eliminates the need for 
> highly stabilized 
a. | voltage supply. 

Py z High Efficiency and 
Stability. 
I Offered in two stock sizes: 
No. 625B—short tube, 60 mm 
8 8 each, $25.00 
No. 625C—long tube, 180 mm 
No. 625-B No. 625-C “© each, $30.00 


Special sizes available on order 
For full information, write 


W. M. Welch Scientific Company 


Established 1880 


Catalog lists over 10,000 items for 
Science Laboratories 


1515 Sedgwick Street Chicago. Illinois.U.S.A. 
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Associates of the Institute 


| eee the most part, the expenses incurred in carrying out the 

work of the American Institute of Physics are defrayed by 
its Founder Societies. The institute exists to do the work which 
these Societies have assigned to it—work toward objectives which 
they have in common. These Societies are being materially aided 
by a significant added contribution from the dues of the Asso- 
ciates of the Institute. The Associates are a group of forward 
looking companies and laboratories who believe it is valuable to 
them and to America to maintain a vigorous advance of physical 
science. In aiding the Institute they are helping physicists to 
secure for themselves and their science the necessary oppor- 
tunities and facilities for effective service to industry and to 
society in general. 


The Associates of the Institute are 


Agfa Ansco Corporation 
Allegheny Steel Company 
Aluminum Company of America 
American Telephone and Telegraph Company 
Bausch & Lomb Optical Company 
Bell Telephone Laboratories 
Crane Company 
Dow Chemical Company 
E. I. du Pont de Nemours & Company 
Eastman Kodak Company 
General Electric Company 
General Motors Corporation 
Gulf Research & Development Company 
Libbey-Owens-Ford Glass Company 
Loomis Laboratory 


Massachusetts Institute of Technology, 
Department of Physics 


Massachusetts Institute of Technology, 
Division of Industrial Cooperation 


Radio Corporation of America 
Research Corporation 
Shell Development Company 
Sperry Gyroscope Company, Inc. 
Westinghouse Electric & Manufacturing Company 


The American Institute of Physics cordially invites other 
interested corporations to become its Associates and will welcome 
their inquiries addressed to,the Director. 


American Institute of Physics 
175 Fifth Avenue, New York, N. Y. 
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MATTER, MOTION AND ELECTRICITY 


A Modern Approach to General Physics 
By H. D. Smyru, Princeton University, and C. W. Urrorp, Allegheny College 


647 pages, 6x9. $3.75 


this outstanding new book, intended primarily for students with a background of 

high school physics, the authors review fundamental principles and definitions, 
while presenting a fresh approach that avoids a repetition of the conventional beginners’ 
course. Modern physics is introduced earlier than usual. A feature of the book is 
the wealth of illustrative examples and stimulating problems. 


“This book seems to be the first attempt in many years at a searching reorganization 
of the materials of physics suitable for a general introductory course of college level, 
and it is one of the first to present a consistent and clean cut organization of the subject 
in keeping with the now well established and commonly accepted modern concepts 
regarding matter and energy.” 


—Professor J. C. Hussarp, Johns Hopkins University 


Send for a copy on approval 


McGRAW-HILL BOOK COMPANY, INC. 


330 West 42nd Street New York, N. Y. 


RAMAN SPECTROGRAPH 


Excellent resolution and permanence of adjustment characterize 
the GAERTNER Raman Spectrograph. Supplied with either 12” 
focus F/3.5 or 20” focus, F/6.3 camera. 


Write for catalogue L2. 


THE GAERTNER SCIENTIFIC CORPORATION 
1212 WRIGHTWCOD AVE. CHICAGO, U.S.A. 
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FIRST SURFACE 


MIRRORS 
by High Vacuum Deposition 


EVAPCO first surface are excellent 
for projec tion mirrors and other uses 
that require an opaque mirror of 
high reflectivity that does not tar- 
nish or give double image. 


DUOLUX semi-transparent for one 
exposure three color cameras, range 
finders, or wherever an accurate 
beam-splitter is needed. Made to 
order for all practical transmission 
~—-reflection ratios. 


EVAPORATED 
METAL FILMS CORP. 
ITHACA, N. Y. 


B -POTENTIOMETER | 


A high-grade, general-purpose laboratory po- 
tentiometer with 3 ranges: 1.6, .16 and .016 volts. 
Extremely accurate and convenient. Used in lead- 
ing standardizing laboratories. 

Described in Bulletin No. 270 which also 
lists Rubicon portable potentiometers. 
Other products: standards of resistance, in- 
ductance and capacitance, Wheatstone and Eel- 
vin bridges, resistance boxes, galvanometers, 

electrometers, permeameters, colorimeters. 


RUBICON COMPANY 


29 North 6th Street Philadelphia, Penna. 


TheNew Signal 
Generator 

| Range 120 KC 
to 120 MC 

| In Fundamentals 


New service instrument is important companion 
to the Rider Chanalyst and Rider VoltOhmyst 


Dollar for dollar, feature for feature, the new RCA 
SIGNALYST.. . is the best buy in Signal Generators. Its 
amazing range is greater than any test oscillator. 
Its accuracy and stability are excellent... Stray signal 
leakage is kept at a minimum... AC operated with 
regulated power supply! It is an instrument that will 
find many applications in the laboratory in $] 0750 
Net 


radio and television work. 


Over 335 million RCA radio tubes have been pur- 
chased by radio users. In tubes, as in parts and test 
equipment, it pays to go RCA All the Way 


RCA Manufacturing Co., Inc., Camden, N. J. 
A Service of Radio Corporation of America 


STABILIZED A-C VOLTAGE 


INPUT 95 130 volts 
Instantaneous action 
A magnetic unit 


115 volts +'4% 
Stabilizes at any load 
within rating 
? No moving parts. Nothing to wear out. 
RAYTHEON Voltage Stabilizer 


APPLICATIONS Broadly it insures stable operation 


of all precision apparatus obtaining its power from 
an A.C. source, for example :— 


Insures constant brilliancy from all types of 
lamps. 


Improves the operation of X-ray equipment. 
Stabilizes all electronic apparatus. 
Write for Bulletin 48-71 JP. 


RAYTHEON MFG. CO. 
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Send Us Your “Specs” to get the 


Don’t worry about that Rhe- Ri ‘Ri 

ostat problem! Ohmite En- ight heostat 
gineers, with their wide 1 
cialized experience, can hel 
you solve it quickly, 
ably, economically! Send in 
your“Specs”—often the right 
Rheostat can be found in the 
wide range of Ohmite stock 
types and sizes. Or, Rheo- 
stats can be made-to-order 
promptly to meet your exact 
requirements. 


OHMITE MFG. CO. 


4888 FlournoySt.,Chicago 


be Kight with OHMITE 


* RESISTORS + TAP SWITCHES 


Wrue for Catalog 17 


An eight pound rock salt prism blank cut from 
a twenty-five pound single crystal 


SYNTHETIC CRYSTAL OPTICS 


For infrared and ultraviolet spectroscopy and for 
polarizing optics... . Single crystals of Lithium 
Fluoride, Sodium Chloride, Sodium Nitrate, and 
Potassium Bromide. . . 


. Write for information. 


THE HARSHAW CHEMICAL CO. 


SPECIAL PRODUCTS DIVISION CLEVELAND, OHIO 
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EPPLEY 
THERMOPILES 


for 


RADIANT ENERGY 


MEASUREMENTS 


3. 


| 


A view of some of the apparatus used for 
calibrating thermopiles at the Eppley Lab- 
oratory, Inc. 


Adequate manufacturing, testing and 
calibrating equipment in the hands of 
our skilled and experienced technicians 
insures satisfaction to the users of Ep- 
pley Thermopiles. 


We are prepared to furnish thermopiles 
of either bismuth-silver or copper-con- 
stantan, mounted in vacuum or air type 
cases. We also furnish thermopiles 
and mountings especially designed to 
meet individual needs. A card will 
bring you our Bulletin R-3 which con- 
tains illustrations, sample curves, com- 
plete descriptive matter, prices, refer- 
ences, etc., with regard to Eppley 
Thermopiles. 


THE EPPLEY LABORATORY, INC. 


SCIENTIFIC INSTRUMENTS 


NEWPORT, 
U.S.A. 
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Write for 
Bulletin 1620-R 
which lists new 

lower prices. 


JAGABI 
RHEOSTATS 


have a graphited 
carbon lubrica- 
tor associated 
with the contact 
brushes. Elimi- 
nates scratching 
and abrasion— 
gives smoother 
operation and a 
better contact. 
Makes Jagabi 
Rheostats better 
rheostats. 


A 
LUBRICATED 

SLIDING 
CONTACT 


JAGABI RHEOSTATS 
are made in 3 sizes 
and in 76 different 
ratings—all regularly 
carried in stock. 


JAMES G. BIDDLE CO. 
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Eyes for the Future look through 


POLAROID* 
Light Controls 


Already in use in variable windows, anti-glare day 
glasses and desk lamps, Polaroid developments fore- 
shadow 3-dimensional color movies, glareless auto 
headlights, many other important new devices. In- 
struction in light polarization has logically become a 
part of every science course, elementary or advanced. 


For school and college science depart ments, Polaroid 
filters come in economical, convenient form, with 
kits for performing many dramatic experiments. 


POLAROID CORPORATION 
720 Main Street, Cambridge, Mass. 


(@) U. S. PAT. OFF. 


INDEX TO ADVERTISERS 


Name Page 


Makers of Microscopes, Microtomes, Colorimeters, Re- 
fractometers, Spectrometers, Balopticons, Photomicro- 
graphic and Microprojection Apparatus and related in- 
struments. Also makers of Orthogon Eyeglass Lenses 
for Better Vision. 


* Jagabi "’ Rheostats; Adam Hilger and Kipp & Zonen 
Optical Instruments; Pointolite Lamps; Electrical 
Testing and Speed-measuring Instruments. 


CENTRAL SCIENTIFIC COMPANY ........ccccccee Cover 4 


Manufacturers of Cenco Physical Apparatus and Instru- 
ments to meet all requirements of University, College 
and High School Physics Laboratories. Specializing in 
high vacuum pumps and development of instruments 
and apparatus for various sciences. 


THE 


Purified Organic Chemicals for research purposes; Plates 

for Photography, Photomicrography, Spectroscopy, Pho- 
tometry, Astronomy; Wratten Light Filters; Cameras 

and Films. 

Standard Cells, thermopiles, pyrheliometers and tempera- 

ture bridges. 


EVAPORATED METAL FILMS CORPORATION .......... v1 


First surface and semi-transparent mirrors by high vac- 
uum deposition. 


GAERTNER SCIENTIFIC CORPORATION 7 Vv 


Spectroscopes, Spectrometers, Spectrographs, Spectropho- 
tometers, Heliostats, Measuring Microscopes, Compara- 
tors, Cathetometers, Reading Telescopes, Interferometers, 
Chronographs, Dividing Machines, etc. 


Consens, COMPARE Cover 3 


Manufacturers of electronic measuring instruments; 
vacuum-tube voltmeters, amplifiers and oscillators; wave 
analyzers, noise meters and analyzers, stroboscopes; lab- 
oratory standards of capacitance, inductance and _fre- 
quency; irapedance bridges, decade resistors and con- 
densers; air condensers and variable inductors; rheostats, 
Variacs, transformers; other laboratory accessories. 
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HaARSHAW CHEMICAL COMPANY ..... 


Polaroid Kit No. 326, $10, 


Polaroid Kit No. 324, $25. 


Write Division 4 for free 
handbook and survey of 
applications. 


Optical Lithium Fluoride grown in single crystal to pre- 


determined size. Industrial Chemicals. 


Leeps & NorkTHRUP COMPANY ......... 


Manufacturers of Galvanometers, Resistors, Bridges, 


Condensers, Inductances, Potentiometers, 


Testing Sets; 


Temperature Measuring, Recording and Controlling Ap- 
paratus; Instruments for Measuring and Controlling 
Conductivity of Electrolytes and Hydrogen Ion Con- 


centrations. 


McGraw-Hitt Book Company, Ine. 


OuHMITE MANUFACTURING COMPANY .. 


Manufacturers of close control rheostats, 
justable power resistors, precision and 


fixed and ad- 


non-inductive 


resistors, attenuators, tapswitches and R.F. chokes. 


PoLAROID CORPORATION 


Polaroid experimental kit for demonstration. Polarized 


light. Handbook. Sheet Polaroid. 


RAYTHEON Mrc. COMPANY ........... 
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Voltage stabilizers and regulators. Electrical measure- 


ment and control instruments. 


RCA MANUFACTURING CoMPANY, INc. 
RCA Oscillators and Oscillographs, RCA 


ment, RCA Ultra-Sensitive DC Meter. 
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Test Equip- 
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Galvanometers, electrometers, potentiometers, Wheatstone 
and Kelvin bridges, resistance boxes, hydrogen ion and 


conductivity apparatus. 


W. M. Wetcu Scientiric COMPANY .. 
Scientific instruments—laboratory apparatus. 


10,000 items. 
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IMPROVED NULL DETECTOR 
for a-c Impedance Bridges 


SENSITIVE . RUGGED e CONVENIENT TO OPERATE 


HE G-R TYPE 707-A NULL DETECTOR uses a one-inch cathode-ray tube in a new non-inductive 


degenerative amplifier circuit, with tuning and phasing networks and sweep and sensitivity con- 
trols. This new detector offers many advantages which include: 


Operation in noisy locations Shows immediately any drift of either or both 
Not affected by strong fields components 
. ~ a Provides positive indication of the direction of 
off-balance for either component as selected 
Separately indicates balance of the resistive Can be calibrated to show the degree of un- 
and the reactive components balance 
Makes possible precise balancing of ¢ither Can be used at all times at maximum sensi- 
component with only moderately close bal- tivity, even with the bridge far off balance 
ance of the other Supplies instantaneous response 
Precise measurement of the steady component Will withstand any overload caused by marked 
can be made while the other varies errati- off-balance, and is instantaneous in re- 
cally covery 


The input impedance of the null detector is 1 megohm. Its sensitivity is 100 uy at 60 cycles and 200 
10 300 uv at 1000 eyeles. Its selectivity is 40 db against the second harmonic. Plug-in units tune the 


inplifier to any operating frequency desired between 20 and 2,000 cycles, with a continuous tuning | 
range of +5% for each unit. 


TYPE 707-A CATHODE-RAY NULL INDICATOR .... © «© © © «© « $195.00 
WRITE FOR BULLETIN 
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For Demonstrations of the 


PIEZO-ELECTRIC EFFECT 


81330 81335 80796 


ROCHELLE SALT CRYSTAL 


\ nearly flawless half crystal, of approximately 350 grams weight and measuring 10 by 8 
em, in which the piezo-electric effect is very large and reversible. By cementing tinfoil 
coatings on opposite faces, alternating emf can be applied to cause it to vibrate. With 
the output of a radio receiver, audible sounds are emitted. 

81330 Piezo-Electric Crystal 


Each $10.00 


PIEZO-ELECTRIC DEMONSTRATION 


Shows production of an electric charge by mechanical distortion of a Rochelle salt crystal 
and the converse effect. When crystal is tapped, the neon bulb flashes. A shock is 
received when terminals are touched with moistened fingers. When neon bulb is replaced 
by a connection to the output of a radio receiver and a paper cone attached to the upper 
end of the crystal mounting, audible sounds are produced. 

$1335 Piezo-Electric Demonstration Apparatus. . Kach $16.50 


-PIEZO-ELECTRIG TELEPHONE RECEIVER 


\ commercial application of the piezo-electric effect which is excellent for laboratory work 


with alternating currents. Can be used for long periods without fatigue and loss of 


sensitivity of the ear. Especially suited for balancing A.C. bridges. It is non-magnetic, 
has high impedance, high current sensitivity, draws negligible current, and vibrates with 
vreat amplitude in response to applied potentials. © Provided with head band and standard 
5-foot telephone cord with pin tips. 


807960 Piezo-Electric Telephone Receiver .... Each $3.50 
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Lakeview Station OFF Cambridge A Station 


LANCASTER PRESS, INC., LANCASTER, PA. 


] 
4 
| 
‘ 
\ 
\ 
\ 
| 


